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ABSTRACT 

An  attempt  is  made  to  determine  the  amplitude 
absorption  and  the  velocity  dispersion  of  a  high 
frequency  plane  sound  wave  propagating  through  a 
dilute  monoatomic  gas.  The  problem  is  formulated  in 
terms  of  the  Boltzmann  collision  equation  with  the 
collision  integral  replaced  by 

=  (f-?) 

^  coll.  t 

where  f  is  the  local  equilibrium  (displaced  Maxwellian) 
distribution  in  the  presence  of  the  sound  wave  and  t 
is  the  so-called  relaxation  time.  Once  the  distribution 
function  f  has  been  derived  from  the  Boltzmann 
equation,  the  linearized  conservation  equations  (number, 
momentum  and  energy)  are  calculated.  The  determinant 
of  these  homogeneous  equations  yields  the  dispersion 
relation,  relating  the  complex  wave  number  k  to  the 
sound  frequency  cd.  Expressions  for  the  amplitude 
absorption  coefficient  a  and  the  phase  velocity  V(gd) 
follow  immediately. 
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Two  special  cases  are  considered  for  t  , 
namely,  (i)  t  is  independent  of  the  instantaneous 
particle  velocity  v  ,  (ii)  t  =  $/v  where  the  mean 
free  path  &  does  not  depend  on  _v  .  For  the  former 
case,  it  is  shown  that,  for  cut  «  1,  the  results 
obtained  for  a  and  V  exhibit  the  same  qualitative 
behaviour  as  the  well-known  relations  of  Wang-Chang 
and  Uhlenbeck.  Our  velocity  dispersion  curve  fits  the 
experimental  data  better  than  theirs  but  our  absorption 
curve  displays  less  agreement.  For  cut  »  1,  it  is 
shown  that  V  becomes  frequency- independent  once  more, 
while  cl  increases  linearly  with  the  frequency. 
Furthermore,  cl  no  longer  depends  on  t  (which  is 
related  to  the  type  and  state  of  the  monoatomic  gas 
being  considered) .  The  latter  case  of  H  independent 
of  v  is  considered  in  an  attempt  to  improve  the  slope 
of  the  absorption  curve  for  low  frequencies  (cut  «  l) . 
It  is  found  that  there  is  some  slight  improvement  but 
not  enough  to  suggest  a  re-examination  of  the  high 
frequency  situation.  Some  suggestions  (such  as  the 
extension  of  the  calculations  to  polyatomic  gases)  are 
made  for  future  investigation. 
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Chapter  1 .  INTRODUCTION 


The  study  of  the  absorption  and  dispersion 
of  sound  in  monoatomic  gases  is  nearly  a  century  old,  yet 
one  still  lacks  a  fundamental  understanding  of  how  these 
phenomena  behave  at  very  high  frequencies. 

The  low  frequency  behaviour  of  the  absorption 
coefficient  a  and  the  phase  velocity  V  is  well  understood 
and  may  be  explicitly  stated,  but  as  G.  Uhlenbeck1  (1963) 
phrased  it,  "what  happens  for  a)  -*  00?"  .  Is  the  collective 
motion,  which  a  sound  wave  represents,  the  same  at  very 
high  frequencies  as  at  lower  frequencies? 

p 

Numerous  investigators,  such  as  H.  Primakoff 
(19^2),  have  suggested  that  for  sufficiently  high  frequencies 
no  sound  wave,  "in  the  ordinary  sense",  will  be  propagated. 

C.  Kittel^( 19^6) ,  in  fact,  states  that  "the  mean  free  path 
is  effectively  a  lower  limit  to  the  wavelength  of  sound 
waves  in  a  gas" . 

Others  have  speculated  that  a  will  go  through  a 
maximum  and  then  decrease,  while  V  approaches  a  constant 
value . 


But  there  are  no  solid  reasons  for  either 
belief.  That  is  to  say,  there  are  no  consistent  derivations 
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which  yield  either  of  the  above  results. 

Various  theoreticians  have  solved  the  problem 
at  low  and  intermediate  frequencies.  Their  predictions 
agree  very  well  with  the  results  of  experiment  at  low 
frequencies,  but  diverge  from  the  experimental  points  as 
the  frequency  increases.  At  very  high  frequencies  their 
mathematical  apparatus  breaks  down  completely,  because  of 
the  inherent  lack  of  generality  in  their  assumptions. 

One  difficulty  that  makes  the  theoretician's 
task  in  this  direction  quite  arduous  is  the  lack  of 
experimental  data  at  very  high  frequencies.  Thus,  one 
cannot  say  that  his  proposed  theory  agrees  or  disagrees 
with  experiment,  but  merely  that  the  results  derived  are 
completely  consistent  with  the  basic  physical  or  mathematical 
assumptions  originally  made. 

The  primary  purpose  of  this  thesis,  then,  shall 
be  to  formulate  a  theoretical  model  which  will  consistently 
predict  the  high  frequency  sound  behaviour  in  monoatomic 
gases . 

As  a  very  brief  side  issue,  it  will  be  pointed 
out  how  this  model  may  be  extended  to  explain  the 
behaviour  of  sound  in  polyatomic  gases  (such  as  oxygen) 
where  there  exist  discrepancies  between  theory  and 
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experiment.  The  machinery  for  carrying  out  this  latter 
programme  will  be  outlined,  but  no  detailed  calculations 
will  be  given  at  this  time. 


4. 


Chapter  2.  EXPERIMENTAL  SITUATION 

2 . 1  Historical  Background 

Before  delving  into  the  theoretical  aspects  of 
the  problem,  let  us  briefly  outline  the  experimental 
work  that  has  been  carried  out  in  the  past,  and  discuss 
the  difficulties  standing  in  the  way  of  further  experiment¬ 
ation  at  higher  frequencies . 

N.  Neklepawjev^  (1911)  was  the  first  to  make 
absorption  measurements  in  air  but  his  results  were 
affected  by  the  carbon  dioxide  content.  The  same  may  be 
said  of  the  experiment  of  D.  L.  Rich  and  W.  H.  Pielemeier^ 
(1925).  In  1925.>  G.  Pierce^  was  able  to  detect  the  change 
of  phase  velocity  with  frequency  for  the  first  time.  He 
experimented  with  air  and  carbon  dioxide . 

A.  van  Itterbeek  and  P.  Mariens^  (1940)  measured 
the  absorption  of  sound  in  Helium  and  other  rare  gases. 

The  observed  absorption  was  four  times  greater  than  the 
value  which  was  predicted  on  the  basis  of  classical 
theory . 


Finally,  it  was  shown  by  van  Itterbeek  and 
Verhaegen^  (1951)  and  by  M.  Greenspan^* 10 J 11  (1950,  1953* 
1956)  that  van  Itterbeek* s  earlier  results  were  adversely 
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affected  by  impurities  and  other  experimental  difficulties 
Their  results,  which  did  agree  with  classically  predicted 
values,  were  more  recently  checked  by  the  Russian 
investigators  I.  I.  Perepechko  and  V.  P.  Yakovlev  (i960) 

We  shall  very  briefly  describe  Greenspan's 
experimental  setup  as  it  illustrates  several  of  the 
difficulties  in  one's  quest  to  go  to  higher  frequencies. 

A  block  diagram  of  Greenspan's  setup  is  indicated  in 
Fig.  1.  The  gas  under  study  is  contained  in  a  cylindrical 
cavity.  Plane  sound  waves  are  emitted  by  a  quartz  crystal 
at  one  end  and  detected  by  another  (movable)  crystal  at 
the  other  end.  This  so-called  double  crystal  interfero¬ 
meter  is  quite  often  used  in  high  frequency  absorption 
work.  In  1950  Greenspan  used  crystals  with  a  fundamental 
frequency  of  1  Mc/s.,  but  in  his  1956  paper  he  worked 
with  11  Mc/s.  crystals.  Rather  than  vary  the  frequency 
oo  in  the  interferometer  he  varied  the  gas  pressure  p. 

For  not-too-low  pressures  (otherwise,  he  would  have  the 
well-known  Knudsen  case),  this  procedure  is  equivalent  to 
varying  the  frequency  as  the  pertinent  quantities  a,  and 
V  depend  only  on  the  ratio  a)  to  p.  That  this  was  true 
was  obvious  from  the  fact  that  both  his  1  Mc/s.  and  11 
Mc/s.  measurements  yielded  the  same  results. 
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Greenspan's  Experimental  Set-up 
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A  detailed  explanation  and  diagram  of  the  physical 
construction  of  the  interferometer  may  be  found  in 
Greenspan’s  1953  paper.  The  detailed  theory  of  this 
double-crystal  interferometer  was  first  described  by  W. 

Fry13  (1949). 

In  Fig.  1,  the  power  oscillator  fed  a  signal 
into  the  11  Mc/s.  source  resulting  in  a  sound  wave  passing 
through  the  gas  to  the  11  Mc/s.  receiver  (which  could  be 
moved  back  and  forth) .  The  received  signal  was  fed  Into 
a  broadband  preamplifier  and  then  mixed  with  a  10. 56  Mc/s. 
signal  from  the  first  local  oscillator.  An  Intermediate 
frequency  of  440  Kc/s.  resulted  and  this  was  passed 
through  a  series  of  wide-band  amplifiers  and  (low-band) 
matching  filters.  The  resulting  overall  band  width  of 
23  Kc/s.  was  determined  by  the  crystal  filter.  The 
intermediate  frequency  signal  was  then  rectified  by  a 
germanium  diode  and  fed  into  a  logarithmic  recorder 
which  gave  a  direct  measure  of  twice  the  absorption  coef¬ 
ficient  a.  Electrical  "cross-talk"  (or  interference) 
was  neutralized  by  taking  the  signal  from  the  power 
oscillator,  mixing  it,  then  feeding  it  through  the 
attenuator  and  phase  shifter  into  the  input  of  the  first 
intermediate  frequency  amplifier. 
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The  distance  between  the  emitter  and  the  receiver 
was  made  sufficiently  large  so  that  the  large  absorption 
of  the  sound  by  the  gas  prevented  the  formation  of  standing 
waves.  The  apparatus  was  not  run  in  the  interference 
region . 


The  phase  change  was  determined  by  comparing  the 
signals  from  the  first  reference  mixer  and  the  first 
signal  mixer  by  means  of  a  commercial  phase  meter.  The 
details  of  such  a  phase  meter  may  be  found  in  a  paper  by 
E.  Meyer  and  G.  Sessler1^  (1957) • 

Greenspan's  results,  for  the  monoatomic  gases 
He,  Ne,  A,  Kr,  and  Xe,  are  indicated  in  Fig.  2.  The  more 
recent  results  of  Meyer  and  Sessler  are  presented  in  Fig. 

3,  but  a  word  of  caution  should  be  voiced  about  their 
seemingly  high-frequency  results.  Unlike  Greenspan, 
they  took  the  emitter-receiver  distance  so  short,  and  the 
gas  pressure  so  low,  that  the  mean  free  path  of  the  gas 
particles  was  greater  than  the  cavity  dimensions.  Thus, 
most  of  the  sound  damping  was  due  to  collisions  of  the 
particles  with  the  cavity  walls,  rather  than  with  each 
other.  That  is  to  say,  they  really  considered  the  Knudsen 
case.  However,  their  low  frequency  data  is  certainly  valid 
and  agrees  with  that  of  Greenspan. 
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2 . 2  Difficulties  and  Future  Work 

This  last  point  in  Sec.  2.1  indicates  one  of  the 
difficulties  which  one  encounters  in  trying  to  go  to 
higher  frequencies .  One  must  overcome  this  by  using 
crystals  with  a  much  higher  fundamental  frequency  than 
those  presently  being  used  in  gaseous  work.  Consequently, 
if  one  is  to  keep  the  emitter-receiver  distance  at  least 
a  few  mean  free  paths  in  length,,  one  must  either  develop 
more  sensitive  detecting  equipment  or  else  use  waves  of 
higher  amplitudes  in  order  to  overcome  the  effect  of 
large  sound  absorption  in  gases.  However,  one  cannot 
make  the  amplitude  too  large  or  else  the  waves  will  no 
longer  be  sinusoidal  but  saw-toothed  in  appearance,  ie . , 
they  will  no  longer  be  sound  waves,  but  rather  nearly 
shock  waves . 

No  ultrasonic  measurements  whatsoever  have 
been  made  above  a  few  thousand  Mc/s.  The  reason  for 
this  may  be  explained  on  the  basis  of  the  following 
example:  Suppose  that  we  want  a  crystal  with  a  funda¬ 

mental  frequency  of,  say,  3000  Mc/s.  As  pointed  out  by 
Kittel^,  this  crystal  must  be  only  10,000  A  thick.  This 
is  much  too  thin  for  one  to  prepare  to  the  optical 

tolerances  required  (3000  Mc/s.  corresponds  to  a  wave- 

-4 

length  in  quartz  of  about  2  x  10  cm.  which  is  only 
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3  or  4  times  the  wavelength  of  visible  light),  so  that  one 
must  use  a  thicker  crystal  driven  at  a  higher  harmonic. 
However,  very  little  energy  will  be  present  in  the  higher 
harmonics  and  if  the  crystal  is  not  uniformly  thick  to 
a  high  degree  of  accuracy  there  will  be  interference 
between  the  different  harmonics.  Furthermore,  quartz  has 
a  high  absorption  coefficient  at  high  frequencies  so  that 
a  great  deal  of  the  energy  is  lost  in  the  crystal,  in 
the  form  of  heat,  and  very  little  energy  will  be  trans¬ 
mitted  to  the  gas. 

To  ensure  plane  waves,  the  transducer  and  receiver 
must  be  very  carefully  aligned. 

Impurities  can  play  havoc  with  the  results,  but 
with  our  modern  methods  of  purifying  gases  this  can 
generally  be  ruled  out. 

Other  practical  difficulties  (such  as  the  need 
for  extreme  electrical  stability)  exist,  but  we  will  not 
mention  them  here.  The  experimentalist  is  referred  to 
Greenspan1  s  papers^ 3  for  further  details. 

Roughly  the  same  sort  of  techniques  may  be  used 
in  investigating  the  behaviour  of  ultrasound  (which  refers 
to  sound  with  a  frequency  >  20,000  c/s.)  in  diatomic 
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and  polyatomic  gases.  For  more  explicit  information,  the 
reader  is  referred  to  the  text  of  K.  F.  Her zf eld  and 
T.  A.  Litovitz  J  (1959)  which  contains  a  summary  of  work 
in  this  direction. 
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Chapter  3.  THEORETICAL  BACKGROUND 

3 . 1  Physical  Ideas  Behind  Absorption  and  Dispersion 

The  physical  Idea  behind  sound  absorption  is 
that  sound  energy  is  absorbed  by  a  medium  because  the 
directed  motion  of  the  sound  wave  is  irreversibly  converted 
into  random  thermal  motion  of  the  constituent  particles. 
This  absorption  depends  on  quantities  characteristic  of 
this  irreversible  process,  namely,  the  static  shear 
viscosity  r\  and  the  static  thermal  conductivity  x*  The 
former  refers  to  the  conversion  of  sound  energy  into  random 
thermal  energy  by  the  process  of  internal  friction  or, 
more  explicitly,  by  the  molecular  transport  of  momentum 
in  the  (negative)  direction  of  the  velocity  gradient  in 
the  medium.  Sound  absorption  due  to  heat  conduction  may 
be  explained  as  follows.  Heat  conduction  causes  energy 
to  flow  from  the  density  crests  (which  are  hotter  than 
the  average)  to  the  density  troughs  (which  are  cooler 
than  the  average),  thus  equalizing  the  temperature.  The 
regions  compressed  by  the  sound  wave  will,  upon  re¬ 
expanding,  return  less  energy  to  the  wave  than  was 
spent  in  compressing  them.  The  same  is  true  of  the 
recompression  of  rarefaction  regions.  The  result  is 
that  the  sound  wave  loses  energy. 
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Qualitatively,  what  do  we  expect  for  the  sound 
dispersion  ( ie . ,  the  variation  of  phase  velocity  with 
frequency)?  As  pointed  out  by  K.  P.  Herzfeld  and  P.  0. 

Rice1^  (1928),  the  effect  of  viscosity,  or  internal  friction, 
is  to  retard  the  movement  of  a  given  wave  form.  Thus,  for 
a  given  frequency,  the  wavelength  A  will  be  larger  with 
friction  than  without.  The  net  result  is  that  the  phase 
velocity  is  increased  by  the  internal  friction  and,  there¬ 
fore,  increases  with  frequency. 

Heat  conduction  causes  the  phase  velocity  to 
decrease  with  frequency.  With  increasing  frequency,  the 
time  allowed  for  heat  conduction  decreases.  However,  the 
amount  of  heat  conducted  at  any  moment  is  proportional  to 

r)2ip  1 

— -  which  is  clearly  proportional  to  .  Therefore,  the 

dx2  A2 

heat  conducted  during  one  period  increases  with  frequency 

due  to  the  increasing  steepness  of  the  temperature  gradient. 

Thus,  the  gas  tends  to  become  more  isothermal  and  less 

adiabatic.  This  suggests  that  V  will  tend  to  decrease 

with  a)  from  Laplace's  (adiabatic)  result  1//2  to 

Newton's  (isothermal)  result  (pf)1^2  • 

Quantitatively,  it  turns  out  that  the  net  effect  of 
these  two  processes  is  for  the  phase  velocity  to  increase 
with  frequency. 

With  this  qualitative  background  behind  us,  let 
us  look  at  previous  attempts  to  quantitatively  solve  the 
problem  of  sound  absorption  in  gases. 
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3 . 2  Navier  -  Stokes  Theory  of  "Classical  Absorption" 

Long  before  any  experimental  work  was  carried 
out,  G.  Stokes^  (1845)  and  G.  Kirchoff1^  (1868) 
calculated  the  effect  of  r|  and  x  on  sound  absorption 
and  dispersion. 

They  considered  a  small-amplitude  plane  sound 
wave  propagating  through  a  gas,  say,  in  the  z-direction, 

viz.,  ei(u)t-kz) 

Writing  k  =  k-[  -  ik2  2  ^  -  ia  ;  a  >  0  (3.1) 

the  plane  wave  may  be  written  as  e~az  •  eicD^  “ 

Thus,  for  a  >  0,  one  has  a  plane  wave  being  progressively 
damped  as  it  moves  along  the  z-axis.  The  quantity  a  is 
called  the  amplitude  absorption  coefficient  and  it  has  the 
physical  significance  that  the  amplitude  of  the  wave  is 
l/eth  of  its  initial  value  (at  z=0)  after  the  wave  has 
travelled  a  distance  l/a  cm.  It  is  obvious  that  the 
absorption  coefficient  for  the  sound  intensity  is  2a. 

I 

The  term  dispersion  refers  to  the  fact  that  the  phase 
velocity  V  may  be  a  function  of  frequency. 

Now  let  us  briefly  outline  the  Stokes-Kirchof f 
derivation  of  sound  absorption  and  dispersion,  which  is 
of  more  than  just  academic  interest  to  us  here. 
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The  general  equation  of  motion  of  the  fluid  Is 


5u.  3E]  . 

p  ar  =  srf 


(3.2) 


What  Is  the  form  of  E^j  (or  of  =  -E^j,  the  pressure 
tensor),  the  stress  tensor?  By  appealing  to  everyday 
observations  such  as:  (l)  a  fluid  at  rest,  or  in  uniform 
motion,  cannot  sustain  shearing  forces,  (2)  for  the  whole 
fluid  in  uniform  rotation,  there  are  no  shearing  stresses; 
and  assuming  that  the  velocity  gradients  involved  are 
small,  one  can  construct  E.  .  or  -P.  .  in  the  form 


-pij  -  E±j  -  -  p5ij  +  2Tieij  +  t div  b  qj 

(3.3) 


where  |  is  called  the  coefficient  of  volume  viscosity 
and 
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1  (Sui 

2  '  c*x  . 


+ 


du 

"<3x 


i) 


Then  E^  +  E^  +  E^  =  -  3p  +  2r)  div  u  +  3£  div  u. 

Since  one  would  not  normally  expect  shear  forces  to  contribute 
to  the  compression,  he  chooses  E-^  +  E ^  +  E^  =  -3p 


ie.. 


2t|  +3^  =  0 


(3.4) 


This  is  referred  to  as  Stokes1  assumption.  It  may  be  shown 
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to  be  identically  true  for  monoatomic  gases  by  using  the 
well-known  kinetic  theory.  However,  it  is  not  generally 
true  for  polyatomic  gases.  This  problem  has  been 
discussed  in  detail  by  L.  Tisza1^  (19^2). 

Then,  from  (3.2),  (3.3)*  and  (3.^)  the  equation 
of  motion  in  one  dimension  is 


du 

P  "St 


n 


(3.5) 


The  continuity  equation  is 


M  +  it  (pu)  =  0  ^3*6) 

One  then  linearizes  these  by  making  the  acoustic  approxi¬ 
mation,  ie  .  ,  one  writes 

p  =  pQ(l  +  A)  |  A  |  «  1 

P  =  PQ  +  A'  |  A'  |  «Eb  (3.7) 

For  example,  for  ordinary  speech  A  ~  10*  .  Then  (3.5) 
and  (3.6)  become 

8A/8t  +  8u/8z  =  0  (3*8) 

pQ(du/dt)  =  -  3a'/3z  +  -g  r) 


(3.9) 
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Using  the  second  Tds  relation  of  thermodynamics ,  one  can 
obtain  the  linearized  heat  flow  equation  (per  unit  mass) 


dT* 
Po  6t 


-  t  e 

o  at 


82T» 

8z^ 


(3.10) 


where  0  =  ^)p 

expansion  and 


is  the  usual  coefficient  of  thermal 


T  =  Tq  +  T«  ;  |  T*  |  «  T  (3.H) 

Since  these  three  equations  contain  four  unknowns,  one 
needs  another  equation.  By  applying  the  chain  law  of 
differentiation  to  the  volume  ^(p,T)  and  then  linearizing, 
one  has 


A 

where  KT 

compressibility . 


=  KtA’  -  0T '  (3.12) 

1  c)V 

=  ---  (-^p)T  is  the  isothermal 


If  one 
the  fluid  in  the 
proportional  to 
equations  (3.8), 


considers  a  plane  wave  propagating  through 
z-direction,  then  A,  u.  A’,  T1  are  all 
ei(oyt-kz)  .  Substituting  this  into 

(3.9)j  (3-10)  and  (3.12)  yields  four 


homogeneous  equations  in  four  unknowns .  For  there  to  be 
a  non-trivial  solution,  the  determinant  of  the  coefficients 
of  these  unknowns  must  be  zero.  From  experiment,  one 
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knows  that  —  «  1  and  ri  and  y  are  small.  Then, 
co  '  A 

on  separation  of  real  and  imaginary  parts  of  k  In  the 
de terminantal  equation,  one  obtains 

“  -  !  ^  +  1  (3-13) 

and  v  =  VQ  =  hHo)1/2  (3.14) 

for  an  ideal  gas.  If  one  uses  the  result  for  monoatomic 
gases  (obtainable  from  kinetic  theory)  that 


15  _ 
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then  one  can  estimate  the  relative  importance  of  thermal 
conduction  to  that  of  viscosity  in  sound  absorption.  In 
fact, 

ath. cond ,/avis .  ”  °-75 


To  this  approximation  there  is  no  sound  dispersion.  If 
one  had  solved  the  so-called  hydrodynamical  equations 
(3.8),  (3.9)*  (3.10),  and  (3-12)  to  the  second  order  in 
T)  and  x >  he  would  have  found  that 
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Po  v' 


p 


b^i(7-3Y)] 

4 


(3-15) 


ie . ,  to  the  next  order  there  Is  sound  dispersion,  but  of 
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course  It  is  quite  small  since  T)  and  x  are  small.  The 
complete  solution  of  equations  (3-8),  (3-9) ^  (3-10),  and 
(3-12)  (to  any  order  of  r\  and  x)  is  called  the 
Navier-Stoke ’ s  solution  and  the  expression  for  a  in 
(3*13) j  the  "classical"  absorption.  Although  one  can 
solve  (in  principle)  the  above  equations  to  any  order  in 
r)  and  X*  the  procedure  is  not  really  consistent.  To 
keep  higher  orders,  one  should  have  included  higher 
order  terms  in  the  stress  tensor  and  in  the  heat  flux 


vector  (x 


d2T 

dx2 


term)  .  This  was  not  done,  so  that  the 


above  results  are  consistent  only  to  the  first  order  in 


aV 


«  1 


q  and  x-  From  (3.13)  it  can  be  seen  that 
amounts  to  ^  «  1  (since  p  =  PQV^)  .  From  the  kinetic 
theory  (apart  from  numerical  factors),  q  ~  p&c,  where  c 
is  of  the  same  magnitude  as  V,  so  that,  essentially, 

1/ A  «  1  where  A  is  the  sound  wavelength  and  &  the 
mean  free  path  (which  is  of  the  order  of  10~^  cm.  in 
gases  at  ordinary  pressures) .  Thus,  the  expressions  for 
the  "classical"  absorption  and  dispersion  are  only 
valid  for  low  frequencies  or  very  long  wavelengths. 


For  monoatomic  gases  at  low  frequencies,  the 
Stokes-Kirchoff ,  or  Navier-Stokes,  theory  agrees  with  the 
experimental  results  (see  Fig.  2)  but  as  the  frequency 
increases  ( ie . ,  H/h  becomes  larger)  the  theoretical  curves 
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diverge  sharply  from  the  experimental  points.  In  Pig. 

4  is  the  curve  obtained  by  solving  the  hydrodynamical 
equations  numerically  for  any  q  and  x-  However ,  these 
Navier-Stokes 1  results  are  for  Maxwellian  molecules  (l/r^ 
interaction  law)  and  were  arrived  at  in  a  different  way 
than  has  been  previously  indicated.  This  will  be  discussed 
shortly . 

3 . 3  Thermal  Relaxation 

For  polyatomic  gases  it  was  found  that,  apart 

from  the  above  high  frequency  limitations,  theory  did  not 

agree  with  experiment.  This  was  not  explained  until  1928 

16 

when  K.  Herzfeld  and  F.  0.  Rice  introduced  the 
phenomenological  concept  of  "thermal  relaxation".  They 
argued  that  the  additional  energy  loss  was  due  to  the 
slow  energy  exchange  between  the  translational  degrees 
of  freedom  and  the  internal  (rotational  and  vibrational) 
degrees  of  freedom.  The  absorption  due  to  "thermal 
relaxation"  may  be  one  hundred  times  greater  than  that 
due  to  the  viscosity  or  heat  conduction.  If  the  exchange 
is  very  fast,  there  would  be  no  sound  wave  at  all.  As 
the  frequency  increases,  the  sound  absorption  increases 
to  a  maximum  and  then  decreases.  At  very  low  frequencies, 
the  external  and  internal  degrees  of  freedom  will  be  in 
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equilibrium  and  the  phase  velocity  therefore  a  constant 
value.  In  fact  it  is  given  by  the  usual  expression 
(^p)1/2  where  7  =  7/5  if  only  rotation  is  considered 
and  7  =  9/7  if  both  rotation  and  vibration  are  included. 
As  the  frequency  increases,  the  internal  degrees  of 
freedom  follow  the  external  degrees  less  and  less.  The 
effective  specific  heat  of  the  system  decreases,  thus 
causing  the  phase  velocity  to  increase.  As  the  frequency 
becomes  very  high  there  is  no  energy  exchange  so  that 
the  energy  remains  with  the  translational  degrees  of 
freedom.  This  means  that  the  phase  velocity  levels  off 
to  a  new  constant  value. 

Since  we  are  not  going  to  give  any  detailed 
calculations  of  sound  absorption  and  dispersion  in 
polyatomic  gases  in  this  thesis,  we  shall  say  nothing 
further  on  the  subject.  For  those  who  are  interested, 
they  should  read  the  excellent  account  in  K.  F.  Herzf eld's 
text ^  (1959)  "Absorption  and  Dispersion  of  Ultrasonic 
Waves" . 

3 • 4  Chapman  -  Enskog  Development 

Let  us  return  to  the  fundamental  problem  of 
how  to  treat  sound  absorption  and  dispersion  in  monoatomic 
gases  as  the  frequency  increases  (ie.,  for  i/ A  approaching 
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unity) . 


The  problem  of  expanding  the  pressure  tensor  and 
heat  flux  vector  to  higher  orders  of  ft /A  was  independently 
attacked  by  S.  Chapman  ^0,21  ^^6,  1917)  and  D.  Enskog22,23 
( 1917,  1921).  Although  their  methods  were  slightly  differ¬ 
ent,  their  results  were  identical.  Chapman's  method  was 

oil 

inspired  by  a  paper  by  Maxwell  (1879)  while  Enskog's 

2E 

results  were  the  culmination  of  earlier  work  by  Lorentz  ^ 

P6 

(1905)  and  Hilbert  (1912).  The  Chapman-Enskog  results 

27  28 

were  derived  in  a  more  sophisticated  manner  by  D.  Burnett  ' * 

(1935)  who  was  the  first  to  use  Sonine  polynomials  in  this 

connection.  He  was  the  first  to  derive  P.  .  in  its  complete 

-  J 

form  to  the  second  order  in  spatial  derivatives  of  the 
macroscopic  variables  (temperature,  etc.)  and  to  the 
second  power  of  the  derivatives.  Chapman  had  been  the 
first  to  accomplish  this  for  the  heat  flux  vector, 
although  his  results  contained  numerical  errors.  Burnett's 
work  laid  the  foundation  for  most  of  the  later  papers  on 
sound  absorption  and  dispersion  in  gases. 

The  Chapman-Enskog  approach  may  be  briefly 
described  as  follows:  For  a  classical,  dilute,  monoatomic 
gas,  one  can  write  down  the  well-known  nonlinear  Boltzmann 
collision  equation  describing  the  time  evolution  of  the 
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gaseous  distribution  function  f(r,v,t)  at  a  point  r 
due  to  two  effects:  (l)  the  streaming  of  particles 
through  the  volume  dr  at  r,  (2)  the  net  gain  of  particles 
in  dr  due  to  collisions.  Viz.., 


-gf  +  v  •  V^f  +  a  •  Vvf  =  /dVj./dflgltg,  9)  [f 'f  x ' -f^]  =  J(f,f1) 

say 

(3.16) 
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acceleration  of  gas  particles  due  to  external 
f orces , 

element  of  solid  angle 

relative  speed  of  two  colliding  particles 
before  or  after  collision 
differential  cross-section. 


The  collision  equation  (3.16)  will  be  discussed  in  detail 
in  the  next  section.  If  one  multiplies  equation  (3.16) 

p 

by  1,  v,  v  ,  and  then  integrates  over  all  of  velocity 
space,  he  obtains  the  continuity,  momentum  conservation 
and  energy  conservation  equations  respectively. 
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(3.17) 


which  are  expressed  in  terms  of  macroscopic  variables, 
provided  that  we  use  the  following  definitions 


p  =  mN  ;  N  =  f fdv  ;  =  <v. >  =  Jv^fdv 


1  _  tt2 


Q  =  2  PU  ,  the  thermal  energy  density 


where  Uh  =  v^  -  <v^> 


\  .  =  p<LhUj>  ,  the  pressure  tensor 


and  D 


1  2^ 

2  p  <U^U  >  ,  the  heat  flux  vector 

l  ,^Uj_  o u  -* . 

"2  ~3x~  ,  the  rate  of  strain  tensor 


Quantities  such  as  .  may  be  written  out  in  full,  once 
the  form  of  f  has  been  determined  in  detail  from  the 
Boltzmann  collision  equation  (3.16). 


The  Chapman- Enskog  method  essentially  consists 
of  solving  equation  (3.16)  for  f  by  perturbation 
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techniques.  One  writes 


(0) 


(1)  4-  ,2.(2) 


f  =  f^'  +  efv^'  +  e  f v  '  +  ... 


(3.18) 


where  the  zeroth  order  distribution  function  f 
just  the  local  Maxwellian  distribution. 


(0) 


is 


f 


m  ,  p 

(0)  _  3/2  -  -2KT  ll  -  Hi 

in27FKT; 


(3.19) 


where  u  is  the  macroscopic  velocity  imposed  on  the 
particles  due  to  the  sound  wave.  If  ( 3 • 18)  is  a 
perturbation  expansion ,  what  is  the  expansion  parameter  e? 
It  turns  out  that  (3.18)  may  formally  be  written  as 


f  =  f(°)[l  +  A1(v)(£V)  +  A2(v)(£V)2  +  . . .]  (3.20) 

where  the  gradient  acts  on  the  macroscopic  variables 
u(r, t ) ,  T(r,  t ) ,  N(r, t)  and  A^,  A^,  .  ..,  which  depend  on 

the  molecular  force  law,  are  determined  from  the  collision 
equation.  One  can  see  that  e  is  the  relative  change 
of  the  macroscopic  variables  over  a  mean  free  path.  In 
(3-20),  one  will,  in  general,  have  second  derivatives 
of  u  and  T  and  also  squares  of  the  first  derivatives. 

However,  a  sound  wave,  by  definition,  does  not  contain 
second  or  higher  powers  of  the  derivatives  of  u,  etc., 
as  the  derivatives  are  taken  to  be  small.  Not  only 
does  this  correspond  to  physical  reality,  but  it  introduces 
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mathematical  simplicity.  For  a  sound  wave  then,  e  ^  Z/A 
so  that  the  Chapman-Enskog  approach  will  break  down  when 
applied  to  high  frequency  sound  propagation.  How  does 
one  actually  go  about  determining  the  forms  of  f^^, 
f^2^,  ...  ?  Instead  of  equation  (3.16),  consider  the 

modified  Boltzmann  equation 

df/dt  +  v  •  Vrf  +  a  •  Vvf  =  j  J(f,f1)  (3.21) 

From  equation  (3.18),  we  have 

J(f,f1)  =  J(f(°)  +  ef^1)  +  e2f(2)  +  ...,  f.^0)  + 

ef  1^ ^  +  e2f  ^2)  +  ...  ) 

=  J(f  (°)^fi(°) )  +  €[  j(f  (1),f1^0))+  J(f  )  ] 

+  e2[  J(f  (^f^0) )  +  J(f^1^,f1^1^) 

+  J(f  ^°^f1^2^ )  ]  +  .... 

(3.22) 

Then,  since  e  is  arbitrary,  we  may  substitute  (3-18) 
and  (3.22)  into  (3.21)  and  equate  equal  powers  of  e. 

This  procedure  yields 

J(f  ,1-^°) )  =  0 

df(°)/8t  +  v  •  Vrf^0^  +  a  •  Vyf^0)  =  J(f^1^,f1^°h  +  J(f^°^f1^1^) 

df^/dt  +  v  *  Vpf^1)  +  a  *  Vyf^1^  =  J(f  ^  ,f1^0^ )  +  J(f  ^  ,f1^1') ) 

+  J(f  ..f.^2) ) 


(3.23) 
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and  so  on.  Consider  the  first  of  these  coupled  equations 


implies  that 


=  0 


f(0)f  (0) 
1 


f(0) 'f  (0) 1 


which  may  be  satisfied  only  if 

f(°)  -  Nf  m  )  3/2  c-fm/2KTH  v-u  |2 
1  “  wi2ttKTj  e 


which  is  consistent  with  what  was  previously  stated.  If 

this  approximation  to  f  is  substituted  into  the 

transport  equations  (3.17)*  one  just  obtains  the  Euler 

hydrodynamical  equations.  The  coefficient  of 

e°  =  1  yields  hydrodynamic  equations  independent  of  e . 

Similar ly,  the  second  equation  in  (3.23)  (corresponding 

to  setting  the  coefficient  of  e  equal  to  zero)  yields 

the  Navier-Stokes  equations  which  contain  terms  linear 

2 

in  e.  The  equation  corresponding  to  e  yields  the 
Burnett  hydrodynamical  equations,  and  so  on. 

Since  f(^)  has  been  determined,  it  is  now 

possible  to  solve  the  second  equation  of  (3.23)  for  f^) 

Once  is  explicitly  determined,  one  may  solve  the 

(o) 

third  equation  for  fv  ’  and  so  on.  In  this  manner,  it 
is  possible  to  obtain  the  perturbation  expansion  for  f 
to  as  many  terms  as  is  desired. 
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Then,  knowing 
q^  to  any  order  of  e . 
order  (Burnett)  in  e, 
situation. 


f,  one  may  determine  and 

It  turns  out  that,  to  the  second 
one  has  for  the  one-dimensional 
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(3.24) 


The  explicit  values  of  the  0’s  and  rj  and  x  depend 
on  the  particular  force  law  chosen. 


Of  course,  the  general  Burnett  result  contains 
second  powers  in  the  derivatives  ( ie . ,  (-^)2  )*  but  we 
have  not  written  these  down  because,  as  stated  previously, 
we  are  only  interested  in  terms  linear  in  the  derivatives 
for  the  case  of  sound  waves. 


In  addition  to  the  papers  already  cited,  the 
reader  is  referred  to  the  account  of  S.  Chapman  and 
T.  G.  Cowling22  (1952)  for  the  details  of  the  fore- 
mentioned  results. 


3 . 5  Application  to  Sound  Propagation 

If  one  then  substitutes  the  explicit  expressions 
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obtained  in  the  previous  section  for  the  pressure  tensor 
and  the  heat  flux  vector  into  the  hydrodynamical  equations, 
one  can  obtain,  carrying  out  the  same  type  of  calculation 
as  previously  in  obtaining  the  Navier-Stokes  solution, 
the  Burnett  solution  to  the  problem  of  sound  absorption 
in  monoatomic  gases.  According  to  C.  S.  Wang-Chang  and 
G.  Uhlenbeck^-*  ^  (19^8),  the  Burnett  solution  for 
Maxwellian  molecules  is 


0 


(3.25) 


where  R  =  p  /r)0). 


o 

Similarly,  the  super-Burnett  solution  (valid  to  eJ)  is 


I:8[157  x  25 


1  12 

'  r5  '  —  • 


24  '  r3 
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+  3  .  i]  +  +  £3  i]  +  x 

2  R  6  R 


0 


(3.26) 


For  R  »  1,  ie . ,  for  Z /A  «  1,  (3-26)  yields 


(3.27) 


V 


VQ  [1  +  1.07(^)2] 

XT 


(3.28) 
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where  V  =  (-^p- ) as  usual, 

o  3m 

If  one  formally  solves  (3-26)  for  R  «  1,  he 
obtains  (to  the  lowest  order  in  R) 

— 2  =  o.58(-£.)V2 

—  =  °-58(^)1/2  (3.29) 

v  ^ 

Although  one  can  formally  obtain  these  latter  results, 
they  are  not  really  consistent  results.  That  is  to  say, 
they  are  not  really  consistent  with  the  assumptions 
employed  in  deriving  them  which  considered  ^  to  be 
small.  For  ^  large  the  Chapman-Enskog  method  will 
not  lead  to  good  convergence . 

Although  it  is  really  not  legitimate  to  do  so 
(for  the  above  reasons)  for  ^  >  1,  one  can  solve  the 
Navier-Stokes,  Burnett,  and  super-Burnett  solutions 
exactly  for  any  R  by  numerical  methods.  These  results 
are  presented  in  Fig.  4.  Of  course,  the  behaviour  of  the 
curves  for  R  »  1  and  R  «  1  is  just  that  predicted  by 
(3.27,  (3.28)  and  (3.29) . 

It  should  be  pointed  out  that  Wang- Chang  and 
Uhlenbeck  were  not  the  first  to  apply  the  Burnett  results 
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on 

to  sound  propagation  in  gases.  K.  F.  HerzfeldJ  (1935)* 

H.  Primakoff^  (1942),  and  H.  Tsien  and  R.  Schamberg^ 

(1946)  made  various  calculations  based  on  Burnett’s  work, 
but  their  results  were  all  affected  by  numerical  errors. 

As  far  as  is  known,  Wang- Chang  and  Uhlenbeck  were  the 
first  to  obtain  reasonable  error-free  results.  Their 
work  represents  the  final  culmination  of  the  Chapman- 
Enskog  procedure . 

3 . 6  13  Moment  Solution  of  the  Boltzmann  Equation 

In  the  last  decade,  various  other  approaches  to 

the  problem  have  been  suggested.  One  of  these  is  due  to 

34 

H.  GradJ  (1949)  who  solved  the  Boltzmann  equation  in  a 
different  manner  than  Enskog.  One  may  look  at  Enskog's 
approach  in  a  different  way.  Clearly,  a  knowledge  of 
the  macroscopic  variables  (N,  u,  T)  is  equivalent  to 
knowing  the  first  five  velocity  moments  of  the  distribution 

p 

function  ( ie . ,  /  1  •  fdv,  /  Vj_fdv,  J  v  fdv).  The  Boltzmann 
equation  then  completely  describes  the  future  of  the  gas 
in  terms  of  these  five  moments,  at  a  particular  time, 
at  each  point  in  space.  At  a  given  instant,  many  functions 
f(r,  v,  t)  have  the  same  first  five  moments  and  each 
leads  to  a  different  distribution  and  therefore  (in  general) 
a  different  set  of  first  five  moments  at  a  later  time. 
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Thus,  the  Enskog  procedure  does  not  yield  a  general 
solution  but  actually  restricts  one  to  a  particular 
class  of  distribution  functions. 

Grad  has  introduced  a  method  of  solving  the 

Boltzmann  equation  which  justifies  Enskog' s  first 

approximation  and  overcomes  the  above  difficulties. 

His  approach  was  to  multiply  the  Boltzmann  equation 

2^4 

successively  by  1,  v,  v  ,  v  ,  v  ,  ...  and  integrate 
over  the  velocities.  The  first  five  moments  yield  the 
usual  transport  equations.  For  these  first  five  moments 
the  collision  integral  is  zero  but,  in  general,  a  given 
moment  equation  involves  a  moment  of  the  next  higher 
order  while  the  collision  integral  contains  the  entire 
distribution  function.  This  set  of  equations  is  equivalent 
to  the  Boltzmann  equation  that  one  began  with.  One  solves 
this  set  by  choosing  a  trial  function  with  n  arbitrary 
parameters  to  represent  the  distribution  function.  It 
turns  out  that  one  may  write  these  parameters  in  terms 
of  the  first  n  moments.  One  can  express  all  other 
moments  in  terms  of  these  lower  moments.  It  is  then 
assumed  that  the  first  n  equations  completely  describe 
the  behaviour  of  the  gas  with  time .  As  his  trial 
function.  Grad  took  the  Maxwellian  distribution  multi¬ 
plied  by  three-dimensional  Hermite  polynomials  (which 


'  ""  •  ;;r.  w  ..  r  ■  «  ips  n  ;  *;;s  "  I of 

o  ’ 

‘  >.  ...'  ©ti:fK 

i:  ai"  :'r  ; xi  -.'i-v  ...  "Ti-.q 

■ 

•.0;-  ■  ■  '  ;K->  .  V.:  ■  :  •  ■  '  •MP 

JI 


36. 


are  orthogonal  with  respect  to  the  Maxwellian  distribution) . 
In  particular,  he  took  thirteen  parameters  in  his  distribu¬ 
tion  function  and  showed  that  the  function  approached  the 
Chapman- Enskog  result  exponentially  in  time. 

Thus,  in  terms  of  this  new  terminology.  Grad 
derived  a  thirteen-moment  -approximation  to  f(r,  v,  t), 
while  Enskog  had  obtained  only  a  five-moment  approximation, 

( It  should  be  pointed  out  at  this  stage  that  the  calcula¬ 
tions  in  this  thesis  will  not  be  as  sophisticated  as  those 
of  Grad.  Only  the  first  five  moments  will  be  considered 
for  reasons  which  will  become  obvious  later  on.)  The 
advantage  of  the  results  of  the  thirteen-moment  approxima¬ 
tion  (besides  the  obvious  one  of  convergence)  is  that  they 
require  a  more  complete  set  of  boundary  conditions  and, 
thus,  one  can,  in  principle,  solve  more  general  problems. 

More  recently,  investigators  have  taken  seven¬ 
teen  moments  in  trying  to  arrive  at  an  expression  for  f. 

As  far  as  this  author  knows,  no  one  has  applied  these 
"more  exact"  distribution  functions  to  the  problem  of 
sound  absorption  in  gases.. 

3 .7  The  Approach  of  Bhatnagar  et  al 


A  completely  different  approach  was  formulated 
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by  P.  L.  Bhatnagar,  E.  P.  Gross  and  M.  Krook^  (195^). 
For  the  collision  integral  they  took  the  following: 


df  1  _  -n(x, t) [f-f ] 

coll.  a(x,t) 


(3.30) 


where  a  (proportional  to  the  inverse  of  the  collision 
frequency  A)  is  velocity  independent.  Here,  n(x,t)  = 

/f dv  and  F  is  the  usual  local  Maxwellian  distribution 
function.  They  wrote  down  the  first  three  moment 
relations  which  define  the  macroscopic  variables  n(x,t), 
u(x,t),  and  T(x,t).  They  then  Fourier  transformed  the 
three  moment  equations,  thus  eliminating  the  spatial 
dependence,  and  setting  up  an  initial  value  problem.  This 
was  suggested  by  L.  LandauJ  (19^6)  who  stated  that  the 
procedure  of  taking  solutions  with  all  the  quantities 
being  of  the  form  ei(kx-oyt)  j_n  Boltzmann  equation 

lacked  rigour.  One  must  also  consider  the  solutions  of 
the  homogeneous  part  of  the  Boltzmann  equation  and  then 
take  the  correct  linear  combination  in  order  to  solve  the 
initial  value  problem.  The  initial  conditions  taken  by 
Bhatnagar  et  al  corresponded  to  compressing  a  gas  slowly 
between  two  walls  (so  that  the  initial  distribution 
function  is  very  nearly  Maxwellian)  and  then  taking  the 
two  walls  away.  Then  two  pulses  will  be  propagated  away 
from  the  origin  in  opposite  directions  and  will  be  damped 
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as  they  move  through  space.  This  is  the  initial  value 
model  that  they  took  for  sound  waves. 

The  moment  equations  were  completely  specified 
when  they  stated  the  Fourier  transform  of  f(v,t)  at 
time  t  =  0.  Following  Landau's  procedure,  they  then 
performed  a  Laplace  transformation,  with  respect  to  the 
time  variable,  on  the  Fourier- transformed  moment  equations, 
thus  eliminating  the  time  dependence.  Because  of  the 
initial  conditions  they  did  not  obtain  homogeneous 
equations  for  the  transformed  variables.  They  then 
solved  the  equations  for  the  transformed  density.  Setting 
its  denominator  equal  to  zero  yielded  the  dispersion 
relation.  Then,  putting  a)  =  -  io^,  they  found 

expressions  for  co-^,  in  terms  of  k  (which  was 

taken  to  be  real)  and  A.  That  is  to  say,  they  found  out 
how  the  two  pulses,  propagating  in  opposite  directions, 
were  damped  out  with  time.  Comparison  with  experiment 
was  arrived  at  by  taking  their  results,  to  the  first 
order  in  A,  exactly  equal  to  the  well-known  Navier-Stokes 
solution,  rather  than  determining  A  explicitly  from 
their  own  model.  Their  dispersion  curve,  obtained  in 
this  way,  was  found  to  lie  between  the  curves  of  the 
Navier-Stokes  and  Burnett  equations.  They  did  not 
calculate  the  amplitude  absorption  for  progressive  sound 
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waves  for  either  low  or  high  frequencies.  Furthermore, 
they  did  not  solve  their  initial  value  problem  for  short 
wavelengths . 

Various  investigators,  such  as  C.  L.  Pekeris, 
Z.  Alterman,  L.  Finkelstein,  and  K.  Frankowski^'  (1962) 
have  considered  the  propagation  of  sound  in  rarefied 
gases  (the  Knudsen  case),  but  not  the  propagation  of 
sound  waves  of  high  frequency  through  a  gas  at  normal 
pressures . 


. 
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Chapter  4.  THE  BOLTZMANN  EQUATION 

Our  Basic  approach  will  be  to  formulate  our 
problem  in  terms  of  the  B'Oltzmann  equation  which  we 
shall  solve  by  replacing  the  collision  integral  by  a 
simple  model.  Then,  we  shall  use  the  distribution  function 
to  deduce  the  general  dispersion  relation  (applicable 
at  any  frequency)  for  sound  propagation  through  mono- 
atomic  gases. 

Since  it  forms  the  basis  of  all  our  calculations, 
let  us  briefly  derive  the  Boltzmann  equation  and  discuss 
its  limitations. 

4 . 1  Physical  Derivation 

Let  f (r,  v,  t)  be  the  phase  density  of  particles 
(of  mass  m)  of  the  gas  with  a  velocity  between  v  and 
v  +  dv  and  coordinates  between  r  and  r  +  dr.  Then 
f (r,  v,  t ) dr  dv  is  the  probable  number  of  particles  in 
the  (6-dimensional)  volume,  or  cell,  dr  dv  of  phase 
space.  f(r,  v,  t)  is  taken  to  be  continuous  and 
analytic.  Thus,  in  each  cell  of  phase  space  being 
considered,  there  must  be  a  large  number  of  phase  points 
(each  point  corresponding  to  a  particle) .  However,  the 
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cells  cannot  be  so  large  as  to  contain  the  whole  gas. 

In  the  absence  of  collisions,  the  particles 

move  In  such  a  way  that  at  a  time  t 1  =  t  +  At  , 

F 

r  r 1  =  r  +  vAt  and  v  -+  v1  =  v  +  AjAt  ,  where  F  Is 
the  external  force  that  Is  acting  on  the  particles. 

At  Is  assumed  to  be  such  that  t  «  At  «  tc  ,  where 
ts  Is  the  time  duration  of  a  collision  and  tc  Is  the 
mean  time  between  successive  collisions.  Thus,  a 
collision  which  begins  during  At  will  generally  be 
completed  during  It  and  no  particle  can  experience  more 
than  one  collision  In  At.  This  implies  that  the 
intermole cular  forces  are  short-ranged.  These  particles 
then  will  occupy  some  other  volume  of  phase  space  dr  *  dv1 
at  t * .  The  probable  number  of  particles  in  this  volume 
will  be  f  (r 1 ,  v ' ,  t ' )  dr 1  dv 1 .  However,  since  we  are 
considering  the  same  set  of  particles, 

f(r>,  y 1 ,  t')  dr1  dv1  =  f(r,  v,  t)  dr  dv  .  (4.1) 

According  to  Liouville’s  theorem  (see,  for  example,  A. 

oQ 

Sommerf eld 1  s-3  (1956),  "Thermodynamics  and  Statistical 

Mechanics")  of  classical  statistical  mechanics,  the 
volume  dr  dv  of  phase  space  will  remain  constant  in 
time,  although  its  shape  may  change.  Using  this  and 
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expanding  the  left-hand  side  of  (4.1)  about  r,  v,  t,  we 
have  (per  unit  time  and  phase  space  element) 


df  ,  df  ,  P  df 
-  '  a?  +  m  '3^ 


0 


(4.2) 


In  deriving  equation  (4.2),  it  has  been  assumed  that  f 
does  not  change  appreciably  In  the  time  interval  At. 


If  collisions  take  place  In  the  gas,  then  not 
all  the  particles  that  left  dr  dv  at  t  will  be  found 
in  dr 1  dv 1  at  t 1 .  The  probable  number  of  particles  that 
are  "lost"  from  the  volume  dr 1  dv’  (or  dr  dv) ,  due  to 
collisions,  is  B  dr 1  dv 1 .  The  probable  number  gained 
(due  to  collisions)  is  -A  dr 1  dvf .  Then,  when  collisions 
are  taken  into  consideration,  (4.2)  must  be  modified,  ie . , 


|£  +  v  .  |£  +  I  .  |£ 

dt  ~  or  m  dv 


A  -  B 


(4.3) 


Now,  what  is  the  form  of  A-B  in  (4.3)?  Consider 
two  particles  with  velocities  v  and  Vj,  respectively, 
before  collision  and  velocities  v1  and  v^ '  after 
collision.  Let  each  collision  between  pairs  of  such 
particles  be  characterized  by  an  impact  parameter  b. 

From  the  conservation  of  energy,  it  may  be  shown  that 
(if  the  intermolecular  force  posesses  spherical  symmetry). 


. 
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|  v  -  Vi  I  =  |  v'  -  1  |  =  g  ,  say.  The  equality  of 

the  impact  parameters  b  before  and  after  collision 
follows  from  the  conservation  of  angular  momentum. 

What  is  the  probability  that  a  particle  with 
velocity  vq  and  impact  parameter  b  will  experience  a 
collision  in  the  time  interval  dt  with  another  particle 
at  position  r  and  with  velocity  v?  If  the  latter 
particle  is  regarded  to  be  at  rest,  then  the  relative 
speed  of  the  former  particle  with  respect  to  the  other 
(both,  before  and  after  collision)  will  be  g.  Assume 
that  the  interaction  force  is  sufficiently  short-ranged 
so  that  one  can  define  a  sphere  of  influence  of  radius, 
say,  s  for  each  of  the  particles  "at  rest" .  Outside 
this  sphere  the  intermolecular  force  is  neglible.  The 
sphere  of  influence  about  each  "rest"  particle  must  be 
shorter  than  the  mean  free  path  of  the  particles  or  else 
one  cannot  really  speak  of  collisions.  For  example,  for 
a  gas  of  rigid  spheres  of  radius  a,  the  sphere  of 
influence  has  a  radius  of  2a.  To  be  consistent,  all 
particles  with  velocity  Vq  and  impact  parameter  between 
b  and  b  +  db  are  regarded  as  point  particles.  These 
particles,  lying  in  the  cylindrical  shell  of  volume 
gbdbdcpdt  (see  Fig.  5)*  will  suffer  a  collision  if  they 
enter  the  sphere  of  influence  of  a  "stationary"  particle. 
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collision 


Fig.  5.  Calculation  of  the  Boltzmann 


Collision  Integral 
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Then,  we  make  the  Stosszahl-Ansatz  that  the  probable  number 
of  collisions  between  the  two  types  of  particles  will  be 


f(r,v,t)dv  •  f (r, vp, t) dvj  •  gbd9drdt 


dN 


and  writing  bdbdcp  =?  I(g,0)dfi  where  I(g,@)  is  called 
the  differential  cross-section  for  the  collision  into  the 
solid  angle  dfi  =  sin@d0dcp  we  have  dN  =  gl(  g,  0)  dfif  f  -^dv^drdt . 
The  exact  relation  between  g,b  and  0  is  determined  by 
the  dynamics  of  the  situation.  Integrating  over  all 
possible  velocities  vj  and  all  solid  angles  dfi,  we 
obtain  the  number  of  particles  that  are  lost  from  the 
velocity  range  v  to  v  +  dv,  per  cubic  centimeter,  per 
second,  due  to  collisions,  viz.. 


Bdv  =  d vff  dfigl(g,  @)ff1dv-L 


(4.4) 


Ie . ,  because  of  collisions,  we  shall  not  find  these 


particles  in  dr 1  dv '  at  t1. 


Similarly,  considering  the  inverse  collision 


process,  we  have 


Adv’=  dv’//d^gl(g,  @)f  'fi'dv^L 


(4.5) 


From  the  Liouville  theorem,  one  has  dv  dv^ 
so  that  equation  (4.3)  becomes 


dv '  dvj ' , 


. 
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H  +  Z  '  H  +  I  •  ll  =  /dv1/angi(g,e)[f'f1>-ff1]  (4.6) 

which  is  the  complete  Boltzmann  integral-differential 
equation.  The  results  expressed  in  (4.4)  and  (4.5)  seem 
to  be  intuitively  obvious,  but  they  really  contain  an 
implied  assumption  that  there  is  no  correlation  between 
the  velocities  and  positions  of  different  atoms.  This  is 
the  reason  that  we  were  able  to  multiply  the  distribution 
functions  together.  This  neglect  of  correlations  is  only 
an  approximation  since  correlations  must  exist.  (See 
reference  50,  page  22) . 

Furthermore,  the  results  of  (4.6)  are  limited 
to  low  densities,  because  we  have  assumed  in  the  derivation 
that  there  are  only  two  particle  interactions.  There  are 
no  three-body  or  higher-order  collisions. 

4 . 2  Properties  of  the  Boltzmann  Equation 

The  most  surprising  feature  of  equation  (4.6) 
is  that  it  is  not  invariant  under  time  reversal.  If  one 
reverses  the  time  direction  so  that  t  -t,  then 
f  (r,  _v,  t)  ->  f*(r,  -v,  -t)  .  Then,  the  f orm  of  the 
equation  for  f*  is  different  than  that  for  f.  The 
Boltzmann  equation  is  irreversible.  This  is  further 


demonstrated  in  the  famous  Boltzmann's  H- theorem.  If 


one  defines 


H  =  / fin  f  dv 


(4.7) 


then  it  is  easily  proved,,  using  the  collision  integral  of 
(4.6),  that  due  to  collisions. 


(— ) 

^  dt '  coll 


^  0 


(4.8) 


=  0,  only  if  f'f-^'  =  ff^,  which,  as  indicated 

earlier,  leads  to  the  local  Maxwellian  distribution  f(°) 
No  matter  what  the  initial  state  of  f ,  the  effect  of 
collisions  is  for  f  to  approach  the  equilibrium 
distribution  f (®) . 


If  one  is  dealing  with  Fermi-Dirac  or  Bose- 
Einstein  statistics  the  term  f'f-^1  -  ff^  would  have 
to  be  replaced  by 

f'f1*(i-f)(i-f1)  -  fqti-f MU-fy) 

and 

f'f  *(l+f)(l+f1)  -  ff1(l+f ') (1+f  ') 
respectively.  That  is  to  say,  the  exclusion  principle 
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must  be  taken  into  consideration. 

4 . 3  Stosszahlansatz  and  Irreversibility 

It  has  been  noted  that  the  Boltzmann  equation  is 
irreversible  under  time  inversion.  If  one  looks  at  the 
Liouville  equation,  which  is  strictly  derivable  from 
classical  mechanics. 


|f  =  {h;d}  (4.9) 

(H  is  the  Hamiltonian  of  a  system  and  D  is  the  density  of 
such  systems  in  T-space.),  one  can  easily  prove  that 
(4.9)  _is  reversible.  The  Liouville  equation  will  yield 
the  streaming  terms 

v  .  D  +  I  .  D 

—  Sr  m  "civ 

of  the  Boltzmann  equation,  so  that  irreversibility  must 
enter  (4.6)  through  the  Stosszahlansatz  about  the  probable 
number  of  collisions. 

Let  us  examine  the  Stosszahlansatz  embodied  in 
equations  (4.4)  and  (4.5)  more  carefully.  The  crux  of 
the  problem  is  contained  in  the  word  probable.  Since  we 
have  a  statistically  large  number  of  particles  in  our 
system,  we  cannot  possibly  consider  the  detailed  motion 


of  each  particle.  Rather,  we  must  talk  about  the 
collective  behaviour  of  groups  of  particles  and  resort 
to  probabilistic  assumptions.  That  is  to  say,  we  cannot, 
for  practical  purposes,  adopt  the  microscopic  viewpoint 
but,  rather,  must  adopt  the  macroscopic  viewpoint.  Thus, 
in  going  from  the  microscopic  viewpoint  of  the  (reversible) 
Liouville  equation  to  the  macroscopic  viewpoint  of  the 
(irreversible)  Boltzmann  equation,  we  must  make  an 
additional  assumption  which  leads  to  the  fore-mentioned 
probabilistic  outlook. 

In  writing  down  (4.4)  and  (4.5) ^  we  took  dr  dv 
rather  than  considering  small,  but  finite,  cells  in  phase 
space,  namely  ArAv  .  If  we  want  to  speak  meaningfully 
of  the  probable  number  of  collisions  in  a  time  interval 
At,  we  cannot  take  the  cells  to  be  arbitrarily  small. 

This  corresponds  to  considering  phase  space  to  be  "coarse¬ 
grained".  Our  level  of  description  has  changed.  We  can 
no  longer  specify  the  state  of  the  particles  within  one 
of  these  finite  cells.  We  cannot  specify  f  at  a  point 
within  a  cell  but  can  only  specify  the  average  value  of 
f  over  a  cell,  viz., 

f  =  A-  -  f  f  dr  dv 

In  considering  the  macroscopic  state  of  a  system,  we 
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cannot  go  below  a  certain  limit  In  specifying  the 
positions  and  velocities  of  the  particles.  Therefore, 
we  can  obtain  irreversible  macroscopic  behaviour  from  the 
reversible  microscopic  viewpoint  by  assuming  that  phase 
space  is  "coarse-grained". 

One  can  explain  the  H-theorem  on  the  basis  of 
the  following  phase  mixing  argument  originally  proposed 
by  Gibbs^-’^  (1903>  1928)  and  emphasized  by  P.  and  T. 
Ehrenfest^  (1912).  Rather  than  considering  a  single 
system,  consider  an  ensemble  of  such  systems.  Let  each 
macroscopic  state  correspond  to  a  region  in  T-space. 

We  are  taking  T-space  to  be  "coarse-grained".  (From  the 
microscopic  viewpoint  each  point  in  T-space  would  corres¬ 
pond  to  a  state.)  Suppose  that  part  of  T-space  corresponds 
to  a  non-equilibrium  initial  state  while  the  remainder 
corresponds  to  equilibrium  states.  By  Liouville!s 
theorem,  this  non-equilibrium  region  will  not  change  its 
volume  but  may  change  its  shape.  Because  of  collision 
processes,  the  T  points  in  the  non-equilibrium  region  will 
tend  to  diffuse  through  T-space.  The  non-equilibrium 
region  will  be  stretched  into  an  irregular,  long, 
fine,  filament.  This  filament  will  eventually  become 
so  fine  that  from  our  macroscopic  viewpoint  all  of 
T-space  will  appear  to  correspond  to  equilibrium  states. 
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Thus,  by  considering  "coarse-grained"  r-space,  one  may 
obtain  a  more  fundamental  understanding  of  the  H- theorem. 

The  pioneering  work  in  attempting  to  go  from 
the  Liouville  equation  to  the  Boltzmann  equation  by  Intro¬ 
ducing  the  additional  concept  of  the  " coarse-graining" 
of  phase  space  was  done  by  J.  Yvon42j,43  (1935,  1937)  and 
M.  Born  and  H.  areen44j4^46^7(  1946,  1947).  J.  Kirkwood48 
(1947)  achieved  irreversibility  by  considering  time  to  be 
coarse-grained.  By  time-averaging,  one  may  smooth  out 
the  sharp  fluctuations  in  the  dynamical  motions  of  the 
system.  N.  Bogoliubov 1 s4^  (1946)  work  follows  the  latter 
path.  He  replaced  the  time  variable  t  by  the  one- 
particle  distribution  function  F^(r,  v,  t)  which 
varies  more  slowly  with  t  than  does  t  itself.  Then, 
for  instance,  the  two-particle  distribution  function 
could  be  written  as  ^(r^,  Vj,r2,  v^Fi  (  t) )  and  so  on. 

The  work  of  the  above  investigators  leads  to  the 
BBGKY  (Bogoliubov,  Born,  Green,  Kirchoff,  Yvon)  hierarchy 
of  equations  which  are  discussed  at  great  length  in  an 
excellent  account  by  T.  Y.  Wu82>  (1961). 

4 . 4  Limitations  of  the  Boltzmann  Equation 


(i).  To  make  the  concept  of  an  instantaneous 
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collision  meaningful,  one  has  assumed  in  deriving  the 
collision  integral  that  the  intermolecular  forces  are 
short-ranged.  For  long-ranged  forces,  such  as  the 
Coulomb  interaction,  one  has  to  omit  the  collision  term 
and  consider  the  particles  to  move  independently  through 
an  average  potential  field  Veff  which  is  the  sum  of 
the  applied  field  and  the  average  field  produced  by  all 
the  particles  in  the  system.  This  leads  to  the  so-called 
Vlasov- Landau-^ 3  ^  (1938,  1946)  equation 


CS?  +  * 


df 

dr 


^eff .  _d  r 
dr  dv 


0 


(4.10) 


More  recently,  R.  L.  Guernsey  and  G.  Uhlenbeck^2  (i960), 

using  the  BBGKY  equations,  obtained  a  more  general 

relation  with  the  Vlasov-Landau  equation  being  the  zeroth 

order  approximation.  Work  in  this  direction  has  also 

been  done  by  Wu  and  Rosenburg^^  (1962)  and  by  I.  Prigogine 

and  R.  Balescu'^-’  ^  ^  (1959*  i960).  A  synopsis 

of  all  the  material  covered  in  the  fore-mentioned  papers 

SO 

may  be  found  in  Wu  1  s  paper. 


(ii).  Another  reason  for  limiting  the  Boltzmann 
equation  to  weakly  interacting  systems  is  that  it  yields 
the  transport  equation  corresponding  to  the  conservation 
of  kinetic  energy.  However,  in  general,  it  is  really  not 
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the  kinetic  energy  that  is  conserved,  but,  rather,  the 
total  energy.  The  Boltzmann  equation  is  only  valid 
when  the  average  potential  energy  is  much  smaller  than 
the  kinetic  energy  ( ie . ,  in  the  weak- interact  ion  limit). 

(iii).  As  pointed  out  earlier,  because  of  the 
assumptions  employed  in  deriving  the  collision  integrals, 
the  Boltzmann  equation  can  only  be  applied  to  dilute  gases . 
For  dense  gases  and  for  liquids  one  would  have  ternary 
and  higher-order  collisions.  The  generalization  of  the 
Boltzmann  equation  to  higher  densities  has  been  carried 
out  in  a  systematic  way  by  Bogoliubov. 

For  a  more  extended  discussion  of  the  limitations 
of  the  Boltzmann  equation,  one  is  referred  to  the  text 
"Quantum  Statistical  Mechanics"  by  L.  P.  Kadanoff  and 
G.  Baym60  (1962). 

Because  it  is  of  no  direct  concern  to  us  here, 
we  have  omitted  any  discussion  of  the  relativistic  or 
quantum  mechanical  aspects  of  the  Boltzmann  equation. 
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Chapter  5.  OUR  MODEL  FOR  THE  BOLTZMANN 


COLLISION  TERM 


5 . 1  The  Relaxation  Model 


Rather  than  considering  the  intermolecular 


forces  explicitly,  we,  instead,  write  for  the  collision 


in  the  Boltzmann  equation 


(f-n 

T(v) 


(5.1) 


dt; coll. 


where  f  is  the  local  equilibrium  (displaced  Maxwellian) 
distribution  when  the  sound  wave  is  propagating  through 
the  gas 

?  -  Mf-  tlf  m  1 3/g  -jm/[2KT(z,t)]}|v-u(z,t)  \2  , 

and  u  is  the  macroscopic  or  flow  velocity  imposed  on 
the  particles  by  the  sound  wave.  t(v )  is,  in  general,  a 
velocity-dependent  relaxation  time. 

5 • 2  Properties  of  the  Model 


The  chief  characteristics  of  our  model  are: 


(i)  The  model  recognizes  the  collective 


behaviour  of  the  gas  particles  when  the  sound  wave  is 


Djf: 


■  1 0  G 
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propagating  through  the  medium.  Particles  with 
velocities  greater  than  the  flow  velocity  will  make 
collisions  which,  on  the  average,  slow  them  down  to  the 
flow  velocity  at  the  point  in  question.  Conversely, 
particles  with  velocities  less  than  the  flow  velocity 
will  tend  to  have  their  velocities  increased  to  the  flow 
velocity. 


(ii)  By  taking  the  collision  term  to  be  of  the 


form  (5.1) j  we  are  "smoothing  out"  the  actual  decay  of 
the  distribution  function  to  equilibrium. 


(iii)  Let  us  write  the  distribution  function 


as 


f(r,v,t)  =  f(0)(v)(l  +  <p(r,v,t)]  (5.3) 


where  |  cp  |  «  l,  and 


_  22 

It  was  first  shown  by  D.  Enskog  (1917)  and  reproduced 
by  L.  Waldmann^l  in  the  Handbuch  der  Physik,  12,  (1958) 
that  the  collision  term 


If  dvl(m  •  gl(g,0)[f >f1’-ff1] 


=  /  K(v,v-l)  [cp(v1)-cp(v)  ]dvx 


(5.4) 


;qo  lq 
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where  K(v,v^)  =  K(v^,v)  is  referred  to  as  the  (symmetric) 

"collision  kernel".  This  must  be  equal  to  our  collision 

term.  As  has  been  shown  by  H.  Grad  (1958),  our  model, 

containing  a  single  relaxation  time,  may  be  crudely 

obtained  from  the  general  Boltzmann  collision  term.  If, 

in  our  model,  we  had  considered  a  spectrum  of  relaxation 

times  t ^  ,  we  would  be  able  to  solve  the  Boltzmann  equation 
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uniquely  (see  reference  1  and  also  Dettmann's  "Mathe¬ 
matical  Methods  in  Physics  and  Engineering",  p.  23^-258, 
(1962)),  but  we  have  considered  only  a  single  relaxation 
time  t.  Our  model  leads  to  the  result  that  the  solution 
of  the  Boltzmann  equation  is  not  unique,  ie .  one  should 
really  take  a  linear  combination  of  all  possible  solutions. 

As  pointed  out  by  Gross^  ( 1957 ) ,  our  model  implies  that 
the  "collision  kernel"  is  cf  a  degenerate  nature.  This 
simplifies  the  analysis  but  has  the  consequence  that  one 
cannot  expect  quantitatively  exact  results.  However,  one 
should  expect  qualitatively  correct  behaviour  on  the 
basis  of  our  model. 

(iv)  For  a  one- component  system,  such  as  a 
monoatomic  gas,  we  require  that  the  usual  conservation 
laws  be  satisfied. 
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Multiplying  the  Boltzmann  equation  by 


and  defining 

<  f  >  =  /^fdv//fdv  =  ^J^/f dv 

one  has 

[N]  +  ^|[N<vz>]  =  0  =  -/1hf-f)  dv  (5.5) 

gf  [N<vz>]  +  g|  [N<vz2>]  =  0  =  -/  — dv  (5.6) 

[N<v2>]  +  g|[N<v2vz>]  =  0  =  -J  v---^-Z)  dv  (5.7) 

(v)  The  Boltzmann  equation,  with  the  collision 
term  represented  by  (5.1).,  is  still  irreversible.  Ie . , 
the  H- theorem  is  still  satisfied.  If  time  is  reversed, 
one  obtains  the  so-called  anti-Boltzmann  equation 
(corresponding  to  ^‘5t)coli  =  — '  ~  ~  )  whose  signifance 

has  been  discussed  in  length  by  E.  R.  Adams  D  (i960), 

E.  G.  Cohen  and  T.  H.  Berlin^  (i960),  and  D.  K.  C. 
MacDonald^  (1961). 
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5 . 3  Further  Specialization  of  the  Model 

We  shall  further  specialize  (5.1).,  for  the 
purpose  of  simplifying  our  problem,  by  considering  two 
special  cases: 

(i)  t  to  be  independent  of  v  .  This  will 
be  discussed  in  Chapters  7  to  11. 

(ii)  t  to  be  of  the  form  t  =  $, /v  where  the 
mean  free  path  £  does  not  depend  on  v. 

This  case  will  be  considered  in  Chapter  12. 

The  former  assumption,  (i),  is  made  to  simplify 
the  integration,  while  the  latter  assumption,  (ii),  is 
made  in  an  attempt  to  refine  the  results .  As  has  been 
pointed  out  by  H.  Grad^,  no  one  has  systematically 
applied  the  relaxation  model  (5-1)  to  the  propagation  of 
sound  (of  arbitrary  frequency)  through  monoatomic  gases. 
Bhatnagar  et  al,  whose  work  has  been  previously  discussed, 
use  a  similar  model  but  their  procedure  and  aim  is 
different  than  ours. 
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Chapter  6.  SOLUTION  OF  THE  BOLTZMANN  EQUATION 


Consider  a  classical,  dilute,  monoatomic  gas  in 
which  a  plane  sound  wave  of  the  form  ei(kz-oyt)  iS 
propagating.  The  complex  wave  number  k  may  be  written 
as 

k  =  k-^  +  ik2  =  ^  +  ia  .  (6.1) 


Since  the  wave  is  periodic  in  the  z-direction,  we  expect 
the  same  periodicity  to  exist  in  all  other  observables. 
Thus,  the  temperature  of  the  gas  will  be 

T  =  TQ  +  T1  (6.2) 


The  macroscopic  velocity  u(z,t)  imparted  to  the  particles 
by  the  sound  wave  is 


u  = 


i(kz-o)t) 

uQ  « 


(6.3) 


The  number  density  will  be 


N 


N0  +  N1 


ei(kz-oyt) 


(6.4) 


The  distribution  function  is 


f 


r  +  r  i(kz-o)t) 
0  1 


(6.5) 
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f 


0 


( — - — ) 

0 ' 2ttKT  > 


3/2 


e-mv2/2KTc 


(6.6) 


Here  f^,  T^  and  are  the  equilibrium  values  in  the 
absence  of  the  sound  wave.  Furthermore ,  by  the  very 
definition  of  a  sound  wave 


VNq  «  1  ;  T^/Tq  «  1  5  fqAo  «  1  (6- 7) 

Assuming  that  there  are  no  external  forces 
acting  on  the  gas  particles,  we  may  write  the  Boltzmann 
equation  for  our  relaxation  model  as 


3f  ,  df  (f-F) 

+  V7  =  -  -i - 

Ot  ^  OZ  T 


(6.8) 


Since  there  is  a  mass  velocity  u(z,t),  the  local 
equilibrium  distribution  f  must  be  of  the  form  of  a 
displaced  Maxwellian  distribution 


f  =  N(z  t)r  m  ,  I3/2  g  tm/2KT( z.  t)]  [(vz-u(z,t))2+v  2+v  2] 
1  -  I'Hz<tH2irK'Il(z,t)  i  e.  y 

(5.2) 

This  local  equilibrium  distribution  will  clearly  have  the 
periodicity  of  the  imposed  sound  wave.  Thus, 

f  =  fQ  +  e1(kz-<ut)  (6.9) 


Vf0  ^  1 


(6.10) 
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Then,  using  (6.5)  and  (6.9)*  equation  (6.8)  becomes 


-icof1  +  ikvzf1 


-(fl-fj) 

T 


or  solving  for  f^ 


(6.11) 


f. 


( l-ioyr+ikTvz) 


so  that 


f  = 


f  0i  ( kz.-cot ) 
f0  +  ( l-io)T+ikTvz) 


(6.12) 


(6.13) 


Now,  making  use  of  conditions  (6.7)  and  (6.10),  we  may 
write 


f 


f0  + 


f0*-N0  +  2Povzuo 


(6.14) 


Where  Po  =  2KT 
neglected  60u02 

identifying  f^ 


0 

in 

in 


In  deriving  (6.l4),  we  have 
comparison  with  unity.  Then, 
(6.l4),  we  finally  have 


on 


fotiS  +  +  (^i  +  PqV2)^] 


-3 


f  =  f  + 

0 


Tl. 

O' 


(l  -  icoT  +  ikxv  ) 


(6.15) 


Thus,  we  have  found,  in  a  very  simple  manner, 
the  distribution  function  describing  the  behaviour  of 
the  gas  as  a  sound  wave  of  frequency  co  passes  through  it. 
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One  can  see  that,  within  the  limitations  of  the  model, 
f  may  be  applied  for  any  arbitrary  frequency.  We  shall 
systematically  use  (6.15)  to  determine  the  behaviour  of 
the  sound  wave  for  any  cot  (assuming  the  gas  pressure 
to  be  held  fixed) . 
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Chapter  7 .  LOW  FREQUENCY  SOUND  WAVES  :  cdt  «  1 


7.1  Dispersion  Relation 


At  low  frequencies  |  layr  -  ikTVz  j  «  1  so  that 
we  may  expand  the  denominator  of  (6.15)  in,  essentially, 
powers  of  t 


r  N-i 

f  =  fo  +  fo^  +  2Povzuo  + 


ei(kz-o)t) 


[l  +  ieDT  -  ikTV„  -  CD2T2  -  +  2kO)T^Vr 


-  1cd^t3  +  3ikflD2T3v  -  3ia)k2T3y  2  +  ik3x3y  3]  (7-1) 


In  (7-1),  we  have  kept  terms  out  to  t3. 

By  multiplying  equation  (6,8)  by  1,  vz  and  v2, 
we  obtain  the  number,  momentum  and  energy  conservation 
equations  respectively 

|t  [N]  +  [N<vz>]  =  0  (5.5) 

|t  [N<vz>]  +  [N<vz2>]  =  0  (5.6) 

|t  [N<v2>]  +  [N<v2vz>]  =  o  (5.7) 
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where  (as  in  Chapter  5)  we  have  defined 


<ip>  =  ftfdv/ffdv  =  jj.fi/  fdv 


Then,  using  (7.1)  in  equation  (5.5) *  we  have 


N 


k  ■  ■-  -2  _  3  11c3t3 


u  [—  +  ikT  -  kcox2  -  ik<D  t-5  -  - - 

olcd  2  6 


o 


-i  v3- 


-nJT3 


1  ik3T  1  k^T^  ,  1  ik' t-3  ,  1  ik3x3 

2  P0a,2  '  *  ^3  B  ^37  3  p  2^2  J 


.  £jL  r  1  ik2x  k2x2  _3  ik2cpT3  ,  _1  ik^T3 

To  2  ^o03  Po  2  Po  2  602oj 


1  k^T2  ,  1  ik6x3  n 

tf  r^"2  +  «  — 3” 3  J 


(7-2) 


for  the  number  conservation  equation.  Employing  this 

l/T0 


result,  we  can  express  f  in  terms  of  uq  and  Tn/T 


only 


f  = 


f  +  f  [u  (—  +  2ikx  -  3hcDX2  4ikcu2x3 
o  o  o  1 0 

J1-2  2  2  lk3^  1  k^2  +  1  ik7T3  +  26  v 

'  2  Poo;2  ft  PoV  8  o?  +  2Povz 


T. 


o 


+  —  Iv  T  +  2P^i(DTV^  -  ■  -  TV^ +3k2T2Vz  +  6ik2CDT3V, 


2  poco 


O' 


CD 


■  *  SS Vz + *  sS Vz  *  2P°ikTVz2  ■  2P°“2t2v: 


65. 


#  V  .  ,lk3,v  +  i  -  ayAV 


4  3 


-  ~ ---  v„  +  4Bnkayr2v  2  -  2B  iaD3T3v  +  6iP  koo2T3v  2 

^o0)  z  K°  Z  O  z  O  z 


6iP0a)k2T3vz3  +  ■ -  vz3  +  2P0ik3T3vz^} 


,  T1  f  1  ik2T  ,  3  k2x2  ,  3ik2coT3  ,  1  ik^T3 

+  t;  1-  2  +  2  —  +  — —  + 1  TV 


4  2 


.6t3 


1  k  T  +i  ±Klj  .  3  +  p  v2 


*  6o2032  8  ^3  "  2  "a 
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1  k3T2  „  2ik3T3  „  ,  1  ik5T3  „  ,  3 

2P7  z‘  2  11  W  Z  2  2 


-  P0ikTV2vz  +  cd2t2  -  Pqcd2t2v2  +  vz2 


+  —  k2T2Vz2  -  PQk2T2V2Vz2  -  3k(DT2Vz  +  2pQkCDT2V2V, 


+  -§  ia)3T3  -  iP_o)3T3v2  -  ikao2T3v^  +  3Pnikao2T3v2 


Vr 


+  ||  iajk2T3vz2  -  3iP0cjDk2T2v2vz2  -  ~  ik3x3vz3 


+  P0ik3T3v2Vz3  J  ]  e1^2"^ 


(7-3) 
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The  reader  will  probably  throw  up  his  hands  in  horror  at 
this  rather  lengthy  expression.  Rather  than  eliminating 
NiAo  stage,  we  could  have  used  (7-1)  to  write 

down  all  three  (homogeneous)  equations  for  number, 
momentum  and  energy  conservation.  For  a  non- trivial 
solution,  we  would  have  obatined  a  3*3  determinant 
relating  k  to  cot  .  However,  as  it  turns  out,  this 
would  have  been  just  as  laborious  as  the  procedure 
adopted  here.  Using  (7.3).,  we  can  now  obtain  two 
homogeneous  equations  (corresponding  to  momentum  and 
energy  conservation  respectively)  in  the  two  unknowns 
uQ  and  T-j/Tq.  Viz., 

Momentum  conservation: 

uQ[-4i(jo3  +  4co^x  +  2ik2x2co  +  k^x^x  -■  ik^x^cDX2  +  4icD3x2 

+  6ik2x2co3x2  -  i  _  i6k2x203^x3  -  k^x^x3 

2  03 

.  4o>6t3  -  j  ^4^  ]  +  ^  [  2ikx2u>2  +  lAtoT 

4  00^  Tq 

-  4kx2<D3x  -  2ik3x^cD2x2  -  ik3x3x2  -  6ikx2oc>^x2 

+  l4k3x4co3x3  -  i  k5x6(Dx3-  i  +  8kx2o35x3]  =  0 

2  4  a) 


(7.4) 
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where  x  =  l/v/P^  •  Energy  conservation: 

u  [2ik  -  4kcDT  +  k  x  —  -  6ikcD2T2  -  ^  XrvT  ■  +  8kcD3i;3 

o  CD  2 

+  l4k3K2o>T3  -  1  ^4^  -  2ik2X2T2  -  I  ] 

^0)0  2  CD 

+  [  I  k2x2T  +  i(uk2x2T2  -  ^  ■%  -I- 

-  21  k2x2o>2 t3  -  35  k4x4T3  .  l  _  31cd  +  3a)2T 

+  31o)3t2  -  3cd4t3]  =  o  (7.5) 

For  these  last  two  equations,  (7.4)  and  (7.5)*  to  have 
a  non-trivlal  solution,  the  determinant  of  the  coefficients 
of  the  unknowns  uQ  and  T]_/T0  must  be  set  equal  to 
zero.  Doing  this  and  collecting  terms,  we  have  to  the 
third  order  in  t 

23k^x^cD2T2  +  48ik4xij-cD3T3  +  2ik6x6  cdt3  -  lOk^x^cD^ 

-  8ik^x^cD3x  -  28k^x^cD^T^  -  120ik^x^cD^T3  +  12cd^ 

+  241cd3t  -  36cd^t2  -  481cd^t3  =  0  (7-6) 

which  is  the  dispersion-absorption  equation  for  cdt  «  1. 
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7 . 2  Absorption  Coefficient  and  Velocity  Dispersion 

We  shall  now  determine  the  amplitude  absorption 
coefficient  o t  and  the  phase  velocity  V  corresponding 
to  different  orders  in  t. 

t°  :  To  this  order,  we  have  from  (7.6) 

-  10k2x2co2  +  12co4  =  0  (7.7) 

Now,  k  =  k^  +  ik2  =  ^  +  ia  -  (6.1) 

For  cot  «  1,  we  assume  that 

k2  «  (7.8) 

Then,  (7*7)  becomes 

-10x2co2  [k-L2  +  2ik1k2]  +  12  oo4  =  0 

The  real  part,  on  being  equated  to  zero,  yields 

V2  =  5x2  .  5^ 

°  3m 

or  V  =  V  =  (22°) 1/2  (7.9) 

3m 

which  is  the  usual  adiabatic  velocity.  Equating  the 
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imaginary  part  to  zero  gives 

a  =  0 

to  the  zeroth  order  in  t.  This,,  of  course,  justifies 
our  assumption  (7.8). 

t1  :  -  10  k2x2ac2  -  8  ik2x2a)3T  +  12  oft 

+  24  icD5T  =  0  (7.10) 

Using  (7.8)  and  taking  k2T  to  be  of  second  order  in 
smallness,  (7. 10)  becomes 

-  10  x^cd^  (k^2  +  2  ik-^kg)  -  8ik-j^x2CD3T 

+  12  oft  +  24  loPr  =  0 

The  real  part  again  yields 

V  2  V  =  (5KTq)1/2  (7.9) 

o  3m 

while  the  imaginary  part  gives 

-  20  x^~o7~  •  k-j_k2  -  8  +  24  cjo^t  =  0 

(7-11) 

Then,  with  (7-9)  and  the  definition  k-^  =  cd/V  ,  we 

have  from  equation  (7. 11) 


/  3f/fc  '  ••  "  'V  .■ 

.  ,  j:.oi.  C'.i.'.  rt'JC 


t  <m[-  o  :■  '  o:  ;  ;  4.,-  ‘  Xn  (3.  Y)  df--*  ‘ 

■  '  }  ,a  aenJ  r ; >. 
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k2  =  «  =  (7-12) 

J  o 

This  result,  (7.12),  justifies  our  assumption  since  t2 
terms  are  neglected  in  this  calculation. 

_t£_:  23  k^x\)2T2  -  10  k^x^co^  -  8  ik^x^a^x 

-  28  k2x2o/^T2  +  12  o/^  +  24  io y’l  -  38  cd^t2  =  0 

(7.13) 

We  now  want  to  keep  terms  out  to  t2.  The  pattern 
suggested  by  the  results  obtained  so  far  is  quite  clear. 

p 

V  will  contain  a  term  in  t  while  a  will  not.  Whan 
we  go  to  the  third  order  in  t,  there  will  be  no  contribu¬ 
tion  to  V  but  there  will  be  some  to  a. 

To  the  second  order  in  t,  (7.13)  becomes 
23  k-^x  o)2t^  -  10  [k-^ 2  +  2ikjk2  -  k22]  x2o2 
-  8i[k]_2  +  2ik]_k2]  x2o)3t  -  28  k]_2x2a/t'T2 
+  12  a/1  +  24  IuPt  -  36  co6t2  =  0  (7.14) 


which  gives 
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l6  k^k2X^co^T  -  28  k-^2x2CD4T2  +  23  k^x^o}2T2 
-  10  ki2x20)2  +  ]_o  k22x2<x>2  +  12  a/*-  -  38  cjd^t2  = 


-  8  k-^x^apT  -  20  k]_k2X^co^  +  24  ar*T  =  0 

As  we  have  just  shown,  k^  contains  no  terms  in 
So,  letting  V2  =  x2  in  (7.16),  we  have  to  the 
order  in  t 


a 


3  cd2t 
5  VQ 


ie .  a  does  not  contain  any  terms  in  t2 


Using  this  result  and  writing 


V  =  V  [1  +  Tdt2]1/2 
o 


we  have  from  (7.15)* 


b  =  7^,  a?  =  1.72  cu2 

(12)(25) 


so  that  for  cdt  «  1 


V  =  VQ[1  +  0.86  cd2t2] 


3  0)2T 

5 


0  (7.15) 

(7.16) 

T. 

second 

(7.12) 

(7.17) 

(7.18) 


and 


a 


(7.12) 
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to  the  second  order  in  t.  We  need  not  have  assumed  V 
to  be  of  the  form  postulated  in  (7-17)*  but  to  arrive 
at  the  same  result  without  doing  so  involves  more  labour. 

_t£  :  Equating  the  real  and  imaginary  parts  of 

equation  (7*6)  equal  to  zero,  we  have  (to  the  third 
order  in  t)  for  the  real  part: 

23  k^x^o>2T2  _  io  k-^2x2a)2  +  10  kr^-x^up- 

+  16  kjkpX^a^T;  -  28  +  12 

-  36  cd^t2  =  0 

which  is  just  (7.15)  again.  For  the  imaginary  part: 

92  k2^k2CJo^r^x^  +  48  o^T^k-^x^  +  2  cDT^kj^x^ 

-  20  cD2k]_k2X2  -  8  ci)3Tk-]_2x2  +  8  a)3xk22x2 

-  58  a/^T^kjkgX^  -  120  ar>T;3k12x2  +  24  gd^t 

-  48  o7t3  =  0  (7.19) 

From  (7.15) j  we  once  more  obtain 

=  VQ[l  +  0.86  CD2  T2  ] 


V 


(7.18) 


,v  •:  .  b  '  t 


'  ;  ■  ..  ;  %  '  •  ’  ■  •  -  '  "  ""  '  ■' 


:  p  i  J >  •  ..  pS 


■  f1  V  <:>£  - 

,  q  h  ■  : 

;  '•  c- .  a  t>'i6iTi  ?o  :xC  c'-!-  - 


73. 


Using  this  result  and  writing 

k2  =  a  =  3  o£r  tl  +  Bt2j  (7.20) 

vo 

we  find,  on  solving  (7.19)  to  the  third  order  in  t. 


B 


-299^ 

900 


-3.33  o? 


so  that 

a  =  |  ^  [1  -  3.33  o)2t2]  (7.21) 

5  vo 

Summing  up  the  results  of  our  low  frequency  calculations, 
we  have  found  that 

a  =  [1  -  3-33  cd2t2]  (7.21) 

5  vo 

V  =  v0[l  +  0.86  cd2t2]  (7.18) 

where  Vc  =  (f  ^)1/2  (7.9) 


We  see  that  as  the  frequency  increases  the 
phase  velocity  increases  while  the  absorption  coefficient 
initially  increases  as  the  square  of  the  frequency  but 
then  more  slowly.  We  want  to  compare  these  low  frequency 
results  (cot  «  1)  with  those  calculated  by  Wang-Chang  and 
Uhlenbeck^0  (1948)  for  Maxwellian  molecules. 
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In  Fig.  6,  we  have  compared  both  our  results  and 

theirs  with  experiment  for  «  1.  It  is  seen  that  in 

Vo 

this  region  our  dispersion  curve  —  displays  better 
agreement  with  experiment  than  does  the  Wang-Chang- 

aVQ 

Uhlenbeck  result.  However,  for  the  absorption  curve  — —  , 

the  Wang-Chang-Uhlenbeck  relation  has  the  "correct" 

(classical)  slope  while  ours  is  about  lb%  too  small.  As 

has  been  previously  mentioned/'  the  fact  that  our  absorption 

curve  does  not  have  the  classical  slope  is  not  unexpected. 

If  we  shift  our  absorption  curve  upwards  on  the  log-log 

scale  of  Fig.  6  so  that  its  slope  coincides  with  that  of 

the  classical  absorption  curve,  we  see  that  the  curvature 

avn 

of  our  ——  curve  is  less  than  that  of  Wang- Chang  and 
Uhlenbeck1 s.  Ie . ,  it  would  be  a  better  fit.  In  Chapter 
12,  we  shall  make  a  calculation  which  will  give  us  some 

av0 

improvement  in  the  slope  of  — —  . 

We  shall  now  solve  the  problem  of  sound 
propagation  for  any  frequency  by  obtaining  the  general 
dispersion  relation  valid  for  any  got  .  For  got  «  1 
we  should  obtain  the  results  already  derived.  Then,  we 
shall  attempt  to  solve  the  general  dispersion  relation 
for  got  »  1.  (All  these  calculations  will  be  done  for  t 
independent  of  v  ,  of  course.) 
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See  p.  56 
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Chapter  8.  GENERAL  DISPERSION  RELATION  FOR  ANY  cdt 

8 . 1  Derivation 

In  this  sect  ion ,  we  shall  determine  the  general 
dispersion  relation  valid  for  any  value  of  o)T.  The 
distribution  function  is  the  same  as  before*  namely 


f 


f  + 
o 


O  O 


(  1  -  iCDT  +  ikTVz) 


(6.15) 


Then*  using  this  distribution  f  in  the  conser¬ 
vation  equations  ( 5  •  5)  -  ( 5 .7)  .>  we  have 


0 


(8.1) 


the  continuity  equation. 


0 


(8.2) 
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the  conservation  of  momentum  equation 


In  NX 


dl- 


tu)  ^  • k  sp;1  %  +  t2P°“  3fT  -  *0*  lg]  uo 


a2i0  o  ai-i  ^2j  T 

PoC0  ^7  +  2 k  ^  +  p°k  -§p7]  t" 


=  0 


o  o 


(8.3) 


the  conservation  of  energy  equation,  where 


+  00 

Im(k,o))  =  /// 


vz*e-PoV‘ 


1  -  io)T  +  ikxv. 


dv  (8.4) 


and 


+00 

Iff 

-00 


v2vzme 


-Pnv2 


1  -  iooT  +  ikTVr 


dv 


dl 


m 


(8.5) 


to  I  . 
o 


In  Appendix  I,  we  have  related  1^,  I^,  and 


-i  (^)3/2  +  —  I 

kr  °o  y  c 


(1.13) 


tJL) 3/2  .  _i 

V  a  )  .  o 


Ykx  P0  y2  o 


(1.15) 


1  ,  w  >  3/2 


(+) 
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where 


Then, 


7  =  kT/(  1-1<jdt) 


(1.3) 


Si- 


3  1  t  ir  s  3/2  i  dl0 

2  P0kT  +  Y  WZ 


(8.6) 


3l2 


z2  l  iIL.\ 3/2  1_ 

2  P^Ykr  (3Q'  '  ^ 


_  9xo 

cy2 


(8.7) 


S2I- 


-15  l  r_L^3//2  +  1 

8Q2kT  Po  7  ^P02 


(8.8) 


Furthermore,  from  Appendix  II,  we  have 


dl 

W, 


/  1  i  1  \  j  1  /77"n3/2 

~  +  Jo  "  7k^  (F") 


Po  V 
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(tm 
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(II. 9) 


Then,  using  (1.13),  (1.15),  (1.17),  (8.6), 
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(8.7),  (8.8),  (II. 8)  and  (II. 9),  equations  (8.1)  -  (8.3) 
simplify  to 


[l 


io)T  -  Qq] 


Ni 

N 
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[21P0 
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<g  % 
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r  1  kr 
L  2  7 


Qr 


PokT  ^ 

+  "TT"  Qo 


0  (8.9) 


[kT  '  4  Q°] 


Nl 

N 


+  [1  - 


2(3 


o  ,  2^o 


7kx 
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On] 


U, 


,  r  i  0  4-  1^° 

+  L —  Q0  +  — s- 
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7  ki 


T- 
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Q  =  (^2) 3//2  I 

^O  K  1 T  '  0 


where 
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If  we  solve  the  continuity  equation  (8.9)  for 
N^/Nq  in  terms  of  uQ  and  T-^/Tq  and  substitute  this 
result  into  equations  (8. 10)  and  (8.11),  we  obtain  for 
the  momentum  and  energy  conservation  equations  respectively. 

f4Pai0>  [;L  _  (l.la)T)  qj  +  2k2  [1  -  Qc]Ju0 

+  kaj  [Qo  +  2!°  Q  -  ^2]  li  =  0  (8.13) 

yd  0  ykx  T0 


and 


-P 
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k 


IX  + 
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%  -  ^2  ] 

J d  Ykr 


u. 


+  c 


11 


~ q:  Qo 


Or 


Pq2  .  3  ,  £0  0  _  i  £°_  _  i  £^_  0  -  p°2  0 
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^r 


Qo  - 


p, 


73kx 


^1 

To 


=  0 


(8.14) 


For  these  two  homogeneous  equations  to  have  a  non-trivial 
solution,  we  must  have 
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Y  kr 


+  ^Q, 


y 


7+ 


fV 


3 

y  kT 
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(8.15) 


Expanding  this  and  collecting  terms,  we  obtain 
on  dividing  through  by  k2  and  writing  (3q  =  l/X2_, 
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17  ,  2  leu 

2  ^3^x2 


k^x^T2"^ 


i  °  P 

k'-x  -  r 


+  3  + 


k2x2T2 


2o)c 


7k3x^' 


=  0 


(8.16) 


Defining 


kXT 

a  -  Yx  -  T^Tmr) 


(8.17) 


(8.l6)  becomes,  on  multiplying  through  by  a\ 
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(8.18) 
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where 


Q, 


'o 


o 


la(  1-Icdt) 


1 


(8.19) 


from  Appendix  II. 

8 . 2  Discussion 

Thus,  we  have  obtained  the  general  dispersion 
relation,  relating  k  to  cdt.  In  principle  this  complex 
equation  (8.l8)  may  be  solved  numerically  on  a  computer 
for  any  value  of  cot.  We  can  write  the  integral  Qq 
as  Qq  =  A  +  iB  where  both  A  and  B  are  real.  Then, 
we  can  separate  real  and  imaginary  parts  of  (8.l8)  to 
obtain  two  extremely  complicated  equations  for  k-^  and 
k^  in  terms  of  cdt.  These  equations,  which  are  not  readily 
amenable  to  simplification,  are  indicated  in  Appendix 
VIII. 


We  can  relate  Qq  ,  as  given  in  (8.19),  to  the 

69,70,71 


well-known  plasma  dispersion  function 
defined  by 


Z  ( i/a ) 
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1  +00  _y^  1 

— Yj2  /  ~ — i  for  Im  (•£)  >  0 

-JT  /  ^  -00  y  -  — 

(8.20) 


and  as  the  analytic  continuation  of  this  for 


Then, 


Q 


o 


— 1 - z 

ia(l-icoT) 


(8.21) 


In  Appendix  III,  we  discuss  the  various  properties  of 
Z (“)  which  shall  be  used  shortly. 

ci 

In  reference  69,  B.  D.  Fried  and  S.  D.  Conte 

(1961)  have  tabulated  Z(— )  only  for  a  narrow  range  of 

a 

the  complex  variable  but  have  indicated  general  methods 
of  extending  the  tabulation.  Thus,  if  one  has  the 
patience  and  time,  he  can  find  the  values  of  k-j_,  (for 
a  given  cdt)  which  satisfy  the  equations  outlined  in 
Appendix  VIII.  In  this  thesis  we  have  not  delved  into 
any  numerical  work,  but,  rather  have  only  considered  the 
two  extreme  cases,  cdt  «  1  (for  which  we  already  have 
results)  and  cdt  »  1. 

In  the  next  chapter  we  shall  find  out  what 
our  general  dispersion  relation  yields  for  cdt  «  1. 

Of  course,  we  should  just  expect  to  get  the  results  of 
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Chapter  7  for  k^  and  k^.  Thus,  this  calculation  will 
serve  as  a  partial  check  on  the  validity  of  the  general 
dispersion  relation. 
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Chapter  9.  SOLUTION  OF  THE  GENERAL  DISPERSION 
RELATION  FOR  <jdt  «  1 


In  Appendix  III,  the  asymptotic  expansion  of 
)  for  |  a  |  «  1  is  derived.  Using  this  result,  we 
have  for  Q. 
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where 


a  =  0  Im(-)  >  0 
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2  Im(|)  <  0 


(9.2) 


For  cut  «  1,  we  note  from  experiment  that 
a  I  «  1  and 
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The  results  that  we  obtain  for  k-^  and  k0 
must  satisfy  these  conditions.  Assuming  that  this  Is 
true,  we  have  here 
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16 


.8 


945  . 10  ,  10,395  .12 
“32  a  +  ”64“  a 


...  ] 


(9.3) 


Q 


0  (1-1cdt)^ 


1 -  d-a^l^-fa6^^ 


1095  .10  .  23,535  -12  1 

~W  a  +  m  a  •  •  •  • J 


(9.4) 


Q, 


3 


(1-1oot)8  2 


[l  -  ^  a2  +  3  -  8  aD  + 


6  4-  Ml  08 


T6 


3807  -10  ,  1251  -12  1 

"■)  CL  *  CL  ••o*J 

32  2 


(9-5) 


It  turns  out  that  If  we  keep  terms  out  to  a^,  we  shall 
then  have  all  the  terms  containing  t3,  t^,  t^,  . 


■  :.f  ■;■>  t  cfc  ;.V  3  /  3  Lf  ■ "1  Offf 


'  '-3  33 


o' 


. . . •• 


.  ft‘.  33!i  ,  :  3.3, 


f 


,  n . :  ■  h  o'  a  one-  i  ,  ,:..e  ©v.s.i  3,3(1 
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With  the  substitution  of  (9-3)  -  (9*5)*  equation  (8.18) 
becomes 


1 

(  I-Icdt)  ^ 


[2  a4  -  3  a6  +  6a8  -  16  al°  +  4lcuT  a* 2 


6icoT  a4  +  12io:T  a6  _  32icoT  a8 
(l-ioox)  (I-Iodt)  (1-Iodt) 


44l  Icdt  „  10 1  ,  1  r  2  cjo2t2 

- -  - - -  a  J  +  - yr  L - 5 

4  ( 1-  Icdt)  ( l-  iooT)  d  ( 1-  iayr)  d 


12  icDT  „ 2  8  cd2t2  „ 2  15  .4  1  „ 2 

- 7T  a  -  - ^  a  -  — -  a  -  - ~  a 

( 1-  1<jdt)  2  ( 1-  Icdt)  2  2  ( 1-  Icdt)  2 


2  oo2t2  a2  +  12  Icdt  a4  +  8  032t2  a4 

(1-Icdt)2  (1-1cot)2  (1-Icdt)2 


+  a6  +  .  1  a4  +  1  OLJlI  .  21  jSL  a6 

(1-icDx)2  2  (I-Igot)2  (  1-Icdt)  2 


14  cd2t2  „6 

"  ci 

(  1-Icot)  2 


105  a8 


1  1  a6 

4  (l-ioox)2 


8 


.  ;  f  '  .i  ,  '  .  r:  :  i! 


.  rikJL 
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9  <d2t2  a6  +  94  Icdt  a8  +  36  oo2x2  a8  +  135  a10 
(l-icox)2  (1-Icjdx)2  (l-ioox)2  4 


9  1  a8  +  249  oo2x2  a8  _  747  icax  &10 

2  (l-icox)2  8  (1-icox)2  4  (1-icox)2 


249  co2x2  10  249  1  10 

2  ( 1-icox)  2  16  ( 1-icox) 2 


1095  co2t2  a10  -j  +  1  [_13_lH_  a2 

( 1-  icox)  2  ( 1-  iojx)  ( 1-  icox)  2 


2  icox  +  4  cd2x2  &2  +  _17  a4  _  2  icox 

(l-icox)  8  (l-icox)2  2  (l-icox) 


2  „  2  13  icox  „4  icox  _  2 

’p  ™  81  ^  q  81 

(l-icox;  2  (l-icox)^  (1-icox)-5 


2  oo2x2  4  _  17  6  ,  icox  2  _  1  „4 

^  G<  U  1  Cl 

(l-icox)d  4  (l-icox)  (l-icox)2 


39  Icdx  o6,3  iQ3x  4  3  co2t2  6 

4  (1-iODT)2  2  (l-iCDT)3  (1-itUT)2 
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'I* 


3  Icdt  4  ,  3  1  a6 

O  c'*  •  o  O  ^ 

^  (l-iODT)  (I-Icdt)11" 


195  Icdt  a8  _  15  Icdt  &6  _  _15  (d2t2  a8 

8  (l- Icdt)  2  4  (1-1o)t)3  2  (1-Icdt)^ 

_  255  a10  +  15  ia^T  a6  _  15  1  a8 

16  4  (i-Icdt)  4  (I-Icdt)2 

,  1385  Icdt  „  10  ,  105  Icdt  _8 

*r  q  a.  "r  ■  q  9. 

16  (1-1cdt)^  8  (1-Iojt)^ 

105  cd2t2  10  _  JL05  Icdt  8 

4  (1-ioDT)2  8  ( 1—  icdt) 

+  125  1  10  _  945  Icdt  10 

8  ( 1-  icdt)  2  16  (1-1cdt)^ 

+  945  Icdt  &10  j  _  5  Icdt  a2 

16  ( 1—  Icdt)  8  (1-1cdt) 


-3  a4- 


(  1-1cdt) 


a2  + 


2  2 

2  CD  T 

(  1-Icdt)  ^ 


0 


(9.6) 
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Collecting  terms,  we  find  that  the  coefficients 
o  2 

of  a  and  a  are  identically  zero.  Thus,  by  keeping 
terms  out  to  a  ^  ,  we  actually  end  up  only  with  terms 
out  to  a^.  Carrying  out  this  programme,  (9.6)  reduces  to 

a^[-  5  -  64  ioDT  +  327  o)2t2  -  810  ior^] 

+  a^[20  icDT  +  24  uX  rX  +  216  iou^T^] 

+  a^[20  aXiX  -  88  iao^T^]  -  24  =  0  (9.7) 


8  10 

The  reason  that  we  do  not  keep  a  ,  a-  ,  . . . 
is  because  these  higher  order  terms  contain  or 

higher.  We  can  now  go  through  the  same  sort  of 
procedure  as  in  Chapter  7.  Rewriting  (9.7)* 


,666 

k  x  T 

(  1-Icdt) 


7> 


-  ^  -  16  iCDT  ] 


4  4  4 

+  ..X  .7  ....  [  5  io>T  +  6  co2t2  +  54  ico^r^  ] 

( 1-Icdt)  ^ 

+  —  -x-  -T  -p  [  5  oXX  -  22  ioo^T^  ]  -•  6  =  0 

(  1-Icdt) 


(9.8) 


& 


:  '  -- 
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Dividing  through  by  and  then  multiplying 

(9.8)  by  (l-icDT)C\  we  obtain  (remembering  that  at  any 
stage  we  keep  terms  only  out  to  t^) 

k^x^x2  [  -  -  16  io)T  ]  +  k^x^  [  5  iooT  +  16  cd2t2 

+  37  ico^x^  ]  +  k^x^a)2  [5-42  icor  -  118  cjo2t2 

+  152  iaPx^  ]  -  6  ao^  +  36  iar^T  +  90  co^t2 
-  120  ia7  t3  =  0  (9-9) 


Assuming  that,  for  cut  «  1,  k^  and  k2  will 
be  of  the  same  form  as  previously  obtained,  we  have 


-  5.  [  k^8  +  6  ik-^5k2  ]  x^t^ 

-  16  io^T^k^x^  +  5  ioor  [  +  4  ik2_3k2 

-  6  k12  k22  ]  x^  +  16  oo2t2  [  k-^  +  4  ik^kg  ] 
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+  37  io)3T3k14x4  +  5  o?  [  kq2  +  2  ik1k2  -  k22  ]  x2 


-  42  1cd3t  [  kq2  +  2  ik1k2  -  k2^  ] 


2  i  2 

1  x 


-  118  go4t2  [  kn  2  +  2  lkxk2  ]  x2  +  152  ia£T3k-,2x2 


-  6  CD4  +  36  Icd^t  +  90  cd3t2  -  120  i(J03T3  = 


0 


(9.10) 


Equation  (9. 10)  may  be  easily  solved  to  the 
first  order  in  z,  yielding  the  expected  results 


T  =  V  -  (2f°)l/2 

0  3m 


(7.9) 


a  =  ^ 


3  co2t 

5 


(7.12) 


To  the  second  order  in  z3  we  expect 


a  = 


3  cd2t 
5  Vo 


(7.12) 


V  [  1  +  br2  ]1//2 
o 


V 


(7.17) 


94. 


This,  indeed,  turns  out  to  be  true.  Solving 

p 

the  second  order  equation  just  yields  b  =  1.72  of~ 
again.  To  the  third  order  in  t,  the  real  part  of 
(9.10)  leads  to  (7.17)*  Using  this  result,  the 
imaginary  part  becomes,  writing 

a  =  3  [1  +  Bt2]  , 

^  v  o 


15  .  aP  .  3  cc2t  .  216  v  6t2 
2  Vq5  5  VQ  125  0 


o 

16  CDTu 


216  6 
- -  GO 

125 


+  5  COT 


1-3.44  co2t2  ]  •  ||  Vo4 


-  30  CDT 


_9  co4fg  36  y  4 

25  VQ2  25  0 


+  64  CD2  T2 


0)3  _  3  (Jo8  T  .  _3§  y 

Vq3  5  Vq  25  ( 


+  37  CO-^  T' 


3^3  .  w 


4 


V 


"4 


^  +  10  CD2  •  — 

25  0  VG 


1  -  0.86  gd2t2] 


-  V  J  '  '■  :  "  ■  '-)/  ■■ 


TI,  <■.  •  '  :  .  !?  ; 
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■3  2 

3  0)  T 
5  V0 


1  +  Btc 


]  ■  -g  V  4 

5  o 


-  42 


3 

CO  T 


o 


[  1  -  1.72  0)2  t2 


6 

5 


V, 


2 


+  42  or*T 


_9 

25 


4  2 

0)  T 


V, 


o 


-  236  CO^T2 


CO 

V. 


3  0)2  T  6  _  2 

5  VG  *  5  o 


+  152  CO^T^  *  CO2  +  36  0)5  T  -  120  CO^T^  = 


0 


This  last  detailed  equation  has  just  been  included  to 
give  the  reader  some  idea  of  how  we  actually  carry  out 
the  calculations .  Solving  for  B  ,  we  have 


~  B  [  972  +  3456  +  3096  +  1944  -  6912 

5  125 


-  6660  +  774  -  10,836  -  2268  +  21,240 


-  22,800  +  15,000  ] 
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or 


B 


2994  2 

-  CD 
900 


-3-33  (B2 


(9-11) 


which  agrees  with  our  previous  low  frequency  results. 

Thus,  expanding  the  general  dispersion  relation 
for  |  a  |  «  1  and  Im  (-i)  >  0  yields 

cl 


a  = 


3 

5  V 

o 


[  1 


3.33  032T2  ] 


(7.21) 


v  =  Vo  [  1  +  0.86  cd2t2  ] 


(7.18) 


for  CDT  «  1 . 

This  gives  us  a  partial  check  on  the  validity 
of  the  general  dispersion  relation.  That  is  to  say, 
we  are  now  reasonably  sure  that  our  algebra  is  correct. 

What  about  the  two  assumptions  that  we  have 
made  in  expanding  Q  ? 

(i)  |  a  |  «  1 


Now, 


kxr  (k]_+ik2)xT 


(9-12) 


(  1  —  iCDT  ) 


(  1-iCDT) 


-  1  v'O i.  ^./oive^q  dw  i  v-  a99«x&*  ilot  % 


I  >>  ro)  f: 


■  (  -:  f  \.j  ,  n.  ■■  j‘  ;  --jj  "  '■  •  .■  J  « .Jt/ocfB  sr w 


■ 


1  )I 
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Since  we  are  considering  cdt  «  1,  we  take  only 
the  first  order  results 


k- 


co 

V 


CD 

V, 


k. 


3 

5 


2 

CD  T 

V 


Then,  for  cdt  «  1, 


(7-9) 

(7.12) 


x(fl  + 


i  *  ¥ 


a  = 


o 


|  V-)d + 


(l  +  cd2t2) 


r  /CDT  3  cd3t3n  .  /  CD2  r 

UV  "5  V  '  [  V 

o  J  o  o 


+  i 


3 

5 


2  2 
CD  T 

V 

o 


)] 


or 


a  = 


(f)1/2  [(cot  -  |  coV)  + 


8  2  2 

5  ~ 


^  ]  (9.13) 


From  (9.13)*  we  see  that  for  cdt  «  1  ( ie .  for  cd  ->  0) 

|  a  |  «  1  confirming  our  first  assumption. 
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(ii) 


We  also  assumed  that 


Ira  (|)  >  0. 


Now, 


Ira  (i) 


=  + 


[kj_  -  kgcDx] 

(k12  +  k22)xT 


cu 

V 


3  03-*  T 


3^3 


o 


V 


o 


■0)e 


4^2 


V, 


9  a)  r 

25  V, 


■]XT 


V 

_ o 

XCDT 


[1 


3  2  2  9  2  2n 

5  ^  T  "  25  ^  T  ] 


<?>1/8 


1  r  24  2  2-t  /  n  n  i.  \ 

—  [1  -  25  CD  T  ]  (9.14) 


As  the  frequency  approaches  zero,  Im  (-^)  +  oo . 

ie .  for  odt  «  1,  Im  (— )  >0  as  assumed.  Therefore, 

a 

the  asymptotic  expansion  that  we  have  used  for  Z(-^) 
to  obtain  the  low-frequency  sound  behaviour  is  correct. 


...  :. 


'  V 
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Chapter  10.  HIGH  FREQUENCY  RESULTS:  orr  »  1 

The  difficulty  in  solving  the  general  dispersion 
relation  (8.l8)  for  gqt  »  1  is  that  no  satisfactory 
experimental  data  exist  for  this  region  which  would 
suggest  the  correct  expansion  to  take  for  Qq .  That  is 
to  say,  we  do  not  initially  know  whether  |  a  |  is  greater 
than  or  less  than  unity  and  whether  Im  (4-)  ^  0  or  is 
less  than  zero.  Without  experimental  data  to  suggest 
what  we  should  expect  for  k^  and  k^  at  very  high 
frequencies,  all  that  we  can  do  is  make  calculations 
based  upon  certain  assumptions  and  then  check  that  the 
final  results  are  consistent  with  these  assumptions. 

Thus,  we  seek  out  self-consistent  high-frequency  solutions 
to  the  problem.  Before  outlining  the  approach  which 
gives  the  only  self-consistent  solution  that  has  been 
found,  we  shall  very  briefly  consider  several  attempts 
which  do  not  yield  physically  acceptable  or  self-consistent 
high  frequency  solutions . 

10.1  First  Approach:  |  a  |  >  1;  Im  (~)  ^  0 

3. 


In  our  first  attempt,  we  assume  that  |  a  |  >  1 

for  cdt  »  1.  For  |  a  |  >  1,  we  have,,  from  Appendix  III, 


,  '  : .  ■  t  •  ■ 

' 
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that. 


z(0  =  iTT 


1/2  e-C‘ 


-  2?[1  -  |C2 


+ 


_4 

15 


c4  ...] 


(III.18) 


In  particular,  for  £ 


1 

a  J 


Z(|) 


Itt 


i/2  _  21 
a 


+ 


i^/2 

a2 


4  i  ,  1  iTr1/2 
3  P+2 


(10.1) 


where  we  have  expanded  the  exponential 
to  the  first  order  in  l/a. 


Then, 


Q 


o 


a(  1-Igdt) 


(10.2) 


2  S 

However,  Qq  and  Qq  contain  only  higher  orders  In 
l/a.  The  general  dispersion  relation  (8.18)  may  be 
written  as 


Qo3  [  2  +  4  lcuT 


( l-iaox) 


2od2t2  #  1 
(I-Icdt)2  s? 


12  iooT  .  1 _ 8  od2t2  1 _ 15 

(1-Icdt)2  a2  (l-icox)2  a2  2 
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1  .  _1_ 

(1-Icdt)2 


+  Q„ 


13  Icdt  <  _1_ 
(1-Icdt)2  a2 


2  lovr 

(  1-Icdt)  3 


1  2  2 

1  +  Z  ^  T 

£?  (1-Icdt)2 


+ 


17 

2 


2  Icdt 
(  1-Icdt) 


2 

(  1-Icdt)  2 


5  lg)T 

(  1-Icdt)  2 


3  - 


1 

(  I-Icdt)  2 


1_  +  2  ^2f2  .  1 

a2  (1-Icdt)3  a5 


0 


(10.3) 


Then,  using  (10.2),  equation  ( 10 . 3)  yields, 
to  the  first  order  in  l/a. 


17 

2 


a(  I-Icdt) 


3  =  0 


(10.4) 


or 


kXT  =  ¥  T1/2 

6 


(10.5) 
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so  that 

klXT  =  "T 

1/2 
7T  / 

or 

> 

II 

,1s 

7 r1/2 

6 

COT 

(10.6) 

Does  (10.6)  satisfy  the  initial  assumptions  that  were  made? 
In  particular ,  what  about  |  a  |  ?  For  cot  »  1* 

[kj2  +  kp2]1/2  xt 
|  a  |  =  — ± - 2 - 

COT 


and  using  (10.6),  this  becomes 


1/2 


a  |  = 

6  COT 


(10.7) 


which  will  yield  |  a  |  <  1  for  cot  »  1  contrary  to 

our  original  assumption.  Thus,  for  cot  arbitrarily 
large,  this  solution  is  not  consistent  and  must  be 
rejected . 


However,  for  cot  ^  1,  we  have  an  interesting 
result.  For  this  range  of  frequencies,  |  a  |  >  1, 

so  that  the  low-frequency  solution  is  consistent.  That 
is  to  say,  we  obtain  a  solution  (namely  (10.6))  that 


' 

b..ioe.  •'•j 
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does  not  display  the  usual  low  frequency  sound-like 
behaviour.  Perhaps,  this  solution  corresponds  to  a 
heat  wave  or  more  likely  it  is  an  extraneous  solution  due 
to  the  crudity  of  the  calculation  carried  out  here.  But, 
at  least,  we  have  shown  that  the  power  series 
expansion  for  Qq  ,  corresponding  to  j  a  j  >  1,  does  not 
yield  any  consistent  solution  for  cot  »  1. 

10 . 2  Second  Approach:  |  a  |  <  1  ;  Im  (i/a)  >  0 

In  our  second  attempt,  we  make  the  following 
assumptions:  |  a  |  <  1  and  Im  (i/a)  >  0  for  <dt  »  1. 

For  this  situation,  we  can  consider  the  first 
few  terms  of  equation  (9-7) 

a2  [  20  cd2/2  -  88  io/x^  ]  -  24  co^x^  =  0 

or  a2  [  5  ~  22  icox  ]  -  6  co2x2  =  0  (10.8) 

In  arriving  at  (10.8),  we  have  really  said  nothing  yet 


about  cox. 


. 
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For  cot  »  1,  (10.8)  gives 


a 


+  /T  (tax)  1/2  r  !  .  _5  _i 
A-''  11  (_i)  V2  44  CDT 


(10.9) 


but 


,  so  that 


a 


1  +  i 


5  _1  ,  _5  .  _1  i 

"44  cot  44  cot  J 


(10.10) 


,  kXT  kXT  t  -\  i  -i  \ 

where  a  =  -  =  -  (10.11) 

(1-icoT)  (-icoT)(l  +  i/cor) 


Then,  on  substituting  (10. 11)  into  (10. 10),  we  obtain 
(for  cot  »  l) 

kXT  =  (cot)3/2  [i  -1]  (10.12) 

=  (k1  +  ik2)  xt 


Taking  the  +  sign,  we  have 


(k1  +  ik2)  x  t 


(cot)3//2  [i  -  l] 


(10.13) 


t  •  : 
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which  leads  to 


k2xx 


k-.  XT 


(o>t)V2 

(mx)3/2 


(10.14) 


(10.15) 


This  solution  is  physically  unacceptable 
because  relation  (10.15)  yields  a  negative  phase  velocity, 
viz . 


v  -  -M  7^72  (1°- 

So,  we  reject  the  solution  corresponding  to  taking  the 
+  sign  in  equation  (10.12). 

Now,  taking  the  minus  sign,  we  have 


(k1  +  lkg)xx  =  -/2I  (tor)3/2  [1-1]  (10.17) 

which  yields 

k2xx  =  - 22  (u^)3/2  (10.18) 


k1XT  =  (cdt)3/2 


(10.19) 
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We  cannot  accept  this  solution  because  it  gives 
us  amplification  of  the  sound  wave,  ie . ,  a  <  o.  A 
further  reason  for  rejecting  this  second  solution  is 
that 


|  a  |  =  (cot)  1//2  >1  for  cjdt  »  1 

contrary  to  our  original  assumption.  This  process  can  be 
carried  out  to  higher  orders  in  a,  but  these  also  yield 
inconsistent  high-frequency  solutions. 

Thus,  once  again  we  have  not  been  able  to  obtain 
any  self-consistent  or  physically  acceptable  high-frequency 
solution.  The  only  possible  procedure  left  to  us  now  is 
to  consider  the  asymptotic  expansion  of  Z(i/a)  with  the 
exponential  term  included.  There  are  two  possibilities 
to  be  considered,  namely  Im  (i/a)  <  0  and  Im  (i/a)  =  0. 

We  shall  consider  the  former  possibility  next. 

10 . 3  Third  Approach:  |  a  |  <  1  j  Im  (i/a)  <  0 

For  cjdt  »  1,  we  assume  that  |  a  |  <  1  and 
Im  (i/a)  <  0.  Equation  (8.21)  gives  the  relation  between 
Q0  and  Z(£),  namely 
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Qr 


1 


°o  ia(  1-Igdt) 

where  £  =  l(1^T)  =  1/ 


z(c) 

(8.21) 

1  /a. 

(II. 2) 

Then,  from  Appendix  III,  the  asymptotic  expansion  of 
Q0  for  Im  (l/a)  <0  is 


Q, 


2t t 


1/2  e1/^  1 

G  ± -  [  1 


(l-icoT)a  (I-Icdt) 


1  fi2  ,  3  a4 
~2  a  ^  Tj;  a 


^  a®  +  ...  ] 


(10.20) 


Let, 


y  = 


1/2  1/a2 

7 t  '  e  ' 


a 


(10.21) 


so  that 


Q, 


_2y _ + _ i_ 


(1-1o3t)  (1-Icjot) 


1  -  a2  +  ^  a4  -  a6  + 


(10.22) 


Q 


4y  +  %  [  1 


(I-Icdt)2  (1-Icdt)2 


■|  a2  +  ^  a21 


a8  +  .  .  .  ]  +  - — - p  [  1  -  a2  +  X  a4 

8  (1-Iodt)2  4 
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-  |  a6  +  ...] 


(10.23) 


Q, 


3  _  By3 


_  +  -.iJLzi.  [  i 

(l-icor)^  (1-Icdt)^ 


a2  +  |  a4 


15  a6  +  ...  ]  +  - [  1  -  a2 


~8 


(  1-Icot) 


3 


+  l  a4  '  |  a6  + 


]  + 


( l-ioor)  3 


[  1 


|  a2  +  3  a4  -  8  a°  +  ...  ] 


A  o  _6 


(10.24) 


Substituting  (10.22)  -  (10.24)  Into  the  general 
dispersion  equation  (8.l8)  yields 


By 


3 


(  1-Icot)  ^ 


.  [  2  a4  +  a 2  ]  +  ..  -lg.-£  [  1  1  a2 

(l-iarr)  (I-Icdt)^ 


3  y  .  15  a6  ]  .  [  2  a4  +  .  1  a2  ]  +  - hi 


(I-Icdt)  (1-1cjot)3 


-[  1  -  a2  +  |  a4  -  |  a6  ]  •  [  2  a^  + 


4  '4  ia)T  a2 


( l-i<DT) 


109. 


+ 


- - - ~  [  1  -  i  a2  +  3  -  8  a  6  ]  •  [  2  a^ 

(  1-Icdt)  3 


+ 


4  lcoT  a2  ]  +  f — +  - LX_  [  i  1  a2 


(  1-Icdt) 


(  1-Icdt)  2  (  1-Icdt)  2 


3  _4  15-6 

71  8  '  '  (1-103T)2 


+  tt  a~  |  +  - - - -  [  1  -  a2  +  ^  a^  -  |  a6  ]  | 


2  12  Icdt  „2  8  cd2t2  „2  15 

cd  ™ *  _  d.  “*  cl 


(  1-Icdt)  2  (  1-Icdt)  2  (  1-Icdt)  ^ 


(  1-Icdt) 


a2  /  + 


-2-Z-  +  - [  1  -  |  a2 


(  1-iCDT)  (  1-iCDT) 


,  3  A  15  a6  1 
+  7f  a 5  J 


13  Iidt  a2  +  2  icDT 

(  1-Icdt)  2  (  1-Icdt)  ^ 


+  4  a2  +  II  a4  _  2  Ict  +  - 2  a2 

(1-1(Dt)2  2  (  1  —  ICDT )  (  1-Icdt) 


5  luJT  a2  -  3  a4 


(  1-Icdt)  2 


2  i  2  (d2t2 

a  + 


(1-1cdt)2  (  1-Icdt)  3 


=  0 


(10.25) 
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Multiply  (10.25)  by  (I-Icdt)^  and  collect  terms, 
keeping  only  terms  out  to  a^.  We  obtain 

y3  [  16  a^  -  16  ioDT  a^  +  32  Iodt  a2  ] 

+  y2  [  8  cd2t2  -  4  a2  -  32  cu2T2a2  -  6  a^ 

+  12  iooT  a^  +  30  CD2x^a^  -  12  a^  +  48  Icot  bP  ] 
+  y  [  4  iaA3  -  2  Icdt  a^  -  2  o)2T^a^  -  8  lo:3T3a^ 

+  a^  -  3  Icot  +  <jD2x2a^  +  17  icc>3T3a^ 

-  12  icor  a^  -  54  o^T^a^  ] 

-  3  co^T^a^  +  aj2T2a^  -  11  ico3T3a^  =  0 


(10.26) 

The  cubic  equation  (10.26)  has  been  solved  for  y 
(see  Appendix  VII),  yielding  the  three  solutions 

yi  =  i  la)T  *  72 


(VII. 32) 
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(VII. 36) 


2  io)T 


(VII. 34) 


In  ( VII . 32, 36, 3*0  we  have  only  written  down  the  first 
non-vanishing  term  on  the  right-hand  side.  If  one  solves 
the  cubic  equation  (10.26)  to  higher  orders,  he  obtains 
a  series  of  the  general  form 


(  )  a°  +  (  )  a2  +  (  )  a^  +  .  .  . 


y 


for  each  solution.  Since  we  have  assumed  that  |  a  |  <1 

(and  so  far  we  have  said  nothing  about  cot)  the  higher- 
order  terms  will  generally  be  progressively  smaller. 

In  any  case,  if  we  kept  the  higher-order  terms  on  the 
right-hand  side  of  ( VII . 32, 36, 3^) *  the  difficulty  in 
solving  the  equations  for  k-^,  would  become  much 

greater.  Even  with  just  the  leading  term  the  calculations 
are  somewhat  messy. 

Let  us  consider  each  solution  y^  separately, 
starting  with  yq . 


From  (VII. 32) 
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1/2  l/a2  ,  , 

7 t  '  e  '  1  .  1 

- - -  -  ^  10)T  *  — 

&  Q 


(10.27) 


or 


4  7r1//2  i 


=  1  e-Va‘ 


COT 


where,  of  course,  a  = 


kXT 


(  1-1cot) 


as  before 


So, 


4  7 r1/2  1  _  (  1-1cot)  -  (  1-1cot)  2/k2x2T2 


COT 


kXT 


i_ 

COT'' 

k. 


-  1cot(  1  +  — )  .  \  2  a  2  2  2 

v  /I''r/  -(1-icot)  /k  x  t 


kn  (1  +  i  ^-)xt 


(10.28) 


Now,  define  r  =  k9/k 


2 7  “1 


(10.29) 


Also, 


(  1-1cot)  ‘ 

2  2  2 

k  X^T 


(1-2  Icot  -  co2t2) 

k12(l+lr)2  x2t2 


[(1  -  co2 t2 )  -  2  Icot]  [l  -  2  Ir  -  r2] 

k^2(l  +  r2)2  x2t2 


r  (co2t2  -  l)(l  -  r2)  +  4  rcoT  i 
k-^2(l+r2)2  x2t2 
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1[  -  2  cdt(  1  -  r^)  -  2r(l  - 


2  2 
Qj  T+ 


n 


k-i  2  ( 1  +  r 


2)2  x2t2 


cp  -  10 


(10.30) 


In  particular,  for  cot  »  1  and  writing  =  co/V 


V* 


<P  =  g,'  g,2  [  d-*2)  +  -  ] 

x  (l+r  )  cot 


(10. 31) 


and. 


0  = 


V£ 


x2( l+r2) 


[  -2(l~r2)  +  2r  ] 


(10.32) 


(JOT 


For  arbitrary  got,  (10.28)  may  be  written  as 


4  7T 


1/2 


GOT 


^ 

(i  +  ) 

CD  T 


V  llzlri  ecp-10  (10_33) 
x  (l+r2) 


or 


V  e 


4  ++(1  -  i.) 

_ _ GDTy  _  _  _ 

cdt(1  +  )  x  (l+r2) 

GC+T+ 


[  (cos0  -  r  sin0  ) 


-  i(r  cos  0  +  sin  0)]  (10.34) 


On  separating  real  and  imaginary  parts,  we  have  two 


equations  in  the  two  unknowns  V  and 

CDT. 


4 


7 r 


1/2 


& 


cot(  1  + 


CD*-  T 


V  _ 

x  (l+rh 


[  cos  0  -  i 


4  7 r 

— - 

CD  T  (l  + 


1/2 


2  2' 
CD  T 


V  e* 
x  (l+r2) 


Dividing  (10.35)  by  (10.36)  gives 


cos  0  -  r  sin  0  1  -  r  tan 

cdt  =  -  =  - 

r  cos  0  +  sin  0  r  +  tan  0 


or 


tan  6  =  1  -  raiT 

r  +  CDT 


or 


1  -  opt  tan  0 
cdt  +  tan  0 


as  functions  of 


sin  e]  (10.35) 


+  sin  0]  (10.36) 


(10.37) 


(10.38) 


(10.39) 


Prom  (10.30)  and  (10-38)  we  have 


) 


■ 
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0  = 


V‘ 


x2(l  +  r2)2 


[  2r(l  -  -5^) 


1a  _  2(1  -  r2) 


2  2 

CD  T  COT 


=  tan'1  (1^22) 
r+CDT 


(10.40) 


or  rewriting  (10.40) 


V‘ 


(i  +  ^2)2  tan'2  4igr> 

[2r(l  -  -J^s)  -  lilzjfl) 


(10.41) 


2  2 
Cj 


CUT 


For  (l0.4l)  to  have  any  physical  meanings  we  must 
insist  that  V2/x2  >  0. 


Now,  multiplying  (10.35)  and  (10. 36) 


together 


_ 16  IT 

oPt^  (l  + 

0)  T 


V2  e2cp 


x2  (l+r2)2 


[  r  cos2  0 


-  r2  sin  0  cos  0  +  sin  0  cos  0 


p 

-  r  sin  0  J 


(10.42) 
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But., 


cos2  0  -  sin^  0  =  cos  20 

sin  9  cos  0  =  75  sin  20 


so  that  (10.42)  becomes 


_ 16  7 r _ 

(l  +  -^-p) 2 

(XT 


y2  e2cp 

x2  (l+r2)2 


[r  cos  20  +  -|(l-r2)  sin  20] 

(10.43) 


We  can  eliminate  V2/x2  from  (10.43)  by  using  (10. 4l). 
Thus,  we  finally  have 


l6  7 r 


3^3 


CJO  T 


(1  + 


032T2 


tan’1^)  exp 


2  tan-^^JCd  -^)(l-r2)  + 


m2  ■ 


[2r(l  -  1  )  - 

(XT  ODT 


r  eos  [2  tan-1(  j;~™T)  ]  +3  U-r2)  sin  [2  tan~1(h™T)l 


r+o)T 


(10.44) 
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We  have  obtained  one  rather  complicated 
equation  in  one  unknown  .( 10 . 44)  gives  r  as  a  function 
of  cdt.  We  shall  attempt  to  solve  (10.44)  for  r  for 
cdt  »  1.  Once  we  have  determined  r  as  a  function 
of  cdt  ,  the  other  unknown,  V,  can  be  readily  found* 
from  (l0.4l).  For  cdt  »  1,  there  are  three  possible 
ranges  of  values  for  r  ,  namely,  r  —  1  ;  r  «  1  ; 
r  »  1. 

(a)  Consider  the  first  possibility,  r  —  1  for 
cdt>>1  .  Then  in  (10.44), 

1-pcdt  ~  _  r 

r+CDT 

for  r  ~  1  and  cdt  »  1.  Using  this  result,  equation 
(10.44)  reduces  to  (for  cdt  »  1) 

16  7r  =  tan~1(-r)  egtan~4-r  H  l-rg)/2r 

cd^  t^  2r 

^r  cos[2  tan”^(-r)]  +  -|(l-r2)  sin[2  tan-'*' ( -r )  ]  J 


(10.45) 


■ 


r 


...  '  ■■  ,  '  v  ,  :  , ‘  ;  o!, 


i..,  su:  '■  .  .:  Itja  ■•£ 
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We  throw  away  all  explicit  cot  dependence  on  the  right- 
hand  side  because  It  turns  out  that  it  is  not  important. 
Now,  expand  tan-1(-r)  about  r  =  +1  and  call 
r  -  1  =  <5  (which  is  a  function  of  cot  in  general). 

Of  course,  6  «  1. 

tan-1 ( -r )  =  ^  ~  ^  6  +  ^  ^  •  (10.46) 

The  multi-valuedness  of  the  inverse  tangent  will  be 
discussed  in  detail  shortly. 

Then, 


cos[2  tan_1(-r)]  =  cos  [-—  -  (6  -  S2)  ] 

=  cos  {■—)  cos  (6  -  62)  +  sin  sin  (6  -  52) 

=  -  1  (6  -  62)  =  -  6  +  &2  (10.47) 


Similarly, 

sin  [2  tan-'L(-r)]  =  sin  (~)  cop  (6  -  6^) 

-  cos  (~)  sin-  (6  -  62) 
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=  -  1  (1  -  \  <52) 

=  -  1  +  62  (10.48) 


Then,,  putting  r  -  1  =  6  and  using  (10.47)  and  (10.48) 

the  right-hand  side  of  (10.45)  becomes 


6  +  62] 
2(1  +  6) 


[(1  +  6)(-  6+  62) 

+  -|(l  -  1  -  26  -  62)(-  1  +  -|  62)] 


6  +  62][-  26 


62]/2(l+6) } 


35  ,  1  s2  _  3ir  .2 
5T  2  6  -  15  6 


(10.49) 


Equating  this  to  the  left-hand  side  of  (10.45)  yields 


l67T 

3  3 
or 


3i  &2 
1 6 


256 


3,3 


3  O)  T' 


(10.50) 


6 


+ 


(10.51) 


or. 


_ 1 6 

/T~ (cot)  V2 


We  choose  the  minus  sign  because  we  expect  from  examining 
the  limited  experimental  data  that  kg/k^  ^  1.  Thus, 

r  =  1+6=1-  16  (10.52) 

J  3  (ojt)3/2 

As  cut  becomes  indefinitely  large,  r  -*■  +  1 .  This 

result  is  also  easily  seen  directly  from  (10.44).  As 

an  example,  for  cot  =  100,  r  =  0.991  and  for 

cot  =  1000,  r  =  0.9991*  From  (10.52),  we  see  that 

1 — rcoT 

for  cot  =  100  the  error  in  replacing  ■  by  -r 

is  less  than  2%.  It  is  also  to  be  noted  (from  (10.52)) 
that  the  approximate  condition  cot  »  1  should  be 
replaced  by  the  more  accurate  limitation  cot  »  4.4. 

Now,  taking  cot  sufficiently  large  (cot  »  4.4) 
so  that  r  =  1  to  a  good  approximation,  equation 


(l0.4l)  yields 
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x/V^  =0.46  (10.54) 

or  VQ/V  =0.42  (10.55) 

Since  r  -*■  +  1  for  cut  very  large,  we  also  have 

av 

— -  =  0.42  (10.56) 

cu 

It  should  now  be  noted  that,  in  arriving  at  (10.53) *  we 
took  tan-1(-l)  =  3^/4.  Mathematically  speaking,  there 

is  no  reason  to  exclude  the  other  angles  which  satisfy 
the  multi-valued  inverse  tangent.  Since  we  physically 

p  p 

expect  V  /x  >  0,  we  cannot  consider  negative  angles 
such  as  -tt/4  for  tan”^(-l).  However,  there  is  no 
reason  for  not  considering  all  positive  angles 

_  1 

corresponding  to  tan-  (-1),  such  as  7^/4,  11tt/4, 

etc.  That  is  to  say,  in  general,  (10.53)  could  be 
replaced  by 


(2n  +  3)  | 


(10.53)  (b) 


where  n  =  0,  2,  4,  6  ....  . 

n  =  0  just  gives  the  "principal"  value  (10.53). 
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Whether  these  higher-order  solutions,  n  =  2,  4,  6  ...  , 
have  any  physical  meaning  is  not  clear.  One  might  expect 
a  discrete  spectrum  of  frequencies  (see  ref.  1,  p.  96) 
but  what  does  it  mean,  if  anything,  to  have  a  spectrum 
of  possible  absorption  and  dispersion  curves?  One  possible 
argument  for  rejecting  n  =  2,  4  . . .  is  that  for  a  very 
large,  but  fixed,  frequency,  the  solutions  corresponding 
to  increasing  values  of  n  display  less  and  less 
absorption  and  more  and  more  (see  (10.57))  dispersion. 

From  the  examination  of  similar  physical  problems,  one 
would  expect  more  and  more  absorption  for  increasing  n. 

Unless  some  strong  argument  to  the  contrary 
sheds  illumination  on  this  question,  we  shall  reject 
the  higher-order  solutions  as  being  physically  meaningless. 
However,  we  should  always  bear  in  mind  that  we  may  be 
wrong  in  rejecting  these  solutions.  The  only  absolute 
way  of  confirming  or  disproving  our  ideas  on  this 
matter  is  for  some  experimentalist  to  provide  us  with 
pertinent  high-frequency  data. 

We  see  that  the  phase  velocity  V  and  oV0/cd 
approach  constant  values.  From  ( 10.53) j  we  have  that 
V  approaches  a  value  about  twice  VQ . 


Even  though  we 


' 


.  .  ■  i  ■  '  : 

' 


10 


. 

- 

' 
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cannot  legitimately  compare  these  results,  (10.55)  and 
(IO.56),  with  the  experimental  data  of  Meyer  and  |3essler 
(as  pointed  out  in  Chapter  2),  we  see  that  our  curves 
bear  a  remarkable  similarity  to  their  experimental 
curves  (see  Fig.  3)*  especially  for  the  velocity 
dispersion  curves.  Furthermore,  a  contains  no  x 
dependence  so  that  the  results  do  not  depend  at  all  on 
which  mono-atomic  gas  is  being  considered. 

Do  (10.55)  and  (10. 56)  satisfy  the  initial 
assumptions  (  |  a  |  <1;  Im  (i/a)  <  0)  that  were 
made?  For  ox  »  1,  (10.55)  and  (10. 56)  lead  to  the 
result  |  a  |  =0,65  which  is.  Indeed,  less  than  unity. 
Thus,  it  is  possible  to  Judge  the  approximate  error  in 
omitting  higher  order  terms  in  a  in  the  asymptotic 
series  for  y.  Each  term  in  the  series  (aside  from 
factors)  is  a2  greater  than  the  previous  one  and  is, 
therefore,  0.42  the  magnitude  of  the  previous  term. 

The  condition  Im  (i/a )  <0  is  also  satisfied 
since  £  =  i/a  =  1.09(l-i)  for  ot  00.  (The 

behaviour  of  the  pole  £  in  the  complex  plane  is 
sketched  in  Fig.  8.)  Thus,  (10.55)  and  (10. 56)  represent 
a  completely  consistent  high-frequency  solution.  It 
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should  be  pointed  out  in  passing  that  the  "rejected" 
solutions  also  satisfy  the  initial  assumptions  namely, 

|  a  |  <  1  and  Im  (i/a)  <  0. 

(b)  So  far  we  have  only  examined  r  ~  1  for  out  »1 
in  considering  equation  (10.27).  Does  (10.27)  have  any 
other  solutions?  Let  us  now  take  r  =  6 (out)  «  1 

for  out  »  1.  Then  (10.44)  becomes 


16 


7 r 


oPt3(1  +  -4-  -p) 

OU^ 


1  n2 


tan-1(  1~<^CUT)  exp 
6+out 


tan'+i^Hd  -  -ip)(l-S2)  +  ^1 

6+OUT  OU  T  tUT 

[26(1  -  -4o)  -  gh-6g)  ] 


...2  _2 
OU  T 


OUT 


[26(1  -  -i— )  -  2<1-6g)  ] 


2  2 
ouc  T^ 


OUT 


fs  cos  [2  tan 


-ljil&DT)]  +  1  (1.S2}  gln[2 

6+OUT  ^  O+OUT 


(10.57) 
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e  = 


,  -1  /1-5cot\ 

tan  (TOT} 


^n-1  (44^) 


tan-1  [  (—  -  6)(l  -  — )  ] 

vcot  '  v  COT' 


,  -1  r  1  s  6  ,  6 

lcdt  t2  cot 


and  expanding  about  zero 


0  =  —  -  6 - +  - - \  [—  -  6]  3 


1 

COT 


032T2  COT  3  COT 


1  .  6  I  .  +  I  63 


COT 


3  co^  t^  3 


(10.58) 


Since,  2  0  = 


—  -  26 

COT 


2  .  1  +  2  ,3 

3  co^  t^  3 


cos  20  =  1  - 


\  [  —  -  26  ]2 
2  COT 


or 


cos  20  =  1 


2  2 
co  T 


+ 


46 

COT 


26' 


(10.59) 
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sin  20  = 


2 

COT 


26-4 


1 

0)3  t3 


+  2  63 

3 


i(—  -  26) 3 
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OT 
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Then  (10.57)  becomes 


16  7 r 


O)3  T3 (  1  +  -7T-o)2 


CO2  T2 


6  -  I 


0)T 


3  0)3t3  3 


+  I  63] 


2[6  - 


0)2T2 


0)T  COT 


26 


exp 


COT 


^  +  2  a3][i  _  _1 
3  Q  3 


0)3  T3 


2  2 
O)^  T^ 


62  +  5 


+  M] 

2_2  cot 


0)  T 


2  [  6  - 


2  2 
0)  T 


i  +  ^  ] 

CJOT  COT 


[  6 - +  — -  263  +  —  -  6  - ^ — —  +  63 


2  2 
00  T  COT 


GOT 


3  3 
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or,  on  cancelling  terms 
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In  arriving  at  equation  (10.62),  it  has  been  assumed 
that  6  (cot)  «  1  for  cot  »  1.  Furthermore,  the 
functional  form  of  6  must  be  compatible  with 
|  a  |  <1  and  Im  (i/a)  <  0.  Also,  a  must  be  a  real 

positive  number  as  well  as  V. 


We  can  rewrite  (10.62)  as 


1 


It  is  seen  that  this  equation  is  satisfied  by 
6  =  1/cot  .  If  one  substitutes  this  into  ( 10 . 4l ) ,  viz. 


x 


V' 


.2 


2 


1/C0T 


and  takes  the  limit  for  cot  »  1  and  6 


i 
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(10.63) 


Clearly  this  cannot  be  accepted  since  it  would  imply 
that  the  phase  velocity  V  is  purely  imaginary.  No 
other  solution  of  (10.62)  can  be  found  of  the  form 
6(cdt)  «  1  for  o)T  »  1  which  is  compatible  with  all 
the  necessary  conditions,  so  we  must  conclude  that  we 
cannot  have  6(cdt)  «  1  for  cjdt  »  1. 

(c)  The  third  possibility  to  consider  is  <5(<jdt)  »  1 
for  cdt  »  1.  Once  again,  we  have 
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sin [2 


(10.44) 


If  6  increases  more  quickly  with  cdt  than 
does  (jot  Itself,  equation  (10.44)  becomes 


i6tt 

3  3 

or 


_  q  .  2  tan_1(-cDT)  [-62  +  ~  ] 

tan  (-cot)  exp  ' 


26  +  2  6^/cdt 


2  [(>  +  —] 
COT 


j^6  cos  [  2  tan-1  ( -got)  ]  -  62  sin[2  tan-1  (  -cot)  ]  J 


(10.64) 


and  for  cdt  >>  1, 

tan-1  (-cot)  —  —■ 

(cannot  consider  negative  angles  for  V2/x^  >  0),  so 
that 


3  1|  4  '  T  031 

64  ^  T  6 


6 


(10.65) 
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This  must  be  rejected  for  two  reasons:  (i)  the 
minus  sign  makes  6  negative  which  Is  Incompatible  with 
a  >  0,  V  >  0  .  (ii)  for  cot  »  1,  6  «  1  contrary 

to  the  initial  assumption  made  here. 

If  6(oot)  <  cdt  ,  (10.44)  yields  for  out  »  1 


16 


IT 


3  3 

CD  T 


tan-1 ( -  6) 


2tan‘1(_6).[_  62]/25 


2  6 


5  cos[2  tan_^(-6)]  -  ^  S2  sin[2  tan-^(-6)]j 


( 10 . 66) 


and  since  we  are  assuming  that  6  »  1 


64 


3cd^t^ 


(10.67) 


Because  of  the  minus  sign  in  ( 10.67).,  we  cannot 
expect  6  to  give  us  a  reasonable  result.  Therefore, 
for  cdt  »  1  we  cannot  have  any  solution  for  which 
6  »  1. 
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Thus,  the  only  consistent  and  physically 
acceptable  high-frequency  solution  of 


1  *  1 
yi  “  if  l£0T  ‘  7: 


(VII. 32) 


Is 


V  /V  =  0.42 

O'  oo 


(10.55) 


aV  /a>  =  0.42 


(10.56) 


What  about  y^  =  -  Icdt  a  and  y^  =  -  lo3T  ? 


(ii)  Consider, 


1/2  l/a2  . 

- - - -  =  -  —  loox 

a  2 


(10.68) 


We  proceed  in  exactly  the  same  way  as  we  did  with  y^ . 
Analogous  to  equation  (10.33) *  we  have 


^6!  (  1  +  .l  ) 


COT 


COT 


§  (  1  +  ir  )  e9'19  (10.69) 


where  cp  and  6  are  defined  by  (10.31)  and  (10.32) 
respectively.  The  real  and  imaginary  parts  of  (10.69) 
give 
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e1^  [  cos  0  +  r  sin  0  ] 


2l T 
CD 


—  [  r  cos  9  -  sin  9  ] 

V 


Dividing  (10.70)  by  (10. 71), 
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r  cos  9  -  sin  9 
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(10.72) 
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which  is  analogous  to  (I0.4l) 


Multiplying  (10. 70)  and  (10.71) »  we  easily 


obtain 
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[2r  ( 1  -  1  )  -  ] 
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exp 
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r  cos  [2  tan'^j1'-)  ]  +  -^(r2-l)  sin[2  tan"  ^  (-”-"--) 

r+toT  2  '  r+tox  ' 


(10.75) 


as  compared  with  (10.44) 


(10.75)  Is  solved  for  r  in  exactly  the  same  manner  as 
(10.44).  For  (jot  >2*  1  it  is  found  that  =  k.2 

again.  Then  (10.75)  has  only  solutions  of  the  type 
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OO 


(10.76) 
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(5/6)  V2 

(n  +  1)V2  .1/2 


(10.77) 


The  "principal"  value  n  =  0  does  not  satisfy  the 
condition  |  a  |  <1  since  |  a  |  =  1.14,  so  it  must 

he  rejected.  As  with  y^,  we  reject  the  higher-order 
solutions  n  =  2,  4,  6,  ...  as  being  physically  meaning¬ 
less  (although,  it  must  be  pointed  out  that  the  higher- 
order  solutions  do  satisfy  |  a  |  <  1  and  Im  (i/a)  <  0) . 


(iii)  For  the  third  solution  of  the  cubic,  namely 


y. 


1/2  1/a' 

t r  '  e  ' 


-  Icdt  a^ 


(10.78) 


so  that 


-hr  (  i  +  5i  )  = 

3cut  cut 
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(  1  +  ir  ) 


5  ep-i@ 


(10.79) 


Except  for  details,  the  calculation  here  is  exactly 
the  same  as  for  y^  and  y^  •  For  this  reason  we  shall 
only  indicate  the  results.  As  it  turns  out,  the  only 
solution  turns  out  to  be  exactly  the  same  as  (10.  j6)  and 
(10.77)  and,  therefore,  is  rejected. 
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Thus,  we  have  only  found  a  single  consistent 
and  physically  acceptable  solution,  namely  (10.55)  and 
(IO.56),  to  the  high-frequency  problem. 

id .4  1  a  |  <  l;  Im  (1/a)  =,  0 

To  complete  our  Investigation  of  possible  high- 
frequency  solutions  let  us  consider  the  only  conditions 
that  we  have  not  yet  studied,  namely  j  a  j  <  1,  Im  (i/a)  =  0. 

From  Appendix  IV,  we  have 


(kjooT  +  k2)  +  i(k-j_  -  k203T) 
(l^2  +  k22)xT 


(IV.1) 


Then,  the  condition  Im  (i/a)  =  0  leads  to 
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(10.82) 
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The  procedure  is  identical  with  that  in  the 
previous  section  but  now  o  =  1  in  the  exponential 
term  rather  than  a  =  2.  For  a  =  1,  we  easily 
find  that 


*1 
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1/2  1/a2 

ir  '  e  ' 
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(10.83) 


y 
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3  •  4 

-  75  icdt  a 


(10.84) 


y^'  =  -  icor 


(10.85) 


Let  us  consider  each  solution  y*  separately. 


(i) 


a  7 r1/2  e1/3" 
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(10.83) 


and  using  (10.82) 
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CO 
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2tt  7  x 


(10.87) 
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One  must  now  investigate  the  nature  of  k2  for  the 
region  cot  »  1.  For  odt  »  1  the  exponential  term 
approaches  unity  if 
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U  =  to 

2  27T1/2  x 

so  this  must  be  the  solution  of  (10.87).  Then 

i  i  cjd2t 

kl  =  k2“T  = 


so  that 


V 


2tt1//2x 

cot 


But  for  this  solution. 


|  a  |  =  k2xi  =  -"i/g  >:>  1 

2ir  ' 


for  cdt  >:>  1  contrary  to  our  original  assumption. 
Thus,  this  solution  for  k2  and  k-^  must  be  rejected 
because  it  is  not  consistent  with  one  of  our  original 
assumptions . 


. 
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(ii)  From  (10.84),  we  have 

V^T5  exP  <k2sW>  “  3E7-  (10.88) 

For  cot  »  1,  this  has  no  solution  since,  in  general, 

1/d*2  £T 

e  7  increases  more  quickly  than  g^  decreases  as 

g  0. 

(iii)  From  (10.85), 

k2XT  exp  (■--~2-g)  =  7T1/2  (10.89) 

kg  X  T 

Once  again,  for  cot  »  1,  this  equation  has  no  solution 
for  the  same  reason  as  in  (ii). 

Thus,  we  see  that  the  pair  of  conditions 
|  a  |  <1  and  Im  (i/a)  =  0  are  not  compatible  for 

consistent  high-frequency  solutions.  Without  going 
through  this  brief,  and  somewhat  crude,  analysis,  we 
could  have  predicted  this  since  Im  i/a  =  0  leads  to 

(10.80) ,  viz . 


k2/k 


1 


6 


_1 

COT 


■  . 


.  ■■ 


.... 


■  r  '  '  ,  •’  ' 

:  :  ■  ''  -<  t .  .y y:;.y  .  .  ■  -i j 


:  . 

i  ii  L  ■ hi  ;  i  '  L©  IC  '  '  "  U'l. 

h  J?:J 


l4o. 


In  terms  of  our  previous  notation.  Except  for  factors, 
(10.83)  -  (IO.85)  are  the  same  equations  as  were 
obtained  for  Im  (i/a)  <  0.  There  it  was  shown  that 
6  =  l/ cot  was  not  a  consistent  solution  for  cdt  »  1. 

10 . 5  Summary  of  High  Frequency  Results 

In  studying  the  high-frequency  sound  propagation, 
we  have  considered  all  possible  combinations  of  initial 
assumptions  because  we  lacked  experimental  data  to 
suggest  the  correct  choice  to  make.  The  possible  combin¬ 
ations  were:  (i)  |  a  j  >  1,  Im  (i/a)  arbitrary; 

(ii)  I  a  |  <  1,  Im  (i/a)  >  0;  (ill)  i  a  |  <  1,  Im  (i/a)  <  0; 
(iv)  |  a  |  <  1,  Im  (i/a)  =0.  It  was  found  that  only 
the  third  choice  (|  a  |  <  1,  Im  (i/a)  <  0)  yielded 

high-frequency  solutions  consistent  with  the  initial 
assumptions  and  also  physically  meaningful. 

A  single,  physically  meaningful,  result  was 
found  for  lot  »  1,  namely 

voA°  =  °-'42  (10.55) 

avo/co  = 


0.42 


(10.56) 
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A  rather  weak  argument  was  presented  for 
rejecting  higher-order  solutions  which  differed  from 
(10.55)  and  (IO.56)  only  in  the  numerical  factor x . 

Thus  for  cdt  -*•  00  ( ie .  for  very  large  frequencies) 

the  phase  velocity  V  approaches  a  constant  value 
which  is  approximately  2VQ.  The  amplitude  absorption 
coefficient  a  is  found  to  increase  linearly  with  the 
frequency  and  does  not  depend  on  the  relaxation  time  t. 


*  It  has  been  recently  pointed  out  by  Dr.  Schlff  that 
the  higher-order  solutions  resulting  from  the  multi¬ 
valuedness  of  the  inverse  tangent  may,  indeed,  have  physical 
significance  and  should  be  investigated  more  closely.  If 
the  multi-valuedness  is  real  perhaps  it  extends  into  the 
intermediate  frequency  region  and  could  be  checked 
experimentally . 

Since  the  writing  of  the  thesis,  this  author  has  come 
across  two  papers  in  plasma  physics  in  which  this  multi¬ 
valued  behaviour  is  present.  The  reader  is  referred  to 
the  following: 

R.  M.  Lewis  and  J.  B.  Keller,  Phy .  Fluids  j5,  1248,  (1962) 
Reference  50,  p.  211-213. 

Unfortunately,  neither  of  these  papers  attempts  to  discuss 
the  physical  implications  of  the  multi-valued  results . 
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Chapter  11.  DISCUSSION  OF  THE  RESULTS  FOR 
t  INDEPENDENT  OF  v 

Let  us  now  sum  up  what  we  have  obtained  for 
t  independent  of  v  .  For  cdt  «  1,  our  results  displayed 
the  correct  qualitative  behaviour.  The  results  agreed 
qualitatively  with  experiment  and  with  the  work  of  Wang- 
Chang  and  Uhlenbeck. 

For  cot  »  1,  we  have  found  ( equations  ( 10 . 55 ) 
and  (10.56))  that  the  phase  velocity  becomes  independent 
of  frequency  while  the  amplitude  absorption  coefficient 
a  increases  linearly  with  the  frequency.  Furthermore, 
a  does  not  depend  on  the  relaxation  time  t.  The 
high-frequency  results  were  obtained  for  cdt  »  4.4. 

It  is  our  belief  that  the  limiting  length  in  sound 
propagation  is  not  the  mean  free  path,  as  put  forward 
by  Primakoff  and  others,  but  the  average  interatomic 
distance  between  gas  particles.  This  means  that  we  cannot 
really  expect  our  results  for  cjot  £  100  (since  the 
average  interatomic  distance  is  of  the  order  of  one- 
hundredth  of  the  mean  free  path)  to  correspond  to 
physical  reality.  For  Helium  at  S.T.P.,  t  =  1.91  x  10-10 

sec.  so  that  cdt  =  100  corresponds  to  a  frequency  of 
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about  8  x  cycles  per  second. 

In  Fig.  7>  we  have  Indicated  the  behavior 
of  both  the  dispersion  and  absorption  curves  for 
cdt  »  1  and  got  «  1  and  have  further  compared  our 
results  with  experiment  and  with  the  results  of 
Wang- Chang  and  Uhlenbeck. 
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Chapter  12.  LOW  FREQUENCY  RESULTS: 

i  INDEPENDENT  OF  v 

Another  approximation  to  the  low-frequency 
problem  is  to  treat  $,  as  a  constant.  As  is  shown  in 
Appendix  VI,  if  one  takes  t  =  i/v,  where  l  is 
independent  of  v,  one  arrives  at  the  classical 

1  C  T/* 

ratio  of  xAl  *  namely,  xAl  =  — 5  —  •  This 
suggests  that  if  we  assume  t  =  f/v  we  might  obtain 
the  "classical"  absorption  (to  the  first  order  in  l ) . 
For  t  of  the  form  suggested, the  distribution  function 
is 
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(12.1) 

Expanding  the  denominator  for  low  frequencies, 
we  obtain  to  the  first  order  in  i 
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+  2PQvzu0  +  2P0lcn£  cosO  uQ 
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Proceeding  in  the  same  manner  as  before*  we 
can  obtain  the  number  conservation  equation 
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With  the  aid  of  equation  (12.3)*  we  eliminate 
Ni/No  from  (12.2)  immediately*  so  that 
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Then,  using  (12.4),  the  momentum  and  energy 
conserving  equations  are 
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For  this  pair  of  equations  to  have  a  non¬ 
trivial  solution. 
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which  yields 
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On  solving  for  the  real  and  imaginary  parts 
of  this  equation  and  letting  k  =  k^  +  ikg  (kg  «  k^), 
we  find  that 
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Prom  Appendix  VI 
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so  that 
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and  a  =  0.625  ^7-  (12.10) 

PVQ 

Thus,  we  see  that  by  taking  t  =  &/v  we  have  obtained 

some  slight  Improvement  in  the  slope  (from  0.600  to 
0.625)  of  our  absorption  curve.  However,  the  result 
still  does  not  agree  with  the  slope  of  0.700  that 
one  obtains  "classically".  Perhaps,  this  is  due  to  the 
expansion  that  we  used.  Certainly,  terms  with  higher 
powers  of  i  in  our  expansion  lead  to  divergent  integrals 
while  the  unexpanded  distribution  function  is  still  finite. 
This  suggests  that,  perhaps,  the  correct  expansion 
procedure  might  yield  the  classical  slope  for  the 
absorption  curve.  Yet,  why  did  we  obtain  the  correct 
ratio  of  xAl  (which  involves  only  the  first  order  in  $)? 
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Chapter  13 .  DIATOMIC  AND  POLYATOMIC  GASES 

Before  concluding  this  thesis,  we  shall  very 
briefly  mention  how  this  "relaxation  model"  may  be 
applied  to  the  problem  of  diatomic  and  polyatomic  gases. 
No  explicit  calculations  will  be  carried  out  here  but 
we  shall  merely  indicate  the  method  of  attack. 

The  qualitative  behaviour  of  sound  propagation 
through  polyatomic  gases  has  already  been  discussed  in 
Chapter  4.  Experimentalists,  such  as  A.  Smuda^2  (1951), 
have  pointed  out  that  the  existing  theory  does  not  quite 
fit  the  experimental  points  in  gases  such  as  Og .  Perhaps 
our  "relaxation  model",  if  carried  out  in  detail,  would 
rectify  this  shortcoming.  One  can  associate  certain 
relaxation  times  with  the  transfer  of  energy  from  the 
translational  to  the  rotational  or  to  the  vibrational 
degrees  of  freedom.  The  rotational  relaxation  time  xr 
in  most  gases  is  of  the  same  order  of  magnitude  as  the 
translational  relaxation  time.  However,  the  vibrational 
relaxation  time  tv  is  considerably  longer.  Associated 
with  the  Tr  and  tv  will  be  internal  energies  Ep 
and  Ey  or  characteristic  temperatures  Tr  and  Tv. 

In  any  case,  we  shall  have  two  new  "macroscopic"  unknowns. 


V  ■  ..  L  .  :  ■■■•■■)  :  '  .  /■  ..  ■ "  I :  l  .dr-  t  i  L 


Thus,  one  shall  need  two  mare  equations  to  supplement 
the  number,  momentum  and  energy  conservation  equations. 


These  two  additional  equations  may  be  taken 
to  be  of  the  form 


dTv,r 

dt 


1 


T. 


v,r 


The  new  distribution  function  f  (N,  u,  T,  ,  T  ,  T  ) 
will  still  satisfy  the  Boltzmann  equation.  The  momentum 
and  energy-conserving  equations  will  be  different  due 
to  the  additional  degrees  of  freedom.  The  five  equations 
in  the  five  unknowns  may  be  solved  to  give  a  general 
dispersion  relation  in  terms  of  tp ,  xv,  .  The 
details  of  this  calculation  are  left  to  some  future  date. 
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Chapter  l4 .  CONCLUSIONS 

In  concluding.  It  should  be  pointed  out  that 
our  methods  led  to  very  interesting  results  for  the 
behaviour  of  sound  as  cdt  ->  oo.  It  was  found  that  the 
phase  velocity  becomes  independent  of  frequency  once 
more  while  the  amplitude  absorption  coefficient  increases 
linearly  with  the  frequency.  The  latter  result  can 
qualitatively  be  explained  on  the  basis  of  the  following 
argument  due  to  Primakoff2  (1942).  As  the  sound  wave¬ 
length  becomes  sufficiently  short,  "a  molecule  can  travel 
without  any  intervening  collisions  from  a  region  where 
the  phase  of  the  sound  wave  (and  so  the  phase  of  the 
molecule's  ordered  velocity),  has  a  particular  value 
to  a  second  region  with  an  altogether  different  value  for 
the  phase  of  the  sound  wave.  Thus,  there  will  occur 
numerous  collisions  between  molecules  having  quite 
different  values  for  their  ordered  velocities;  as  a 
result  of  these  collisions,  the  ordered  or  "mass" 
motion  of  the  molecules  will  be  transformed  into  random 
thermal  motion  and  the  propagation  of  the  sound  wave 
violently  affected."  However,  unlike  Primakoff,  it  is 
thought  by  this  author  that  sound  will  propagate  for 
wavelengths  shorter  than  the  mean  free  path  but  not 
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shorter  than  the  average  Interatomic  distance.  It  is 
felt  that  for  wavelengths  shorter  than  the  average 
interatomic  distance,  the  concept  of  a  sound  wave,  or 
some  equivalent  collective  mode,  is  meaningless. 

For  cot  «  1,  the  calculations  were  somewhat 
improved  by  taking  i  to  be  independent  of  v  rather 
than  t.  However,  our  absorption  curve  still  did  not 
have  the  "classical"  slope. 

Several  suggestions  were  made  in  the  thesis 
as  to  future  experimental  and  theoretical  work. 
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APPENDIX  I 


Evaluation  of  integrals  of  the  form 
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We  are  interested  here  in  relating  1^, 

I^  to  I0 .  First  let  us  briefly  investigate  IQ .  We 
shall  discuss  the  properties  and  various  forms  of  I0 
more  thoroughly  in  Appendix  III. 
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so  that. 


47T 


T(l  -  Icot)  ° 


/”  v  tan-1  (yv)e-^°''/ ~  dv 


or  I 


o 


47T  roo  ,  -1/  V  V2  , 

t —  f  v  tan  (tv)  e  °  dv 
kx  J  o  v  ' 


(1.5) 


For  |  Tv  |  <  1,  (1.5)  may  be  expanded  by  using 

the  well-known  result 
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As  a  check  on  this,  we  have  directly  from  the 
definition  of  I 
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Q-^ovz2. 


7T 


/ 


dVr 


-oo  1  +  iyv. 


P0(l  -  ioyr) 
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Assuming  that  |  yv  |  is  small  on  the  average 


I 


o 


(f)3/2 

Ho 


(  1  -  l03T ) 


1 


which  i s  the  same  as  previously. 

A  further  justification  of  this  expansion  is 

that  the  integrand  is  largest  around  v_  =  0  where 

£ 

|  yv  |  is  small.  For  v  large  the  contribution  from 
the  integrand  is  negligible.  It  will  be  shown  later 
(in  Appendix  III)  that  the  asymptotic  expansion  naively 
derived  above  holds  only  in  a  certain  region  of  the 
complex  plane  and  one  must  take  the  correct  analytic 
continuation  in  other  regions. 


Now,  consider  equation  (1.5) 


I 


o 


(1.5) 


77 

From  the  Integraltaf el  of  W.  Grobner  and  N.  Hofreiter 
(1961),  we  find  that 


dx 


7T 


2 


e 


A 


[1  -  $■(/!)] 


(1.7) 
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where 


4>(  JT)  =  —  SJK  e'X"' 

/if  o 


for  the 


Re  /T"  >  0. 


Integrating  (1.7)  by  parts  yields 


-Ax' 


•o 


1  +  x' 


.oo  _  7\-v-  _  1 

-dx  =  2A  /  x  e  tan  x  dx 

?  Jo 


so  that  r  x  e~^x  tan-^  x  dx  =  ~  e^  [l 

4a 


Using  this,,  (1.5)  "becomes 


I 

o 


T  eP°/'y2  [1  -  ®(PoVy)] 


PokT 


where 


®002A) 


_2  Ao2/^  --X2 

A  o 


e  dx 


2jVl  fl  e-(P0A2)y2 
t/7  o 


(1.9)  may  also  be  written  in  the  form 


TT 


2 


o 


(1  -  1o>t)P  y 


e^o/Y  erfc  (f3Q2/Y) 


(1.8) 


4>(/A)] 


(1-9) 


(I. 10) 


(1. 11) 


It  is  easily  shown,  as  a  check,  that  for 
1 

Re  (P02/t)  >  0  (I. 11)  yields  the  same  expansion  as 

(1.6).  We  shall  now  relate  I,,  In  and  L  to  I  . 
v  '  12  3o 


(2) 


Iff 


+c 


vz  e 


-Pov‘ 


1  -  Icot  +  ikTVr 


dv 


(1.12) 


v3  e~^°v  cos 0  slnQ  d^dv 

1  -  icDT  +  ikxv  cos& 


.~f+l  v3  u  e_f3°v2 

-2 *  Lf 


„co 
'o  _i 


1  -  Icot  -  ikrvu 


dudv 


3 

i  i  o  ‘  O 


-2tt  j 00  j+1  v-3  u  e~Ho 

(l  -  Icjdt)  0  -1  l  -  iyvu 


dudv 


Let  y  =  l  -  iyvu  or  u  = 


Then,  dy  =  -  iYvdu  and  for  u  =  -1,  y  =  1+iyv 

u  =  +1,  y  =  1-iyv 

With  this  transformation 


/ 


+1 

-1 


u  _  1  fl-iyv 

du  0  p  J 

1  -  iyvu  y  v  1+iyv 


-  1)  dy 


-1  1-iyv  i_i<w 

-p-p  [log  y  -  y]  =  -p-p  [  log (  -)  +  2iyv 

y  l+lyv  ydv  l+iyv 
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2tt 


(1  -  Icot)  Y" 


o 


v  e 


-Pov 


C-21YV  -  log(-^^)Jdv 


but  l°s(~irrYv)  =  -  2  1  tan-1(  yv) 


2ir 


(l  -  1<dt)y2  ° 


/  v  e  ^°v  [-21YV  +  21  tan'^Y^]  bv 


-  4tt1 


2  -P„v‘ 


4tt1  r00  __  _  -f3„v2 


—  /  v  e-K°  dv  +  — 

(1  -  iaox)  Y  °  ( 1  ~  1o>t)y2  ° 


/  v  e_po 


tan-2(Yv)  dv 


or 


xi  “  71  <^)3/2  +  1  i0 

1  kx  P0  Y  ° 


(1.13) 


(3)  I0  =  /// 


+  00 


2  -PqV2 

vz  e 


”°°  1  -  Icot  +  ikTVr 


dv  (1.14) 


=  27 T  f  f 

Jr\  Jr 


00  17  v^  e~^°v  cos2d  sind  dddv 


o  •'o 


1  -  icsT  +  ikv  cos Q t 


27 T  r°°  r+l 


4  ,2  -Pnv2 


/  / 


v  u  e 


(l-iayr)  o  -1  1  -  ±YVU 


dudv 
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Now,  f 


+1  2  l-iyv  . 

- -  du  =  /  — - (^  -  2  +  y)  —1.  -  ■ 

1+iYv  (i7v)d  y  (-iyv) 


-11-  l^vu 


— T~^  [log!”"-)  +  21Tv] 


- -  l,  v  e"P°v2  [21TV  +  logN—^y-)  ]  dv 

iY+l-iorr)  0  1+1YV 


_ 47T 

72(1-1cjdt) 


47T 


-  J 

y^(l-io)T)  ° 


v  tan"^(Yv)  e”^°v 


dv 


1 


Ykx 


(+)3/2  -  4 


Y 


2  Jo 


W  I3  =  /// 


+  00 


3  -Pov 
v  e  ° 
z 


dv 


-00  1  -  IcjDT  +  IkTV 


-2tt  f°°f+1  e~^°v  u^v^ 
(l-iorr)  Jo  j-l  1  -  iyvu 


(1.15) 


(1.16) 


dudv 
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J 


+1  _u£du_  =  r 1_lyv  _ rl  o  _|_  Oy  v2i  dy 

_1  1-iyvu  i+iTV  (iYv)3  y  3  3y  y  J  ( -iyv ) 


1  _ /  1-lyv  ^  21  ,2  1 

yv 


7*7  1°s(i+iyv)  '  y3*3  +  3 


21  21  ,  2  1 

4  4  tan  "  3  3  +  3  yv 

y  v  y-'V  ' 


I 


47Tl 


(l-ioDT)y^  ° 


i  R  2 

/  v  tan”  (yv)  e”  °v  dv 


,  47ri  r°°  2  -Pnv2  ^ 

H - - -  I  v  e  °  dv 

J  n 


y2 3(  1-iayr) 


47ri  r  00  4  -f3  v^  ^ 

J  v  e  0  dv 


3y(  1-1odt) 


3 


1  (^)3/2  -  -i-(^)3/2  -  i 


2  {W~ 

ikxy  po 


2krP0  po 


Y 


o  I 

3  o 


(1.17) 
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APPENDIX  II 


^2  T 

Relating  — £  ,  - 2  to  I  . 

sp0  ap02 


!0  =  Iff 


+  00 


.-PoV‘ 


1  -  io)T  +  ikTVr 


dv 


(I.D 


7 T  +0°  e-PoVZ 

ft  **  1  -  i(JDT  +  ikxv  dvz 

Pq  -00  Z 


Let  60v72  =  y2  or  B  2v 


o  z 


=  y  and  also  put  x  =  1/(3q2 


Then 


VH572  / 


+  CO 


-y 


6  O/ c-  -oo  1  _  1Q3 T  +  ikXTy 
o 


dy  = 


77x3 


+oo  -y2 


I 


ikxT  -oo  y  + 

1KX  T 


dy 


or 


2  +oo  -y^ 

I  =  —  /  - - —  dy 

0  ikx  _0°  y  -  t, 


(II. 1) 


where  £  = 


/1-iCQTN 
V IkXT  ' 


(  1+CDT  ) 

kXT 


(II. 2) 


Strictly  speakings  this  representation  is  defined 
for  the  Im  £  >  0  while  for  Im  £  <  0,  one  must  take 
the  analytic  continuation.  The  properties  of  this 

74  75 

representation,  variously  known  as  the  Hilbert  transform  *  " 
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69  70  71 

of  the  Gaussian  or  the  plasma  dispersion  function  *  ’ 

will  be  discussed  In  Appendix  III. 


We  may  write 


— - — -s  =  1  /  dt  valid  for  Im  £  >  0 

y  -  C  o 


so  that  I 


7TX 


+00.00  -lty+lt^-y2 


r  -  I  f  e 

kT  J0 


dtdy 


Carrying  out  the  y  Integration  and  letting  y  =  —  s 

v 


/+“  e-^y-T2  dy 


-00 


i  ,  its 

X2  r+“  e-  — 


1  2 


-c 


ds 


1  (-1)2 
ri  e"  2  ^ 


i  -t2/^ 
e  / 


Here  we  have  just  taken  the  inverse  Fourier  transform. 


I 

o 


++£ f  -  t£>  dt 

k  T  o 


(II. 3) 


It  can  be  shown  that  this  representation  is  valid  for 
all  £.  We  may  write 

L(  £ )  =  /°°  elt^  t  /4  dt  (II. 4) 

o 


■  ; ' 
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so  that. 


3/2  2 

TT'J/  X 

lCT 


L(0  = 


7 r 


3/2 


L(C)  (II. 5) 


We  wish  to  determine  dl  /df3  an^  /d(3 


o'  "o 


“o/  ^  o 


dl 


o 


-7T3/2 


o  P02kT 


L(C)  + 


7T 


3/2 


c)L 

p0kT  ap0  ac 


„3/2  „.3/2,  Sr 

^-js—  *■(?)  +  p^-  — 

PQ  kt  2^o^kT  ^ 


Now,  (  L(0 


1  /”  dt  e-4  /4  •  — (elfc^) 
1  Jo  dt 


and  then  integrating  by  parts 


C  L(c)  -  i  e-t2/4  eitc  I  “  .  1  ro  (=f)  e-t2/^  +  “C 


dt 


=  -  I  1  /”  it  >  +  “C 

i  2  Jo 


dt 


1 

2  ^ 


i  -  4  L«(0 


so  that 


=  2  [i  - 


L'(C) 


?L(C)] 


(II. 6) 


3  '  ■  ■  • ' 
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Si 


o 


o 


•!■(?)  +  ^^7  [i  -  £  L(C)] 
f>o  1"  p0  kT 


and  since  C  =  t— 


i 

Yx 


(II. 7) 


Sp 


1_ 

P, 


0  Ykx  ^o; 


±-(T-^/2  +  1  j 

y2  O 


Yo,3/2  5i 
or  (— >  3F 


o 


o 


(_  1_  +  1  )(£°)3/2  T  1 

Po  7T  '  ' 


■o  YkT 


(II. 8) 


Similarly, 


S2I 


Sf3 


o 


P. 


2  o 


I  - 


1  Si 


£  + 


PQ  Sp0  2Y1ctP0 


3  ( 7T  \  3/2  ,1  Olo 

—  +  7  sp; 


1  I.  +  -Yr  I_  +  - - -  t?)3^2 - Y-  I 


Po2  °  ^o2  °  7kTS0 


y2s  0 

o 


+ 


3  ,7^  3/2  1 


2YkrP0 


1  (2L_)3/2 


Y2S0  °  Y^t 


+  -=T7  I 

Y+  o 


7+  1  Jo  +  WfiZW 


5 _ /  7r  \  3/2  _  1  /7T  \3/2 


Y  ■  ~  ^  i^o  '  ^  7^kx  ^ 


[_£ _ 2 _  +  li  (£o)3/2  I 

Po2  y2p0  7  t 


+ 


2yk;TPo  Y^kr 


(II. 9) 
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APPENDIX  III. 


write 


As  we  formally  proved  in  Appendix  II,  we  may 


"o'0  -  IT.  C  fh  ^  (IIa) 

where  £  =  l-3:  -t.  (II. 2) 

kxx 

B.  D.  Fried  and  S.  D.  Conte^  have  defined  the  so-called 
plasma  dispersion  function  as 

i  +oo  p-y2 

Z(C)  =  ^770  /  — -  dy  (III.1) 

-00  y  -  C 

for  Im  £ >  0  and  as  the  analytic  continuation  of  this 
for  Im  0.  Therefore,  we  have 

7 r3/2  x2 

Io(0  =  ~ - —  Z  (?)  (III. 2) 

ikx 


for 

Im 

£>  0 

and  as 

the 

analytic 

continuation  of  this 

for 

Im 

o. 

Since 

*0 

and  Z 

are  simply  connected 

we  shall  study  the  properties  of  Z.  It  should  be 


- 
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pointed  out  that  Z,  as  given  by  (Ill.l),  occurs  in  many 
plasma  physics  problems.  The  interested  reader  is  referred 
to  such  papers  as  those  by  S.  Rand^  (i960),  K.  M.  Watson 
et  al^  (i960)  and  E.  A.  Jackson^  (i960).  The  problems 
investigated  range  from  "Satellite  Wake  in  the  Ionosphere" 
to  "Drift  Instabilities  in  a  Maxwellian  Plasma" . 

70 

The  first  thing  to  be  noted1 ^  about  Z  is 
that  Z(£*)  is  not  the  analytic  continuation  of  Z (£). 

By  expanding  (y  -  as  a  Fourier  integral 

viz.  — 1—  =  i  /“  e-lu(y-^  du 

y  -  C  ° 

and  then  carrying  out  the  y  integration  as  in  Appendix 
II,  we  have 

Z(0  =  i  /°°  elu£"(u/2)2  du  (III. 3) 

o 

Now,  letting  y*  =  iu  +  £  so  that 

dy'  =  -|  idu  and  for  u  =  o,  y '  =  £ 

u  =  00,  y1  =  ioo 


' 
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(III. 3)  becomes. 


z(C)  = 

2  e2C(y-0+(y'-C)2  dy, 

or  Z(0  = 

ioo  y  •  2_  p 2 

=  2  /  ey  dy* 

_p2  ioo  vi2 

=  2 e  f  ey  dy  * 

C 

and  letting  y' 

1  =  is 

z(C) 

_j»2  co  _  g  2 

=  2ie  ^  /  e  ds 

-1C 

or  Z(C) 

»2  1C  c2 

=  216"^  /  e  ds  (III. 4) 

—  00 

This 

alternate  representation  is  valid  for 

both  signs  of 

Im  £.  These  last  few  representations  have 

been  used  by  L. 

,  Landau^  (1946)  and  by  A.  Vlasov^0  in 

their  study  of 

plasmas.  The  analytic  continuation  of 

(ill.l)  for  Im  £  ^  0  may  be  found  by  using  (III. 4) 

If  we  let  £  =  ift  =  !(■-"---) 

krx 


so  that  Q 

=  (1“layr)  =  1  _  Po2  (III. 5) 

kxT  yx  y 

where 

Y  -  d-Sr)  (I'3) 

then. 

n2  oo  2 

Z  =  2ie^  /  e-s  ds  (ill. 6) 

Q 

and  for  the  real  part  of  Q,  >  0,  this  may  be  written 


as 

Z  =  i  fw  erfc  0Q2/y) 

so  that 

I  =  - — -  e^o/V  erfc((3  2/y)  (I.ll) 

°  (l-iO)T)P0Y 

which  is 

just  what  we  found  in  Appendix  I.  This  confirms 

the  fact 

that  the  series  expansion  obtained  there  holds 

only  for  Re  Q  >0  and  that  for  Re  £7  <  0  we  must  take 
the  analytic  continuation  in  the  manner  prescribed  above. 

Now  that  we  have  discovered  several  representa¬ 
tions  for  Z  ,  let  us  investigate  its  properties. 

(l)  Asymptotic  Expansion:  Z(£)  exhibits  a  Stokes 


phenomenon.  Ie .  different  asymptotic  expansions  are 
required  in  different  portions  of  the  complex  plane . 


. 
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The  Stokes  phenomenon  is  a  characteristic  of  the  asymptotic 
expansions  of  analytic  functions. 


Lemma:  For  |  £  |  -*■  o°  ,  Z(£)  has  the  asymptotic  expansion 


Z(C)  -  l7ri/2  0  e-c2  .  1U  +  _1_  +  +  ...  ] 


=  IT1/2  a  -  z 


n=o  it 


l/2^2n+l 


(III. 7) 


where  o 


0  Im  C  >  0 

1  Im  C  =  0 

2  Im  £  <  0 


Proof:  Various  proofs  of  this  lemma  have  been  given  by 

O  i  O  Q 

M.  Moshinsky  (1951)*  V.  N.  Faddeyeva  and  N.  M.  Terentev0^ 
(195^-)*  F.  Berz^  (1956),  G.  Beck  and  H.  Nussenzveig^ 
(i960),  B.  D.  Fried  et  al ^  (1961)  and  finally,  T.  H. 

Stix^  (1962).  Our  proof  closely  parallels  the  latter. 


Consider  the  representaion 


00 


I 


Cl 


(III. 8) 
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We  first  assume  that  Re  Q  >  0  which,  since  £  =  if2, 
corresponds  to  Im  £  >  0 .  To  find  the  asymptotic  expansion 
of  Z  ,  we  choose  as  path  of  integration  the  straight  line 
y  =  Im  Q=  constant  >  0  parallel  to  the  real  axis. 


le . 


Z  =  2ieu  / 


00+ iy 

fi=x+iy 


-s2 

e  ds 


Now  letting  s^  =  t  ds  =  t”^^  dt 


00  p 

—  R 


/  e“  ds 

a 


2  J^2 


h?  y 


-i/2  e-y 


e  J  dy 


and  integrating  by  parts 


i  r  e-y 
~  2  L~  ^172 


y  7  fi2 


+  (-^ 


1  [  e _ 

2  L  a 


+  (~  |)(-D  ^2  f 

y 7  ^2 


1  /  3  \  r00  e  ~  * 

2  2^  /  372  dy 

fr  * 


-y 


1  -n2 

—  e 

2Q 


1  - 


2QC 


+  |  n  e^C 


f 


-y 


nd  y 


572 


dy  ] 


+  (-  V 

2 


( 2n-l) ! ! 


1 

fi2n 


+  Rn(ft)  ]  (III. 9) 


.  • 
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where  Rn(ft) 


(-  |)n  (2n+l) ! ! 


„  fi2  f”  e‘s2  ^ 

n  s2n+2  d 


(III. 10) 


Now,  to  prove  the  validity  of  the  asymptotic  expansion, 
let  us  investigate  the  remainder  term  Rn(ft) . 


Now 


efi2  / 


-  s 1 


ft 


,2n+2 


ds  = 


ft* 


oo+  iy 

r 

x+iy 


-  s 


,2n+2 


ds 


Let  s  =  t 1  +  iy  so  that  we  have 


oo  ft2-(t'+iy)2 

/  - - n — UTo -  dt 1 

Jx  2n+2 


oo  x2- 1 1 2+2iy(x-t 1  ) 

/  - - FT-r^ -  dt* 

J  x 


,2n+2 


Now  t*  >  x  and  x  >  0  ,  so  that 


ft2  r00  e  s  d  <  r°°  dt 

"  h  ,2n+2  ds  L  2n+ 2 


(2n+l) 


Then 


Rn|  <  2-n  (2n-l)>!  \-i-  | 


n 


2n 


( III. 11) 


and  as  |  £  |  =  |  ft  |  -►  oo,  |  Rn  |  ->  0  thus  proving  the 

validity  of  the  asymptotic  expansion.  We  have  proved 
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Z  -  -  [  1 

Cl 


1  I  3 
2  Cl2  4^? 


..  ] 


Re  Cl  >  0  ,  or  since  C  = 


z  -  -  i  [  1  +  -ig  +  A  +  ...  ] 
c  2tr  4c 


( III. 12) 


for  Im  C  >  0  . 


Now  suppose  that  Im  C  <  0  or  Re  ft  <  0 


q2  oo  q  2 

Z  {Cl)  =  2ieu  /  e~s  ds 

Cl 


where  Re  Q  <  0  =  -  x  ,  say,  and  x  >  0 . 

Let  us  integrate  again  along  a  line  parallel 
to  the  real  axis.  Then  we  may  consider 


■y-2  oo  _«2 

Z( -x)  =  2iex  /  e~  ds 

-x 


y2  o  c;2  °°  q? 

=  2iex  [  /  e"  ds  +  /  e~  ds  ] 
-x  o 


0.  x2  r  yir  r°  -s2  ,  -j 
=  2ie  L  -^2“  +  J  e  ds  J 


-x 


l8l. 


=  2iex2  [  4T +  JX  e's2  ds  ] 

2  o 

=  2iex2  [  +  q  -  /“  e-s2  ds  ] 

-y  2  00  c.  2 

=  2ie  [  yW  -  fx  e~s  ds  ] 

v2  x2  00  _  g 2 

=  2i/i T  eA  -  2ieA  /  e  ds 


Z(-x)  =  2i/7r  ex2  -  -  [  1  -  — +  .  .  .  ] 

x  2x^  4x^ 


Z(y)  =  21YY  ey2  +  -|  [  1 - ^  +.  .  .  ]  y  >  o 

y  2 jd  4y^ 

or  changing  from  the  dummy  variable  y  to  Q, 


Z(Cl) 

or  Z(?) 


02  i  1 

2iv/7f  eu  +  -  [  1 - -  + 

«  2fi2 


.  ]  , 


2it7r  e 


-  ^2  1 


1 


_  ±  r  i  +  _ 

£  2£2 


+  ...  ] 


Re  ft  <  0 

(III. 13) 


for  Im  £  <  0  . 
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What  about  the  case  Im  £  =  0? 


Then  £  =  x,  a  real  number 


_*2  lx  _<»2 

Z(x)  =  21e  /  e  ds 


from  (hi. 4) 


=  2ie  f 


2  o+ix  2 

-Xc  p  -s  , 

e  ds 


-oo+10 


(since  this  is  what  is  really  meant  here) 


-x2  °  -<*2  lx 

=  2ie  [  /  e  ds  +  /  e  ds  ] 
-  oo  io 


Let  s  =  iy  in  the  2nd  integral 


Z(x)  =  2ie'x2  [  M  +  i  /X  e^2  dy  ] 

d  o 

=  i/jf  e”x  -  2e”x  /  e^  dy  (III.14) 

and  we  are  interested  in  expanding  the  integral  for 
I  £  I  =  I  x  |  -  oo  . 

Now, 

_y2  x  v2  X  v2-x2 

e  /  e^  dy  =  /  ey  dy 

o  o 
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Let  y2 


Then 


x 

/ 

o 


2  2  2  2  2  /""2  2 
—  ~  _  -  —  i/x  -y 


-x  =  —  s  ,  or  x  -  y  =  s  ,  or 


1  1 


27^5  <-2y>  dy  -  ds 


2-»2)V2 


or 


-(X  -s') 


dy  =  ds 


v2_x2 

ey  x  dy  = 


■°  e's2s 


-  /  ^=4  ds  =  / 


x  se 


—  s 


^  '/x-S2  0  /x-s2 


ds 


1  rx  se 


—  s 


=  -  / 

x  o 


ds 


/I  - 


T2 


1  s2  ,  in-  il  (Z£)2 


1  f00  0„_s2  r  .  i  B“  , 

x  /o  [  1  +  2  x2  +  2!  -  X2 


(-  i)(-  |)(-  f)  s2  3 

+  — ^ — § — —  +  ...]  <js 

3!  x2 


2  .00  _  0  2 


-  /  se"s"  [  l+--^o  +  --Tr  +  —  -^7:  +  ...]  ds 
x  Jo  2x2  8?  16? 


I  [  1  +  _J_  +  3  +  _15  +  ..I 
x  2  4x2  8x^  l6x^ 
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1 

2x 


[  1  + 


_I_  +  3  +  15 

2x2  4? 


T5 


+ 


(III. 15) 


for  x  -*•  00  so  that 

Z(x)  =  iff  e'x2  -  i  [  1  +  ^  +  ...]  (III. 16) 

for  Im  £  =  0. 

Thus,  we  have  proved  our  lemma.  A  more  rigorous  proof 
is  outlined  in  B.  D.  Fried,  Gell-Mann  et  al^1.  As  a 
matter  of  interest,  J.  Rosser^  (1948)  has  calculated 
the  first  twenty  terms  in  the  asymptotic  expansion  of  Z, 
viz . 


Z 


ITT 


1/2 


o  e 


-C‘ 


_i_  +  _3 _  +  15 

2£2  (2£2)2  (2£2)3 


105  ,  945 

(2C2)4  (2C2)5 


10.395  +  135,135  +  2,027,025 

(2£2)6  (2?2)7  (2£2)8 


34,459,425  +  654,729.075  +  13,749,310,575 
(2^2)9  (2?2)1°  (2J;2)11 


316,234,143,255  +  7,905,853,580,625 
( 2£2) 12  (2£2)13 


'  ■■  ' 

•*  -)s 

>  ,  ,  ■  '■■’  ..  .  . 

•- 


■  .  ■  ;  ..  ,/  "  ■  ,.c;'  '  r  nj;  .  ,  u;  :  1  a.; 

'•o  -  :•  v  ;  ruij  nl  m<i&$  .^n&wcf  tonil  sciS 
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213,^58,046,676,875 

(2?2)14 


6,190,283,353,629,375 

(2C2)15 


+  191,898,783,962,510,625 

(2^2)16 


6,332,659,870,762,850,625 

(2?2)17 


221,643,095,476,699,771,875 

(iF)13 

8,200,794,532,637,891,559,375 

(2^2)19 


(III. 17) 


(2)  Power  Series: 


2(0  - 


•  1/2  -£‘ 
itt  7  e  s 


-  2C  Cl 


— 

3  ;  15  C 


8 


C6  + 


105 

(ill. 18) 


e’t2  -  !  I  rl/2/{n  +  U 

n=o  ^ 


(hi. 19) 


which  converges  for  all  £ . 


l8'6. 

Proof:  Consider  the  representation 

-t2  -s2 

Z  =  2ie  ^  /  e  ds 

-00 

Formally,  we  have 

_r2  o  _s2  _  /*2  if  «2 

Z  =  2ie  ^  /  e  b  ds  +  2ie  ^  f  e  8  ds 

-  oo  O 

and  letting  s  =  iy£ 

Z  =  It1/2  e<2  -  2?  e<2  /  e ^  dy 

o 

=  It1/2  e*?2  -  2?  e"^2  f1  [1+  £2y2  +  +  ...]  dy 

o  ^  • 

=  i7Tl/2  e‘?2  -  2£(l  -  ?2  +  |^)(1  +  ji2  +  ^  c4  +  .  .  .) 

=  iir1/2  e'?2  -  2?  [  1  -  |  C2  +  2;  C4  -  .  .  .  ] 

for  any  £,  thus  completing  the  proof. 

(3)  Z  satisfies  the  differential  equation 

Z*  =  -2(1  +  £ Z)  (III. 20) 


for  all  £. 


187. 


Z(0)  =  1 TT1//2  (III. 21) 

Numerical  values  of  Z(£)  are  readily  generated  with  a 
computer  by  Integrating  the  differential  equation  and 
using  Z(0) . 
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APPENDIX  IV 

Movement  of  the  pole  in  the  complex  £  plane. 

The  position  of  the  pole  of  the  plasma  dispersion 
function  Z (£)  in  the  complex  £  plane  is  given  by 


•  _  i  _  i(  I-Icdt) 

a  (k1+ik2)xT 

=  ( i+cpT)  (kj-ikg) 

(k-L2+k22)xT 

(kjCDT  +  k2)  +  i(k-|_  -  k2o)T) 

or  C  =  - o - o - - 

(k-j_  +  k2  )  xt 


(II. 2) 


(IV.1) 


For  a) t  «  1,  we  have  found  (to  the  first  order  in  t) 
that 


k-,  =  ~ 


a 


=  k. 


3 

5 


(7.12) 


where  V, 


(-5KTovl/2  ^  so  that  (IV.1)  becomes 
3m  ' 


On  p 

03  T  .  3  03^  T\  .  .  /03  3  03^  T 

■V  5  Vj  VV_  "  5 


5  V, 


C  = 


03 


+ 


9  03 


4 


V02  '  25  V  2 
o  o 


)  XT 


:  ■  "  .  ■ 


r 


< 
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Therefore, 

q  V  +  ±V°(1  -  £  co2t) 

50  “To  2I  (IV-2> 

(1  +  ^9  cd2t2)  x 

0  V 

as  cot  0  ,  £  — — h  ioo  =  1.45  +  ico 

(IV.  3) 

Similarly,, 

for  cot  »  1, 

k-.xT  =  k2xx  =  0.46  co  (10.55),  (10.56) 


so  that 

c  = 

(0.46  co2t2  +  0.46  cot)  +  1(0.46  cot  -  0.46  co2t2) 

2  (0.46)2  co2  t^ 

and  as  cot  00 

£  -  1.09  -  1.091  (IV. 4) 

These  results  are  indicated  in  Fig.  8.  The 
exact  path  traced  by  the  pole  as  the  frequency  goes 
through  intermediate  values  (cot  ~  l)  has  not  been 
determined  here.  It  has  been  symbolically  indicated  by 
the  dotted  line  in  Fig.  8. 


*  •  '  '  ■ . 
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Im  £ 


1 . 45  t  loo 
A 


0)T  “*■  0 


I  2 


Re  £ 


/ 


0)T  00 


1.09  -  1.091 


Pig.  8. 


Movement  of  £  =  i/a  with  cot  . 
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APPENDIX  V 


Calculation  of  r)  and  x  for  T  Independent  of  v. 

(i)  Static  viscosity  (ri)  ;  The  static  viscosity,  tj, 
is  measured  when  no  sound  wave  is  passing  through  the 
gas,  but  a  macroscopic  velocity  gradient  is  imposed  on 
the  fluid.  Ie .  no  macroscopic  temperature  gradients, 
etc.,  are  present. 


If,  initially,  no  velocity  gradient  is  present 


—  +  v  — 

dt  dz 


If  a  macroscopic  velocity  in  the  x  -  direction. 


A 


u 


u(z)i. 


is  present 


o 


and  assuming  that  PQu2  is  small  and 


N 


N0(l  +  A(z)) 
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Then, 


Let, 


Also, 


Then, 


f  =  fQ  [1  +  A(z)  +  2P0vxu(z)]  (V.l) 


5f  Sf 

"St  +  vz  Si 


(f-f) 

T 


f  =  f  +  fn 

O  1 


fQ  +  hu(z) 


(V.2) 


f  =  f  +  F, 
o  1 


the  Boltzmann  equation  becomes 


df 


v. 


1  


 fl  +  fl 


dz 


(no  explicit  t  dependence) 


or  f . 


df- 

fi  ■  'r'z  ¥ 


or  hu  = 


?1  -  '"z11  Is 


hu  [l  +  v  T 
2 


I  .  in  ] 

u  5z 


f 


1 


hu 


_ £i 

1  +  vzT 


i 

u 


du 

dz 


193. 


hu 


fQ  [A  +  2PQvxu(z) ] 


1  +  vzT 


T  5u 

u  ’  "5z 


using  (V.l).  Expanding  the  denominator,  we  have 


hu  =  fQ  [A  -  AVZT  •  i  -|^  +  2P0VXU  -  2Pqvxvzt  ] 

1 

The  second  term,  involving  A  •  —  *  ,  is 

of  the  second  order  in  smallness  so  it  may  be  dropped. 

(It  integrates  to  zero  anyways  in  the  following  calculations). 
Then, 


f  =  f  +  hu 
o 


f 


o 


[1  +  A  +  2f30vxu  -  2P0vxvzx 


ciu 

dz 


] 


(V.3) 


The  pressure  tensor  becomes 


m  Iff  vxvzf  dI  = 


2m^0T 


du 

•Si 


iff 


v  2v  2f 
X  z  o 


dv 


-  m  N0  t  c>u 
2f30 


(V.4) 


' 


..  ■,  L'i  ;  at 


19^. 


m  N  t 

n  =  — —  =  n0kt0t  =  px 

2  P0 

where  p  is  the  static  pressure. 

Thus,  we  finally  have 

n  =  PT  (V.5) 


Since  we  assume  that  t  is  a  quantity  character¬ 
istic  of  the  gas  under  study  rather  than  of  the  process 
to  which  the  gas  is  subjected,  it  must  be  the  same  as 
the  t  derived  in  connection  with  the  absorption  and 
dispersion  of  sound. 


(ii)  Thermal  conductivity  (x) :  To  measure  the  static 

thermal  conductivity,  x>  we  impose  a  macroscopic 

dT 

temperature  gradient  on  the  gas. 

oz 


Then, 


f  (N  ,T  ) 
o  o  o' 


m  2 

f  =  N(__m  )3/2  e"  2KT 

27rKT  ( z ) 


T  V  t 

'  o 
N  4  N 

o 


where 
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Again, 


f  = 


f  - 


l£i 

Sz 


Write  =  g(v)[T  -  T  ] 


so  that 


f  =  f  +  g(T  -  T  )  =  f 

o  &v  o' 


St 

.Tg  -r- 
oz 


Thus, 


g(T  -  TQ) 


f  -  f, 


v„  t  St 


(l  2 — I — ill  ) 

(T-T0)  3Z 


Now 


,  f  =  fo  [1  +  A+P0v2  (-Z-2)  -|(^2)] 


g(T-T  )  =  f  [A  +  P  v2  (^1^2)  _  i 

O'  O  O  '  rp  '  O  '  rn  '  ■ 

^O  iO 


[1  - 


V„  T 


z  *  St 


(T-Tn)  Sz 


] 


f  =  f  [1  +  A  +  P  v2  (^1^2)  -  i  (^zlo) 

O  O  X  rp  /  0  V  rp  > 

■Lo  2  1Q 


-  Avzt  £T 
(T-T_)  ^ 


^ov  VZT  ST  ;  3  vzT  ST  ] 

T  Sz  2  T  Sz 

o  o 


(V.6) 


.  . . .  '  l~  . .  •••  .  -  . . 


C 


Now, 


or 


where 


Then, 


which 
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m  ///  f  v2vz  dv  =  -  2x  H  (V.7) 


-mp0N0T  Jr  _  35  i  3  _  mNQT  _  &T  _ 

dz  4  (3q^ 


m 

o  "  2KTQ 


Sz  8p  3  2  t 

o  ( 


■  -*£ 


X  =  5  |  •  N  KT  t 

^  m  00 

or  X  =  5  |  •  PT  (V.8) 

X  =  5  -  (V.9) 

T|  m 

is  different  than  the  classical  result,  namely 


I  =  15  K 

r\  4  m 


(V.10) 


'  • 
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APPENDIX  VI 


X  and  r|  calculated  for  %  Independent  of  v. 


(i)  Calculation  of  y:  From  Appendix  V,  we  have 


T-T^  o  t-T, 


f  =  f  [1  +  A  +  P0v2  (1-1°)  -  |  (±1±0) 


o  '  rp  /  p  '  m 
^o  o 


2  i  St  ,  3  i  St 

PoV  zT  T0  Sz  '  2  z  T0  $>z 


+  smaller  terms  ] 


Now,  m///f  v  v  dv  =  -  2x 


ST 

Sz 


(VI. 1) 


defines  x  •  Letting  t  =  H /v  (VI. 1)  becomes, 

using  f. 


ft  «  No  St  r  r  r  3  2  -Pov' 

-mpj  —  -r-  JJJ  VJV_  e 


°  T0  Sz 


dv 


+  |mN^  &  2  e.poV2 


2  Tq  Sz 


dv  = 


or 


-4mN0i5  x3  2mN  A  3 


T. 


o  t r 


T72  +  T 


o 


"T72 

Tf  ' 


=  -  2X 
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or 


X 


2KNq^ 

PQV2  W 


(VI. 2) 


(11)  Calculation  of  ti  ;  From  Appendix  V,  we  have  for 
t  =  4/v 


vxv 

f  =  fQ  [1  +  A  +  2P0vxu  -  2P0  —  I  ^  ] 


Then  mfffv  v  fdv  =  - 

X  Z  —  tct 


8  mV 


1  8u 


15  772  '  JJ7?  5, 


=  -  r|  ~  defining  4 

o  z 


Therefore  *  r\  =  — 


8  mNo  1 

15 


(VI. 3) 


From  (VI, 3)  and  (VI. 2),  we  have 


^  ^  -  (VI.  4) 

rj  4  m 

which  is  just  the  "classical"  ratio  of  shear  viscosity 
to  thermal  conductivity.  The  "classical"  ratio  differs 
from  the  experimental  results  by  no  more  than  2%  for  all 
monoatomic  gases. 


a ",:v 


.  .  I 


■  5X3‘''.nV 
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APPENDIX  VII 


Solution  of  the  Cubic  Equation  (10.26). 


y^  [  l6a^  -  l6ici)Ta^  +  32io)Ta^  ] 

+  y2  [  8cd2t^  -  4a^  -  32oo^T^a^  -  6a^  +  12ia>Ta^ 


+  30o3^T^a^  -  12a^  +  48io)Ta^  ] 


+  y  [  +  4100^x3  -  2icjDxa^  -  2o)^x^a^  -  8io)^T^a^ 

+  a^  -  31(^>Ta^  +  tu^T^a^  +  17iao^T^a^ 

-  12icDTa^  -  54co2T2a^  ] 

-  3^T^a^  +  OD^T^a^  -  llicD^T^a^  =  0  (10.26) 


where 


A/2  e1/32 


y 


a 


(10.21) 


■  -  .  .  ->  ;  :  :  ■ 


"  e'  T 
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We  shall  solve  (10.26)  for  the  three  solutions 
yx,  ~$2>  ^3  ln  standard  way.  (10.26)  Is  of  the 

form 


Ay3  +  3By2  +  3Cy  +  D  =  0  (VII. l) 

We  can  put  It  into  the  form 


z3  +  3Hz  +  G  =  0  (VII. 2) 

where  z  =  Ay  +  B  (VII. 3) 

H  =  AC  -  B2  (VII. 4) 

G  =  A2D  -  3 ABC  +  2B3  (VII. 5) 

If  one  forms 

p  =  1  [_g  +  /g2  +  4h3  ]  (VII. 6) 

q  =  1  [ -G  -  Jg2  +  4H3  ]  (VII. 7) 


then  W  +  W  Is  a  solution  of  (VII. 2) . 

4 

Then,  from  (10.26),  we  have,  keeping  terms  to  a 


■ 
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=  [32icoTa2  +  l6a^  -  l6icoTa^]y 


+  [  §  cd2t2  -  |  a2  -  3|  a)2T2a2 


-  2a 


4 


+  4l<jDTa^  +  10(jD2T2a^  ] 


(VII. 8) 


H  =  [32ia)Ta2  +  l6a^  -  l6io)Ta^]  •  [-^  1o)3t3 

-  Icnra2  -  o)2T2a^  -  ia^T^a^  +  ~  a^ 

-  IcDTa^  +  o>2T2a^  +  io^a4] 

.  [  |  o,2t2  _  ^  a2  .  32  (D2T2a2  _  2  a4 
+  4  io)xa4  +  10  o)2T2a^ 


+  a)4T4a2  -  64  ico3T3a4  .  544  ^4^4 


+  o)2 -  6|  ^4  .  16  a4  +  64  ^^2 


(VII. 9) 


. 


y 
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Then, 


+  oo12T12a2 

729 


9,830,400  12  12  4 

— — =L—1 -  a)  t  a 

729 


2,752,512  ^lO^LO^ 
729 


I8M32  ia3llTlla4 
243 


262,144  12  12  983,040  8_8„4 

- - -  a)  t  -  - - -  a)  t  a 

729  729 


786^432  m10x10a2 
729 


(VII. 10) 


2  o  2  4 

A  D  contains  no  terms  in  a  ,  a  ,  a 


-3 ABC  =  (-AC)(3B)  =  +  oo6x6a2 


1024  444  6656  664 

— ~ —  cjo  t  a  -  — *r—  00  t  a 


(VII. 11) 


2B3  =  2(B2)(B)  = 


+  ^  ^  ,4T4a2 


-  o)6x6a2  +  ^4T4a4  +  1536  1(1)5T5a4 


+ 


60^672  Q)6T6a4  +  168  a)2T2a4 


27 


27 


(VII. 12) 


.  .V',a> 

'  ' 
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Then, 


G  =  A2D  -  3 ABC  +  2B3 


3012  m6T6a2  +  m  +  m  ^T6a4 


+  23SH  „6t6  .  1536  aj4T4a2  +  1536  ^5^4 
27  27  9 


768  2  2,4 

+  J —  co  t  a 
27 


(VII. 13) 


Then, 


G  =  + 


11,010,048  12.12.4  ,  1,048,576  12.12 

— - - - - cd  t  a  +  — - —  /,v  * 


729 


729 


CD  T 


+  3,932,160  a)8T8&4  _  6, ,291,456  ^12^2 


729 


729 


,  11,010,048  10  10.4  3, 145,728  10J.0.2 

H - - - - -  CD  t  a  — - -  CD  t  a 


729 


729 


+  3^^728  lallTna4 
243 


(VII. 14) 


G2  +  4h3 


-  28,311,552  ^12^4 

729 


(VII. 15) 


n/g2  +  4h3  = 


5320 

27 


icD^T^a2 


(VII. 16) 


■  '\v  >  . . 1 


. '  •  .  ; r:* 


VA'  Si’I  , 
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=  197  lo)6T6a2 


Then  p  =  ±  [  -G  +  n/g2  +  4h3  ] 


(VII. 6) 


+  oA6a2  -  3|*  o>*A4  -  3|*  m6T6a4 


.  512  ffi6T6  +  168  4  4a2  .  768  to5T5a4 

27  27  9 


-  a)2T2a4  +  98.5  io)6T6a2  (VII. 17) 

At  this  point,  it  should  be  noted  that  the  results 
obtained  so  far  are  really  only  consistent  to  a2 
since  if  we  wanted  p  consistent  to  a^  we  would 
have  kept  terms  in  a^  up  to  the  stage  where  we 
calculated  /g2  +  4h3  .  This  arises  because  lower  order 
terms  cancel  exactly.  Therefore,  co nsistent  to  the 
second  order  ( ie . ,  a2) 


P 


±526  6,6a2  .  o)6 t6  +  1|§  *4,4^ 

27  27  27 


+  98.5  i036T6a2 


(VII. 18) 


■  ji  8oY 


■  '  l  /  U  -  ■  '  '  ■":*!  :  T  1  ■  .:>S<  ;.•'/>  JU  - 

-  '  '•  !  ■  !  '  t:v  ;b  1 ::  eor:lr, 

w  ^  Bu<  -  :  ,  "  q  nl  ,  •  -©  .  S  v  : 
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Similar ly. 


=  -i  [  -G  -  n/g2  +  4H3  ] 


(VII. 7) 


+  ^  o)6T6a2  -  512  ,6,6  +  7|  ,4,4,2 


-98.5  ia)6T6a2 


(VII. 19) 


Now  the  general  solution  of  the  modified  cubic  (VII. 2) 
is 

z  =  3/p-  +  3 /5“  (VII.  20) 

3yp"  =  -  co2t2  +  GD2x2a2  +  -^  a2  +  4.6l6  icjo2T2a2 

(VII. 21) 

3/q~  =  -  co2t2  +  <jo2T2a2  +  4  a2  -  4.6l6  io)2x2a2 

D  3  3 

(VII. 22) 

z  =  -  cd2t2  +  co2T2a2  +  a2  (VII. 23) 


is  a  solution  of  the  modified  cubic  equation  in  z. 
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The  complete  solution  is  obtained  by  multiplying  ( VII. 23) 
by  the  cube  roots  of  unity.  If  we  let  R  =  3/p~  and 
S  =  3/q  then  the  three  solutions  are 


Z1 

=  R  +  S 

(VII. 24) 

z2 

=  cdR  +  cd2S 

(VII. 25) 

z3 

=  cd2R  +  odS 

(VII. 26) 

where 

0) 

i 

II 

1  CO 

li 

>— «|CVi 

+ 

i — 1  |CVI 

1 

II 

1  +  0.. 

866  i 

CD2 

=  - 

-1-0 

.866  i 

Then 

(VII. 25) 

and  (VII. 26)  become 

z2 

= 

-  -|(R  +  S)  +  -iyT  i 

(R  -  S) 

(VII. 27) 

z3 

= 

-  1(R  +  S)  -  \JY  i 

(R  -  S) 

(VII. 28) 

•  • 

Z1 

= 

l6  ?  2  .  16  222 

-  — ^  - ^  ^  & 

+  |  a2 

(VII. 23) 

Z2 

= 

8  2  2  8  2  2  2 
+  -^  or  a)  t  a  - 

1  a2  * 

8  o2T2a^ 

or 

Z2 

= 

,8  22  32  222 
+  ^  0)  t  -  ^  cd  t  a 

-  1  a2 

3 

(VII. 29) 

. 


. 

■:  '2.1 
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' 
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Similarly 


z 


3 


,8  2  2 

+  T*" 


+ 


16 

3 


o)^ 


(VII. 30) 


Thus,  we  have  found  the  complete  solution  of  the  modified 
cubic  equation,  consistent  to  the  second  order  in  a. 


Then  the  solution  of  our  original  cubic 


will  be 


16  2  2  ,16  222,8  2 

-  —  a)  t  +  —  a)  t  a  +  a 


[  32  i(DTa2  +  l6  a^  -  16  icuxa^  ]  y. 


+  [  |  cd2t2  -  1  a2  -  22  a,2T2a2  ] 


y- 


1,  1,1,  1  1,  li 

^  icJOT  •  —  2  +  -g  10DT  -  g  -  2  iO)T  - 


B  GOT 


(VII. 31) 

consistent  to  a°.  If  we  keep  only  the  leading  term 
(which  is  the  largest) 


Vl 


1  . 

Tj;  1GDT 


(VII. 32) 


. 
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Similarly, 

y2  =  o  (vii. 33) 

y3  =  -  \  Icdt  (vii. 34) 

y^,  in  our  calculation  presented  here,  is  only  consistent 
to  the  order  shown,  y2  is  zero  because  we  have  only 
calculated  it  to  a0.  For  y^  ,  we  shall  shortly  show 
that  the  first  non-vanishing  term  is  the  a^  one.  We 
could  prove  this  by  carrying  more  terms  in  the  cubic 
but  this  is  extremely  laborious.  To  get  terms  in  a^ 

Q 

for  y2  we  have  to  keep  terms  out  to  a°  in  the  original 
cubic  equation.  This  procedure  was  initially  tried  but 
ended  in  chaos.  The  numerical  work  became  hideous.  So 
abandoning  this  approach,  we  arrived  at  what  y2*s 
first  non- vanishing  term  must  be  by  the  following  sort 
of  argument : 

Consider  only  the  linear  equation  in  y, 

y  [  4  ior^  -  2  icDTa2  -  2  oo2T2a2  -  8  iopT^a2 

4  4  4 

+  higher  orders  in  a  J  -  3  cd  t  a 
+  [  |  cd2t2  -  11  io)3T3  ]  a6  =  0  (VII. 35) 


.  i  >o- 

v  '  o-:.oq,s  :::,i (•;:!  uot:-.o£G£. 
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The  leading  term  of  this  solution  is 

y  =  -  icoT  (VII. 36) 

Then,  consider  the  quadratic  equation  in  y. 

y^  [  +8  co2  t2  -  4  a2  -  32  cD2r2a2  -  6  a^  +  12  icDTa^ 

+  30  cD2T2a^  ]  +  y  [  +  4  icD-^T^  -  2  icDTa2 

-  2  ao2T2a2  -  8  icD^T^a^  +  a^  -  3  icora^ 

+  cu2T2a^  +  17  icD^T^a^  ]  _  3  ^  x^a^  =  0 

(VII. 37) 


This  quadratic  equation  can  be  solved  in  the 
usual  manner,  yielding  the  following  solutions  (keeping 
only  the  leading  term  here) 

yx  =  -  \  ioyr  (VII. 38) 


y2 


icDTa 


4 


(VII. 39) 


. 

o  ^utt  ■'  ..  ■  :s.y 

. 

.  '  :  :■■■  '  :  s.  .  '■:./<> I  c:r)  .> 


• 
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As  we  have  shown  before  (in  solving  the  cubic  directly) 
(VII. 38)  Is  a  solution  of  the  cubic  equation  In  y 
(to  the  order  of  a  specified).  Perhaps  y2,  which 
satisfies  both  the  linear  and  quadratic  equations  in  y3 
satisfies  the  cubic  equation  also.  By  examining  (10.26) 
in  the  light  of  this  proposal,  it  is  immediately  obvious 
that  (VII. 39)  is  the  third  solution  of  the  cubic 
equation.  As  far  as  this  third  solution  is  concerned, 
the  cubic  term  is  not  important.  To  find  the  other 
solution,  i  icjoT  •  — p  ,  we  have  to  solve  the  cubic 
exactly. 

In  conclusion,  the  leading  terms  of  the  three 
solutions  of  the  cubic  equation  are 

y,  =  i  Icdt  •  (VII. 32) 

a 

y2  =  -  -jf  icoxa^  (VII. 36) 

y^  =  -  \  ia>T 


(VII. 34) 


;C  le  ■  ;  1  r  ;  /, 
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APPENDIX  VIII 


Possible  Numerical  Evaluation  of  the  General  Dispersion 

Relation 


From  equation  (8.19)*  the  general  dispersion 
relation  is 


^o 


3 


2  k^X^T^  ,  4  icJDT  k2X2T2  1 

- _  + -  .  - -  j 

(l-ia)T)4  (1-1o3t)  (I-Iojt) 


+  Q 


o 


2 


2  cjd2t2  _  12  icoT  _  k2x2r2  _  8 

(1-1(JDt)2  (l-idOT)2  (l-i<DT)2  (l-i(DT)2 


k2x2T2  _  15  .  kW4  _  1  #  k2x2T2 

(l-icox)2  2  (l-icDT)  (l-ioor)2  (I-Icdt)' 


+  Q„ 


13  io^T  .  k2x2x2  +  2  iopT  +  4  ip2t2 

(1-1cdt)2  (I-Ioot)2  (1-Icot)^  (l-iconr)2 


2  2?  ii  h.  l\. 

k  x  t  +  IT  .  k  xht  _  2  iooT 

(1-Igdt)2  2  (l-iayr)  (1-Icjdt) 

2  2  2  2  2  2 
2 _  <  k  x^t  j  5  iooT  k  x^t 

( 1-  Icjot)  2  ( 1-  io)T)  2  ( 1-  icoT )  2  ( 1-  Icdt)  2 


. 


. .  ...  ...*■ . .  .  — 


. . .  '  V 


v  .  - -  .  .  .... 
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4  4  4 
3  k4x  t 


k2x2T2 


_ _  _  _  .  _  +  2  o2t2 

(  I-Icjdt)  ^  (  1-Icjdt)  2  (I-Icjdt)2  (I-Iot)^ 


=  0 


(VIII. 1) 

Let  £  =  kxT  and  then  multiply  (VIII. 1) 

by  (I-Icdt)^  and  collect  terms 

Q03  [  2 ^  +  4  IcjDT  j£2  +  4  cd2t2  _^2  ] 

+  Qq2  C  2  co2t2  -  4  ico^T-^  -  2  o/^t^  -  12  i cjot  y;2 

-  8  cd2t2  y2  -  -  y2  ] 

+  Q0  [  13  io)T  _^2  -  4  cd2t2  +  4  cd2t2  y2  +  y^ 

+  6  ico^T^  +  2  +  2  y2  ] 

-  5  io)T  y2  -  3  /  -  y2  +  2  od2t2  -  2  ico^T^  =  o 

(VIII. 2) 

Now, 

y  =  (kj  +  ik2)xT 

y2  =  (k-j_2  +  2  ik]_k2  -  k22)x2r2 


. 


■ 


",  «  &  j 


. 
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yp  =  (k-j^  +  3  ik-^kg  -  3  k-^2  -  ik23)  x3x3 
£4  =  (k-j4  +  4  lk13k2  -  6  kx2k22  -  4  ik-^g3 


Then,  (VIII. 2)  becomes 

Qo3  [  2k^4x4x4  +  8  Ikj_3k2X4x4  -  12  k-]_2k2^x4x4 

-  8  ik-|_k23x4x4  +  2  k24x4T4  +  4  icox  *  k^2x2x2 

-  8  aoTk-^x2'!;2  -  4  icox  *  kg2x2x2  +  4  oo2x2  •  k-^x^2 
+  8  ia)2T2k1k2X2T2  -  4  CJD2T2k22X2T2  ] 

+  Qo2  [  2  co2x2  -  4  iao3x3  -  2  co4x4  -  12  iooTk12x2T2 


+  24  o)Tk-j_k2x2r2  +  12  ia)Tk22x2T2  -  8  co2x2k22x2x2 


-  16  i(jQ2x2k]_k2x2x2  +  8  co2x2kg2x2x2  -  k-^4x4x4 


-  30  ik13k2x4T4  +  45  k-,2kQ2x4x4  +  30  ik-,kQ3x4T4 


1  2 


'12 


13  k  4x4x4  -  k-,2x2x2  -  2  ik-.k.o^x2  +  ka2x2x2  ] 


2  2 


12" 


' 
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+  Q0  [  13  lo)xk12x2T2  -  26  0)Tk1k2x2T2  -  13  io)Tk22x2T2 

-  4  co2t2  +  4  o)2T2k12x2x2  +  8  lcD2T2k1k2x2T2 

-  4  (xf  T2k22X2T2  +  kj^X^T^  +  34  Ikj^kgX^T^ 

-  51  k12k22x^T^  -  34  Ikjkr^x^T^  +  k^x^T^ 

+  6  Icd-^t^  +  2  CJD^T^  +  2k-j^2X2T2  +  4  ik^k2X2T2 

-  2  k22x2T2  ]  -  5  la)Tk]_2x2T2  +  10  ooTk^k-px2!;2 
+  5  icoTk22x2T2  -  3k]_^x^T^  -  12  lki^k2X^T^ 

+  18  k^2k22x^r^  +  12  ik1k2^x^T^  -  3  k2^x^x^ 

-  k]_2x2T2  -  2  lk]_k2x2x2  +  k22x2T2  +  2  cjo2t2 

-  2  ia^T3  =  0  (VIII.  3) 

We  can  write  the  integral  Qq  as 

Q0  =  A(k1,k2,a)T)  +  iB(k1,k2,cjDT)  (VIII. 4) 

Q  2  =  A2  +  2  iAB  -  B2 

o 


t 


'p 
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Qq3  =  A3  +  3  iA2B  -  3  AB2  -  IB3 

Using  these  results  and  equating  the  real  and 
imaginary  parts  of  (VIII. 3)  equal  to  zero,  we  have 

(i)  for  the  real  part: 

[A3  -  3  AB2]  [2k14xZ|TZ|  -  12  k-L2k22x4T4  +  2  k^x^T4 

-  8  o)Tk2k2x2T2  +  4  oy^T^ki^x^T^  -  4  o^r2k22x2T2] 
+  [3  A%3  -  B3]  [-  8  k13k2x^x^  +  8  k1k23x^T^ 

-  4  CDTk^2X2T2  +  4  0)Tk22X2T2  -  8  CJO2  T2kj_k2X2  T2  ] 

+  [A2  -  B2]  [2  co2t2  -  2  oo^t^  +  24  ooTk^k2x2T2 

-  8  C02T2k^2X2T2  +  8  (J02T2k22X2T2  -  k-L^x^T4 

+  45  kj2k22x^T^  -  k2^x^T^  -  k-^x2^  +  ^2X2T2  ] 

+  2  AB  [4  o)3t3  +  12  o)Tki2x2T2  -  12  <DTk22x2T2 

+  16  (D2T2kjk2x2i;2  +  30  k^i^x^T^  -  30  kj_k23x^T^ 

+  2  k-^k2x^T^  ] 
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+  A  [  -  26  CDXk^k2X2X2  -  4  002x2  +  4  CD^T^k-j  2X2X2 


-  4  0)2T2k22x2T2  +  -i|  kx4x4x4  -  51  kx2k22x4T4 


17  ,,  4V4^4  ,  o  ...4-4 


+  k^x^x*  +  2  o^x*  +  2  k  2x2x2  -  2  k  *x*x*  ] 
2  2  1  2 


2V2^2 


+  B  [  -  13  a)Tkn2x^T2  +  13  a)Tk0^x2T2  _  g  T2^-^^x2 T2 


-  34  k1^k2x4T4  +  34  k1k2^x4T4  -  6  ao3x3 


-  4  k1k^x2x2  ]  +  10  a>Tk1k2x2T2  -  3  k^4x4x4 


+  18  k^2k^2x4T4  -  3  k£4x4x4  -  k-^2x^T2  +  k^2x2x2 


2  ? 

+  2  oTr  =  0 


(VIII. 5) 


(11)  for  the  Imaginary  part 


[A3  -  3AB2]  [Skj^k^^x4  -  8  k-^kg-^x^T4  +  4  ank^x2^ 


-  4  <jDTk22x2T2  +  8  a)2T2k1k2x2T2  ] 


+  [3A2B  -  B3]  [  2k-)  4x4t4  -  12  k22k22x4T4  +  2  k04x4x4 


-  8  a)xk1k2x2x2  +  4  co2x2k12x2x2  -  4  a)2x2k22x2x2  ] 
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+  [A2  -  B2]  [-  4  CD^T^  -  12  CD  Tk12X2T2  +  12  CDTk22X2T2 

-  1 6  cD2T2k1k2x2T2  -  30  k-^k^x^  ^ 

+  30  k-^k^x^T^  -  2  k-^k2x2T2  ] 

+  2  AB  [2  cd2t2  -  2  cd^t^  +  24  cDTk1k2x2T2  -  8  cD2T2k12x2T2 
+  8  cD2x2k22x2T2  -  •”  k^x^T^  +  45  k-^2k22x^T^ 

-  —  k^x11^  -  k^2x2x2  +  k22x2x2  ] 

+  A  [13  cDTkj2x2T2  -  13  a)Tk22x^T^  +  8  cD2T^k1k2x2r2 
+  34  k^kgX^T^  -  34  k^^x^x^  +  6  cd^t^ 


+  4  k^k2x^T2  ] 
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+  B  [-  26  GDTk^k-pX2!;2  -  4  (JD2T^  +  4 


-  4  CD2T2k22X2T2  +  k^X4”!4  -  51  k12k2^X^T^' 


+  ~  kg^X^T^  +  2  CD^T^  +  2  k12X2T2  -  2  k22X2T2  ] 


-  5  CD  TkX2X2T2  +  5  CDTk22X2T2  -  12  k^k^T4 
+  12  k]_k2^x4T4  -  2  k]_k2x2T2  -  2  cd^t^  =  0 


(VIII. 6) 


The  forms  of  A  and.  B  are  easily  established. 


Q, 


1 


ikxrrr 


+  o°  e-Z^ 

172  L„  Z  -  i/a 


(8.20) 


7 r 


1  r-^o 
172  ^-00 


-7:2 


1  -  i(DT  +  ikXTZ 


dz 


1  f” 


-z2 


77^72  -00  j  _  lcDT  +  i(ic1+ik2)  tzx 


dz 


7 T 


1  f+°°  e-z2 

l/O  f  - 

'  "°°(1-1<2XTZ)  +  ^kjXTZ-CDT) 


dz 


■ 


t 


>i 


■!  ■£ 
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/ 


(l  -  k'2XTz)  -  i(k;LXTZ  -  (jot) 

p  “  “  " 

(l  -  k^XTz)  +  (k^XTZ  -  cot)  ^ 


dz 


(VIII. 7) 


Therefore, 


-z2 


A  = 


1  j.+0°  (l  -  k2XTz)  e' 

t r1/2  -°°  (l  -  k0XTz)2  +  (k,XTZ  -  cot)2 


dz 


(VIII. 8) 


-*2 


B  = 


-1  j+0°  (kiXTZ  -  cot)  e' 

“°°  (l  -  k2XTz)^  +  (k-j_XTZ 


-  cot) 


dz 


(VIII. 9) 


Thus,  equations  VIII. 5  and  VIII. 6  are  two 
rather  complicated  equations  for  k^  =  ^  and  =  a 

in  terms  of  (jot.  For  a  given  value  of  cot.,  one  should, 
in  principle,  he  able  to  solve  these  equations  exactly 
on  the  computer  for  the  corresponding  k^,  k^  . 

Since  Qq  is  simply  connected  to  the  plasma 
dispersion  function  Z (£),  it  is  possible  to  evaluate 
the  terms  A  and  B.  As  pointed  out  earlier.  Fried 
and  Conte^  (1961)  have  tabulated  Z (£)  for  a  limited 


,  '  .  ■'  •  -■■!  ■  c  J  ■  i: ;  ;;  :v  :  qu  t.  a  at  »s0p  sou  i  u 
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range  of  £.  Their  work  is  based  on  the  earlier  developments 

O  nr 

by  K.  A.  Karpov0  (195^-)  who  investigated  the  function 
.  .  _z2  rz  x2 

cd(z)  =  e  J  e  dx  in  the  complex  domain.  The  method 
of  solving  equations  (VIII. 5)  and  (VIII. 6)  might  consist 
of  varying  k^  and  k2  for  a  given  odt  until  both 
equations  are  satisfied.  For  low  frequencies  this 
procedure  should  not  be  too  laborious  because  the  approx¬ 
imate  values  of  k-^  and  k2  are  known  from  experiment, 
but  for  very  high  frequencies  much  more  work  will  be 
involved  because  there  is  no  experimental  data  to 
suggest  the  approximate  values  of  k^  and  k^.  In 
this  latter  case  one  would  have  to  close  in  on  the  values 
of  k^  and  k2  in  a  systematic  way. 


